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Global Solutions of Shock Reflection by 
Large- Angle Wedges for Potential Flow 

By Gui-Qiang Chen and Mikhail Feldman 
Abstract 

When a plane shock hits a wedge head on, it experiences a reflection- 
diffraction process and then a self-similar reflected shock moves outward 
as the original shock moves forward in time. Experimental, computa- 
tional, and asymptotic analysis has shown that various patterns of shock 
reflection may occur, including regular and Mach reflection. However, 
most of the fundamental issues for shock reflection have not been un- 
derstood, including the global structure, stability, and transition of the 
different patterns of shock reflection. Therefore, it is essential to establish 
the global existence and structural stability of solutions of shock reflection 
in order to understand fully the phenomena of shock reflection. On the 
other hand, there has been no rigorous mathematical result on the global 
existence and structural stability of shock reflection, including the case 
of potential flow which is widely used in aerodynamics. Such problems 
involve several challenging difficulties in the analysis of nonlinear par- 
tial differential equations such as mixed equations of elliptic-hyperbolic 
type, free boundary problems, and corner singularity where an elliptic 
degenerate curve meets a free boundary. In this paper we develop a rig- 
orous mathematical approach to overcome these difficulties involved and 
establish a global theory of existence and stability for shock reflection by 
large-angle wedges for potential flow. The techniques and ideas devel- 
oped here will be useful for other nonlinear problems involving similar 
difficulties. 

1. Introduction 

We are concerned with the problems of shock reflection by wedges. These 
problems arise not only in many important physical situations but also are 
fundamental in the mathematical theory of multidimensional conservation laws 
since their solutions are building blocks and asymptotic attractors of general 
solutions to the multidimensional Euler equations for compressible fluids (cf. 
Courant-Friedrichs [16], von Neumann [49], and Glimm-Majda [22]; also see 
[4, 9, 21, 30, 44, 45, 48]). When a plane shock hits a wedge head on, it ex- 
periences a reflection-diffraction process and then a self-similar reflected shock 
moves outward as the original shock moves forward in time. The complexity 



2 



GUI-QIANG CHEN AND MIKHAIL FELDMAN 



of reflection picture was first reported by Ernst Mach [41] in 1878, and ex- 
perimental, computational, and asymptotic analysis has shown that various 
patterns of shock reflection may occur, including regular and Mach reflection 
(cf. [4, 19, 22, 25, 26, 27, 44, 48, 49]). However, most of the fundamental issues 
for shock reflection have not been understood, including the global structure, 
stability, and transition of the different patterns of shock reflection. There- 
fore, it is essential to establish the global existence and structural stability of 
solutions of shock reflection in order to understand fully the phenomena of 
shock reflection. On the other hand, there has been no rigorous mathemati- 
cal result on the global existence and structural stability of shock reflection, 
including the case of potential flow which is widely used in aerodynamics (cf. 
[5, 15, 22, 42, 44]). One of the main reasons is that the problems involve 
several challenging difficulties in the analysis of nonlinear partial differential 
equations such as mixed equations of elliptic-hyperbolic type, free boundary 
problems, and corner singularity where an elliptic degenerate curve meets a 
free boundary. In this paper we develop a rigorous mathematical approach to 
overcome these difficulties involved and establish a global theory of existence 
and stability for shock reflection by large-angle wedges for potential flow. The 
techniques and ideas developed here will be useful for other nonlinear problems 
involving similar difficulties. 

The Euler equations for potential flow consist of the conservation law of 
mass and the Bernoulli law for the density p and velocity potential 

(1.1) d tP + div x (pV x $) = 0, 

(1.2) d t $ + i|V x $| 2 + i(p) = K, 

where K is the Bernoulli constant determined by the incoming flow and/or 
boundary conditions, and 

i'(p) =P(p)/p = c 2 (p)/p 
with c(p) being the sound speed. For polytropic gas, 

p{p)=Kp r , c 2 (p) = /C7/9 7-1 , 7 > 1, k > 0. 
Without loss of generality, we choose k = (7 — so that 

i{p)=p^\ C (p) 2 = (7-1)^-1, 
which can be achieved by the following scaling: 

(x, t, K) — »■ (ax, a 2 t, a~ 2 K), a 2 = kj/(j - 1). 

Equations (1.1)— (1.2) can written as the following nonlinear equation of second 
order: 

(1.3) 8 t p{K - d t <S> - i|V x $| 2 ) + div x (p(K - fy* - i|V x $| 2 )V x $) = 0, 
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where p(s) = s 1 /^ -1 ) = for s > 0. 

When a plane shock in the (x, i)-coordinates, x = (xi, X2) £ R. 2 , with left 
state (p, Vx^) = (pi,«i,0) and right state (po,0, 0),ui > 0,po < Pi-, hits a 
symmetric wedge 

:= {\x2\ < x\ tan 9 W , x\ > 0} 

head on, it experiences a reflection-diffraction process, and the reflection prob- 
lem can be formulated as the following mathematical problem. 

Problem 1 (Initial-Boundary Value Problem). Seek a solution of 
system (1.1) -(1.2) with K = p^~ , the initial condition at t = 0: 



(1.4) (p,$)|t=o 



(po;0) for I X2 1 > x\ tan 6 w ,x\ > 0, 

(pi,-uixi) /or xi < 0, 



anc? i/ie slip boundary condition along the wedge boundary dW: 

(1.5) v$-Hmy = 0, 

where v is the exterior unit normal to dW (see Fig. 1). 
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Figure 1: Initial-boundary value problem 



Notice that the initial-boundary value problem (1.1)— (1.5) is invariant 
under the self-similar scaling: 

(x, t) —> (ax, at), (p, <3?) — > (p, <S>/a) for a / 0. 

Thus, we seek self-similar solutions with the form 

p(x,t) = p(£,7i), §(x,t)=til>(£,rj) for (£, V )=x/t. 
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Then the pseudo-potential function ip = ip — ^(^ 2 + rf ) satisfies the following 
Euler equations for self-similar solutions: 

(1.6) div (p Dip) + 2p = 0, 

(1-7) i|£M 2 + ^ + ^- 1 =^-\ 

where the divergence div and gradient D are with respect to the self-similar 
variables (£,77). This implies that the pseudo-potential function <p(£,rf) is gov- 
erned by the following potential flow equation of second order: 

(1.8) div (p{\D<p\ 2 , <p)D<p) + 2p(\D<p\ 2 ,<p) = 
with 

(1.9) p(\D<p\ 2 ,<p) = M- 1 -<p-l\D<p\ 2 ). 
Then we have 

(1.10) c 2 = cWD^^pl' 1 ) = ( 7 - l)^- 1 - \\D^\ 2 - ip). 

Equation (1.8) is a mixed equation of elliptic-hyperbolic type. It is elliptic 
if and only if 

(1.11) IDtplKcdDtfMpl- 1 ), 
which is equivalent to 

(1.12) \D<p\ < c*(^p ,7) := ^-^^{Pl 1 -^)- 

Shocks are discontinuities in the pseudo- velocity Dip. That is, if £l + and f2 _ := 
Q \ f2+ are two nonempty open subsets of Q C R 2 and S := dQ + n fl is a C 1 
curve where Dip has a jump, then ip £ W^(Q) n C 1 (0 ± U5)n C 2 (0 ± ) is a 
global weak solution of (1.8) in Q if and only if tp is in W^^°(Q,) and satisfies 
equation (1.8) in J7 ± and the Rankine-Hugoniot condition on S: 

(1.13) [p(\Dip\ 2 ,ip)Dip-u] s = 0. 

The continuity of ip is followed by the continuity of the tangential derivative of 
ip across S, which is a direct corollary of irrotationality of the pseudo-velocity. 
The discontinuity S of Dip is called a shock if tp further satisfies the physical 
entropy condition that the corresponding density function p(\Dip\ 2 , ip) increases 
across S in the pseudo-flow direction. We remark that the Rankine-Hugoniot 
condition (1.13) with the continuity of <p across a shock for (1.8) is also fairly 
good approximation to the corresponding Rankine-Hugoniot conditions for the 
full Euler equations for shocks of small strength, since the errors are third-order 
in strength of the shock. 
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The plane incident shock solution in the (x, /j)-coordinates with states 
(/), V x f ) = (po,0,0) and (pi,ui,0) corresponds to a continuous weak solution 
(p of (1.8) in the self-similar coordinates (£,?/) with the following form: 

(1.14) ipo ^ 7] ) = A^ +r] 2 ) for 

(1.15) <Pi&v) = -l(S 2 + V 2 )+ui(Z-So) for £<£ , 
respectively, where 



(1.16) & = w J 2 w;-f ) = -aa- >0 

V Pi-Po Pi- Po 

is the location of the incident shock, uniquely determined by (po, pi, 7) through 
(1.13). Since the problem is symmetric with respect to the axis n = 0, it suffices 
to consider the problem in the half-plane r\ > outside the half-wedge 

A := {£, < 0,77 > 0} U {n > £tan6 w , £ > 0}. 

Then the initial-boundary value problem (1.1)-(1.5) in the (x, ^-coordinates 
can be formulated as the following boundary value problem in the self-similar 
coordinates (£,77). 

Problem 2 (Boundary Value Problem) (see Fig. 2). Seek a solution ip 
of equation (1.8) in the self-similar domain A with the slip boundary condition 
on the wedge boundary dA: 

(1.17) D<p-v\ dA = 
and the asymptotic boundary condition at infinity: 

n 10 s _ jVo for £ > £ ,V > £ tan 0w, , tf2 , 2 

(1.18) <p <p = < , . . n when ^ + i] z 00, 

/or(<^ ,?]>0, 

where (1.18) holds in the sense that lim ||<p — ^||c(A\s K (o)) = 0- 

ii— +00 

Since 931 does not satisfy the slip boundary condition (1.17), the solution 
must differ from 991 in {£ < £0} n A, thus a shock diffraction by the wedge 
occurs. In this paper, we first follow the von Neumann criterion to establish a 
local existence theory of regular shock reflection near the reflection point and 
show that the structure of solution is as in Fig. 3, when the wedge angle is large 
and close to 7r/2, in which the vertical line is the incident shock S = {£ = £0} 
that hits the wedge at the point Po = (£o> £0 tan^), and state (0) and state 
(1) ahead of and behind S are given by <po and <p\ defined in (1.14) and (1.15), 
respectively. The solutions if and ipi differ only in the domain P0P1P2P3 
because of shock diffraction by the wedge vertex, where the curve P0P1P2 is 
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Figure 2: Boundary value problem in the unbounded domain 

the reflected shock with the straight segment PqPi- State (2) behind P$P\ can 
be computed explicitly with the form: 



the constant velocity U2 and the angle S between PqP\ and the £-axis are 
determined by (0 W , po, p±, 7) from the two algebraic equations expressing (1.13) 
and continuous matching of state (1) and state (2) across PqPi-, whose existence 
is exactly guaranteed by the condition on (9 W , po, p±, 7) under which regular 
shock reflection is expected to occur. 

We develop a rigorous mathematical approach to extend the local theory 
to a global theory for solutions of regular shock reflection, which converge to 
the unique solution of the normal shock reflection when W tends to ir/2. The 
solution ip is pseudo-subsonic within the sonic circle for state (2) with center 
(u2, 112 taxi0 w ) and radius C2 > (the sonic speed) and is pseudo-supersonic 
outside this circle containing the arc P\P± in Fig. 3, so that if2 is the unique 
solution in the domain PqP\P±, as argued in [9, 45]. In the domain Q, the 
solution is expected to be pseudo-subsonic, smooth, and C 1 -smoothly matching 
with state (2) across P1P4 and to satisfy ip v = on P2P3; the transonic shock 
curve P1P2 matches up to second-order with PqP\ and is orthogonal to the 
£-axis at the point P2 so that the standard reflection about the £-axis yields 
a global solution in the whole plane. Then the solution of Problem 2 can be 
shown to be the solution of Problem 1. 




W 1 



which satisfies 



Dip ■ v = 



on <9An > 0}; 
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Incident 
shock 




/////// P 3 

Figure 3: Regular reflection 



Main Theorem. There exist 9 C = 9 c (po,pi,j) € (0,7r/2) and a = 

a (Po,Pi,l) G (0,1/2) such that, when 9 W G [0 C , 7r/2), there exists a global 
self-similar solution 

2 



with 



p{-K,t) = { P l~ 1 -$ t -^|V x $| 



x 



for — G A, t > 



of Problem 1 (equivalently, Problem 2) for shock reflection by the wedge, which 
satisfies that, for (£,77) = x/i, 



(1.20) ^ 



9?0 for £ > £0 and r/ > £ tan 

</?i for £ < £0 an d above the reflection shock P0P1P2, 
V? 2 in P0-P1-P4, 



<p is C 1 ' 1 across the part P1P4 of the sonic circle including the endpoints Pi 
and P4, and the reflected shock P0P1P2 is C 2 at Pi and C°° except Pi. More- 
over, the solution tp is stable with respect to the wedge angle in W^(A) and 
converges in W^ c (A) to the solution of the normal reflection described in §3.1 
as 9 W — ► it/2. 

One of the main difficulties for the global existence is that the ellipticity 
condition (1.12) for (1.8) is hard to control, in comparison to our earlier work 
on steady flow [10, 11]. The second difficulty is that the ellipticity degenerates 
at the sonic circle P1P4 (the boundary of the pseudo-subsonic flow). The 
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third difficulty is that, on P1P4, we need to match the solution in Q with ip2 
at least in C 1 , that is, the two conditions on the fixed boundary P1-P4: the 
Dirichlet and conormal conditions, which are generically overdetermined for an 
elliptic equation since the conditions on the other parts of boundary have been 
prescribed. Thus we have to prove that, if ip satisfies (1.8) in 0, the Dirichlet 
continuity condition on the sonic circle, and the appropriate conditions on 
the other parts of dQ, derived from Problem 2, then the normal derivative 
Dtp ■ v automatically matches with Dif2 ■ v along P1P4. We show that, in fact, 
this follows from the structure of elliptic degeneracy of (1.8) on P1-P4 for the 
solution ip. Indeed, equation (1.8), written in terms of the function u = tp — tp2 
in the (x, y)-coordinates defined near P1-P4 such that P1P4 becomes a segment 
on {x = 0}, has the form: 



plus the "small" terms that are controlled by it/ 2 — 6 W in appropriate norms. 
Equation (1.21) is elliptic if u x < 2x/(7 + 1). Thus, we need to obtain the 
C 1 ' 1 estimates near -P1-P4 to ensure \u x \ < 2x/( r y + 1) which in turn implies 
both the ellipticity of the equation in and the match of normal derivatives 
Dip-v = Dip2-v along P1P4. Taking into account the "small" terms to be added 
to equation (1.21), we need to make the stronger estimate \u x \ < 4x/(3(7 + 1)) 
and assume that n/2 — 6 W is appropriately small to control these additional 
terms. Another issue is the non-variational structure and nonlinearity of this 
problem which makes it hard to apply directly the approaches of Caffarelli 
[6] and Alt- Caffarelli- Friedman [1, 2]. Moreover, the elliptic degeneracy and 
geometry of the problem makes it difficult to apply the hodograph transform 
approach in Kinderlehrer-Nirenberg [28] and Chen-Feldman [12] to fix the free 
boundary. 

For these reasons, one of the new ingredients in our approach is to further 
develop the iteration scheme in [10, 11] to a partially modified equation. We 
modify equation (1.8) in 0, by a proper cutoff that depends on the distance 
to the sonic circle, so that the original and modified equations coincide for p 
satisfying \u x \ < 4x/(3(7 + 1)), and the modified equation J\fp> = is elliptic 
in f2 with elliptic degeneracy on P1P4. Then we solve a free boundary problem 
for this modified equation: The free boundary is the curve P1P2, and the free 
boundary conditions on P1P2 are ip = p\ and the Rankine-Hugoniot condition 



On each step, an "iteration free boundary" curve P1P2 is given, and a so- 
lution of the modified equation J\fp = is constructed in 0, with the boundary 
condition (1.13) on P1P2, the Dirichlet condition <p = if2 on the degenerate 
circle P1P4, and Dip ■ v = on P2-P3 an d -P3-P4. Then we prove that p> is 
in fact C 1 ' 1 up to the boundary P1P4, especially \D(ip — 992)] < Cx, by us- 



(1.21) 




in x > and near x = 



(1.13). 
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ing the nonlinear structure of elliptic degeneracy near P1P4 which is modeled 
by equation (1.21) and a scaling technique similar to Daskalopoulos-Hamilton 
[17] and Lin- Wang [40]. Furthermore, we modify the "iteration free bound- 
ary" curve P\Pi by using the Dirichlet condition ip = ip\ on P\Pi- A fixed 
point <p of this iteration procedure is a solution of the free boundary problem 
for the modified equation. Moreover, we prove the precise gradient estimate: 
\u x \ < Ax/ (3(7 + 1)) , which implies that ip satisfies the original equation (1.8). 

Some efforts have been made mathematically for the reflection problem 
via simplified models. One of these models, the unsteady transonic small- 
disturbance (UTSD) equation, was derived and used in Keller-Blank [27], 
Hunter-Keller [26], Hunter [25], Morawetz [44], and the references cited therein 
for asymptotic analysis of shock reflection. Also see Zheng [50] for the pressure 
gradient equation and Canic-Keyfitz-Kim [7] for the UTSD equation and the 
nonlinear wave system. On the other hand, in order to deal with the reflec- 
tion problem, some asymptotic methods have been also developed. Lighthill 
[38, 39] studied shock reflection under the assumption that the wedge angle is 
either very small or close to tt/2. Keller-Blank [27], Hunter-Keller [26], and 
Harabetian [24] considered the problem under the assumption that the shock 
is so weak that its motion can be approximated by an acoustic wave. For a 
weak incident shock and a wedge with small angle in the context of potential 
flow, by taking the jump of the incident shock as a small parameter, the nature 
of the shock reflection pattern was explored in Morawetz [44] by a number of 
different scalings, a study of mixed equations, and matching the asymptotics 
for the different scalings. Also see Chen [14] for a linear approximation of shock 
reflection when the wedge angle is close to tt/2 and Serre [45] for an apriori 
analysis of solutions of shock reflection and related discussions in the context 
of the Euler equations for isentropic and adiabatic fluids. 

The organization of this paper is the following. In §2, we present the po- 
tential flow equation in self-similar coordinates and exhibit some basic prop- 
erties of solutions to the potential flow equation. In §3, we discuss the normal 
reflection solution and then follow the von Neumann criterion to derive the 
necessary condition for the existence of regular reflection and show that the 
shock reflection can be regular locally when the wedge angle is large. In §4, 
the shock reflection problem is reformulated and reduced to a free boundary 
problem for a second-order nonlinear equation of mixed type in a convenient 
form. In §5, we develop an iteration scheme, along with an elliptic cutoff tech- 
nique, to solve the free boundary problem and set up the ten detailed steps of 
the iteration procedure. 

Finally, we complete the remaining steps in our iteration procedure in 
§6-§9: Step 2 for the existence of solutions of the boundary value problem to 
the degenerate elliptic equation via the vanishing viscosity approximation in 
§6; Steps 3-8 for the existence of the iteration map and its fixed point in §7; 
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and Step 9 for the removal of the ellipticity cutoff in the iteration scheme by 
using appropriate comparison functions and deriving careful global estimates 
for some directional derivatives of the solution in §8. We complete the proof 
of Main Theorem in §9. Careful estimates of the solutions to both the "al- 
most tangential derivative" and oblique derivative boundary value problems 
for elliptic equations are made in Appendix, which are applied in §6— §7. 

2. Self-Similar Solutions of the Potential Flow Equation 

In this section we present the potential flow equation in self-similar co- 
ordinates and exhibit some basic properties of solutions of the potential flow 
equation (also see Morawetz [44]). 

2.1. The potential flow equation for self-similar solutions. 

Equation (1.8) is a mixed equation of elliptic- hyperbolic type. It is elliptic if 
and only if (1.12) holds. The hyperbolic-elliptic boundary is the pseudo-sonic 
curve: \Dip\ = c*(</?, p , 7). 

We first define the notion of weak solutions of (1.8)— (1.9). Essentially, we 
require the equation to be satisfied in the distributional sense. 

Definition 2.1 (Weak Solutions). A function p G W^(A) is called a 
weak solution of (1.8) -(1.9) in a self-similar domain A if 

(i) Pq 1 - <P - fl^M 2 > a.e. in A; 

(ii) (p(\Dp\\plp(\Dp\\p)\Dp\) G (L\ oc (A)f; 

(iii) For every ( G C£°(A) ; 

J (p(\Dp\ 2 ,p)Dp ■ D( - 2p(\Dp\ 2 ,p%) c%dr, = 0. 

It is straightforward to verify the equivalence between time-dependent 
self-similar solutions and weak solutions of (1.8) defined in Definition 2.1 in 
the weak sense. It can also be verified that, if ip G C 1,1 (A) (and thus p is twice 
differentiable a.e. in A), then p is a weak solution of (1.8) in A if and only if 
p satisfies equation (1.8) a.e. in A. Finally, it is easy to see that, if A + and 
A~ = A \ A+ are two nonempty open subsets of A C R 2 and S = <9A + n A is a 
C 1 curve where Dp has a jump, then p G W^(D) n C 1 (A ± U S) n C 1 ' 1 (A ± ) 
is a weak solution of (1.8) in A if and only if p is in W^^°(A) and satisfies 
equation (1.8) a.e. in and the Rankine-Hugoniot condition (1.13) on S. 

Note that, for p G C 1 (A ± U S), the condition p G W^(A) implies 



(2.1) 



[p] s = 0. 



GLOBAL SOLUTIONS OF SHOCK REFLECTION BY LARGE-ANGLE WEDGES 11 



Furthermore, the Rankine-Hugoniot conditions imply 

(2.2) - [<Pr,][p<P V ] = on s 

which is a useful identity. 

A discontinuity of Dip satisfying the Rankine-Hugoniot conditions (2.1) 
and (1.13) is called a shock if it satisfies the physical entropy condition: The 
density function p increases across a shock in the pseudo-flow direction. The 
entropy condition indicates that the normal derivative function <p u on a shock 
always decreases across the shock in the pseudo-flow direction. 



2.2. The states with constant density. When the density p is 
constant, (1.8)-(1.9) imply that ip satisfies 

Ap + 2 = 0, ^\D(p\ 2 + (p = const. 

This implies (Aip)^ = 0, (A<p) v = 0, and (tp^ + l) 2 + ip| = 0. Thus, we have 
iptt = -1, tp£ V = 0, (p m = -1, 

which yields 

(2-3) ^,v) = ~\(f + r] 2 ) +ai + bn + c, 

where a, b, and c are constants. 

2.3. Location of the incident shock. Consider state (0): (p ,u ,v ) = 
(po,0, 0) with po > and state (1): (pi,ui,vi) = (pi,ui,0) with p\ > po > 
and ui > 0. The plane incident shock solution with state (0) and state (1) 
corresponds to a continuous weak solution ip of (1.8) in the self-similar coordi- 
nates (£,??) with form (1.14) and (1.15) for state (0) and state (1) respectively, 
where ^ = £o > is the location of the incident shock. 

The unit normal to the shock line is v = (1, 0). Using (2.2), we have 

Pi ~ Po c . n 
ui = Co > 0. 

Pi 



Then (1.9) implies 

7-1 7-1 1 I 7-, |2 1 Pi ~ PO tl 

Pi ~ Po =-~\D<Pi\ -<Pi = ~ —Co- 

z z p 1 

Therefore, we have 



(2.4) ui = (pi-p )\ 



pI-pI 



and the location of the incident shock in the self-similar coordinates is £ 
£o > ui determined by (1.16). 
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3. The von Neumann Criterion and Local Theory for Shock 

Reflection 



In this section, we first discuss the normal reflection solution. Then we 
follow the von Neumann criterion to derive the necessary condition for the 
existence of regular reflection and show that the shock reflection can be regular 
locally when the wedge angle is large, that is, when 9 W is close to ir/2 and, 
equivalently, the angle between the incident shock and the wedge 



(3.1) 

tends to zero. 



a ■= vr/2 - 9 n 



3.1. Normal shock reflection. In this case, the wedge angle is ir/2, 
i.e., a = 0, and the incident shock normally reflects (see Fig. 4). The reflected 
shock is also a plane at £ = £ < 0, which will be defined below. Then U2 = 
V2 = 0, state (1) has form (1.15), and state (2) has the form: 

(3.2) ^ 2 (£,r?) = -i(£ 2 + 7T 2 ) +m(£- £ ) for £ G (£,0), 

where £o = p\U\/(p\ — po) > may be regarded to be the position of the 
incident shock. 





/ 






Reflected 


(2) 




Location of 


shock 






incident shock 


(1) 




\ \ \ \ 


> 



e. Co t 

Figure 1: Normal reflection 

At the reflected shock £ = £ < 0, the Rankine-Hugoniot condition (2.2) 
implies 



(3.3) 



£ = -^^<0. 



92 ~ Pi 



We use the Bernoulli law (1.7): 



P6 



7— l 7— 1 



Pi 1 + t; u i ~ ^i£o = pl 1 + - Co) 
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to obtain 

/o A \ -7-1 7-1 , 1 2 , Pl u l 

(3.4) p] = p? +- Ul + - -. 

* 92 — Pi 

It can be shown that there is a unique solution p2 of (3.4) such that 

92 > Pi- 

Indeed, for fixed 7 > 1 and p±, u\ > and for F{p2) that is the right-hand side 
of (3.4), we have 

lim F(s) = p7 _1 + —Ui > pi 1 , hm F(s) = 00, 

2 

F'(s) = —, PlU \ 2 <0 for a > pi. 

Thus there exists a unique p~2 £ (pi,oo) satisfying pi 1 = F(p2), i.e., (3.4). 
Then the position of the reflected shock £ = £ < is uniquely determined by 
(3-3). 

Moreover, for the sonic speed C2 = \J {l — 1)P2 _1 °f state (2), we have 

(3.5) iei < c 2 . 
This can be seen as follows. First note that 

(3.6) pi' 1 - 9T l = P(p2 -pi), 

where (3 = (7 — l)pl~ 2 > for some p* G (pi,p~2)- We consider two cases, 
respectively. 

Case 1. 7 > 2. Then 

(3.7) 0<( 7 -l)pr 2 </3<(7-l)pr 2 . 
Since /3 > and 92 > 9i, we use (3.4) and (3.6) to find 

P2 = Pi + j£ («i + y it? + 16/3pi) , 

and hence 

4/?Pi 



(3.8) £ - — - 

ui + \ju{ + 16/Spi 

Then using (3.7)-(3.8), p 2 > Pi > 0, and u\ > yields 

4/?Pi 



«i + V u i + i 6 / 5 / 9 ! 

Case 2. 1 < 7 < 2. Then, since 92 > Pi > 0, 

(3.9) 0<( 7 -l)pr 2 </3<(7-l)pr 2 . 
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Since f3 > 0, (3.8) holds by the calculation as in Case 1. Now we use (3.8)-(3.9), 
p2> P\> 0, u\ > 0, and 1 < 7 < 2 to find again 

Id < VPP~i < \/(7 - 1)PT _1 < ^(7 - 1)pV = P2- 
This shows that (3.5) holds in general. 

3.2. The von Neumann criterion and local theory for regu- 
lar reflection. In this subsection, we first follow the von Neumann crite- 
rion to derive the necessary condition for the existence of regular reflection 
and show that, when the wedge angle is large, there exists a unique state (2) 
with two-shock structure at the reflected point, which is close to the solution 
(p~2,U2,V2) = (p2,0,0) of normal reflection for which 9 W = tt/2 in §3.1. 

For a possible two-shock configuration satisfying the corresponding bound- 
ary condition on the wedge r/ = £ tan 9 W , the three state functions ipj , j = 
0,1,2, must be of form (1.14), (1.15), and (1.19) (cf. (2.3)). 

Set the reflected point Pq = (£0, £0 tan 6 W ) and assume that the line that 
coincides with the reflected shock in state (2) will intersect with the axis r] = 
at the point (£, 0) with the angle 6 S between the line and r] = 0. 

Note that <£>i(£, 77) is defined by (1.15). The continuity of ip at (£, 0) yields 

(3.10) W&ri) = ~\{i 2 + ri 2 ) + U2i + V2V+ («i(f-£o) -U2I). 
Furthermore, ip2 must satisfy the slip boundary condition at Po : 

(3.11) V2 = U2t&n8 w . 
Also we have 

(3-12) I = £0-6^. 

tanfc's 

The Bernoulli law (1.7) becomes 

(3.13) pi' 1 = pl~ l + -{u 2 2 + v\) + («i - u 2 )i - u^q. 

Moreover, the continuity of ip on the shock implies that D((p 2 —ipi) is orthogonal 
to the tangent direction of the reflected shock: 

(3.14) (u 2 - ui,v 2 ) • (cos# s ,sin# s ) = 0, 
that is, 

cos^cos^ 

(3.15) u 2 = ui -— . 

cos(^ - e s ) 

The Rankine-Hugoniot condition (1.13) along the reflected shock is 

[pDip] ■ (sin# s , — cos6* s ) = 0, 
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that is, 

(3.16) pi(ui - f)sin0 s = p2{u 2 Sm ^ s „ — - £sin0 s ). 

COS t'ly 

Combining (3.12)-(3.16), we obtain the following system for (p2,0 s ,£): 

(3.17) (i - £ ) cos + £ sin 6> w cot S = 0, 

7-1 , ufcos 2 s ui sin W sin S ~ 7 _i 
(3 - 18) P2 + 2cos^-0 s ) + cos^ - fl.) = ' 

(3.19) (ni cos S tan(0 s — # w ) — ^sin6* s )/92 — Pi(^i — sin0 s = 0. 

The condition for solvability of this system is the necessary condition for the 
existence of regular shock reflection. 

Now we compute the Jacobian J in terms of (p 2 ,0 s ,S t ) at the normal 
reflection solution state (p 2 , |,£) in §3.1 for state (2) when 9 W = it/2 to obtain 

J = -£o((t - i)Pl~ 2 {P2 - pi) - ml) < 0, 

since p2 > pi and £ < 0. Then, by the Implicit Function Theorem, when 9 W is 
near n/2, there exists a unique solution (p2,0 s >£) close to (p~2, §,£) of system 
(3.17)-(3.19). Moreover, (p 2 ,9 s ,0 are smooth functions of a — it/ 2 — 9 W G 
(0, 0i ) for 0i > depending only on po,Pi, and 7. In particular, 

(3.20) \p 2 - p 2 | + k/2 - a | + || - f| + |c 2 - 02 1 < C0, 

where C2 = \J (j — l)pj _1 i s the sonic speed of state (2). 

Reducing 01 > if necessary, we find that, for any G (0, 01), 

(3.21) f < 

from (3.3) and (3.20). Since 9 W G (vr/2 - 0i,tt/2), then S G (7r/4,3vr/4) if 01 
is small, which implies sin^ > 0. We conclude from (3.17), (3.21), and £0 > 
that tan0„, > tan0 s > 0. Thus, 

(3.22) tt/4 < 9 S < 9 W < vr/2. 

Now, given 9 W , we define <p 2 as follows: We have shown that there exists 
a unique solution (p2,0 s ,£) close to (p~2, §,£) of system (3.17)-(3.19). Define 
U2 by (3.15), V2 by (3.11), and <p2 by (3.10). Then the shock connecting state 
(1) with state (2) is the straight line S± 2 = {(£,??) : PiiCv) = ^(C^)}; 
which is £ = r]Cote s + i by (1.15), (3.10), and (3.15). Now (3.19) implies that 
the Rankine-Hugoniot condition (1.13) holds on S12. Moreover, (3.11) and 
(3.15) imply (3.14). Thus the solution (9 S , p 2 ,u 2 ,v 2 ) satisfies (3.11)-(3.19). 
Furthermore, (3.17) implies that the point Po lies on S12, and (3.18) implies 
(3.13) that is the Bernoulli law: 



(3-23) pT-i + h D<p2 \2 + (p2 
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Thus we have established the local existence of the two-shock configuration 
near the reflected point so that, behind the straight reflected shock emanating 
from the reflection point, state (2) is pseudo-supersonic up to the sonic circle 
of state (2). Furthermore, this local structure is stable in the limit 6 W — > 7r/2, 
i.e., a — > 0. 

We also notice from (3.11) and (3.15) with the use of (3.20) and (3.22) 

that 

(3.24) \u 2 \ + H < Co. 

Furthermore, from (3.5) and the continuity of p 2 and | with respect to 9 W on 
(tt/2 — cri,7r/2], it follows that, if a > is small, 

(3.25) ||| < ca. 

In §4-§9, we prove that this local theory for the existence of two shock 
configuration can be extended to a global theory for regular shock reflection. 

4. Reformulation of the Shock Reflection Problem 

We first assume that <p is a solution of the shock reflection problem in the 
elliptic domain f2 in Fig. 3 and that <p — (p 2 is small in C 1 (r2). Under such 
assumptions, we rewrite the equation and boundary conditions for solutions of 
the shock reflection problem in the elliptic region. 

4.1. Shifting coordinates. It is more convenient to change the coordi- 
nates in the self-similar plane by shifting the origin to the center of sonic circle 
of state (2). Thus we define 

(€,V)new ■= (£,V) ~ (u 2 ,V 2 ). 

For simplicity of notations, throughout this paper below, we will always work 
in the new coordinates without changing the notation (£, 77), and we will not 
emphasize this again later. 

In the new shifted coordinates, the domain O is expressed as 

(4.1) n = B C2 (0)n{rj> -v 2 }n{f(rj) <Z<rj cot W }, 

where / is the position function of the free boundary, i.e., the curved part 
of the reflected shock r s /j OC /% := = f(r])}. The function / in (4.1) will be 
determined below so that 

(4.2) ||/-/||<CVt 

in an appropriate norm, specified later. Here £ = l(rj) is the location of the 
reflected shock of state (2) which is a straight line, that is, 

(4.3) l{rj) = rj cot 6 s + i 
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and 

(4.4) i = £-u 2 + v 2 cot e s < 0, 

if a = ir/2 — 9 W > is sufficiently small, since u 2 and v 2 are small and £ < 
by (3.3) in this case. Also note that, since u 2 = v 2 cot 9 W > 0, it follows from 
(3.22) that 

(4.5) £ > £■ 

Another condition on / comes from the fact that the curved part and 
straight part of the reflected shock should match at least up to first-order. 
Denote by P\ = (£1,771) with 771 > the intersection point of the line £ = l(rf) 
and the sonic circle £ 2 + rj 2 = c|, i.e., (£1,7/1) is the unique point for small 
o > satisfying 

(4.6) Kvi) 2 + vi = cl Zi = i(m), vi >o. 

The existence and uniqueness of such point (£1,771) follows from — c 2 < £ < 0, 
which holds from (3.22), (3.25), (4.4), and the smallness of u 2 and v 2 . Then / 
satisfies 

(4.7) /(t?i) = i( m ), f( m ) = Hm) = cot 9 a . 

Note also that, for small a > 0, we obtain from (3.25), (4.4)— (4.5), and l'(rj) = 
cot 9 S > that 

(4.8) -c 2 <£<£<£i<0, C2 _ ||| > ^zJil > 0. 

Furthermore, equations (1.8)-(1.9) and the Rankine-Hugoniot conditions 
(1.13) and (2.1) on T s h oc k do not change under the shift of coordinates. That 
is, we seek (p satisfying (1.8)-(1.9) in 17 so that the equation is elliptic on ip 
and satisfying the following boundary conditions on T s h oc k'- The continuity of 
the pseudo-potential function across the shock: 



(4.9) <p = <pi on T shock 
and the gradient jump condition: 

(4.10) p(\D(p\ 2 , <p)Dip ■ u s = piDtpi ■ u s on T sho ck, 

where u s is the interior unit normal to Q on T s hock- 

The boundary conditions on the other parts of d£l are 

(4.11) p = v?2 on Tsomc = dnn 8B C2 (0), 

(4.12) ip v = on T wedge = dtl n {77 = £tan0 w }, 

(4.13) <p v = on 9fin{j) = -u 2 }. 
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Rewriting the background solutions in the shifted coordinates, we find 

(4.14) <po(S, n) = + v 2 ) - (u 2 t + v 2 v) - \ql 

(4.15) nit, n) = + V 2 ) + («i - u 2 )i - v 2V - -q\ + Ul (u 2 - 

(4.16) v ) = -\{e + r, 2 ) - l -q 2 + ( Ul - u 2 )| + Ul (u 2 - Co), 

where q\ = u 2 , + v\ . 

Furthermore, substituting £ in (4.4) into equation (3.17) and using (3.11) 
and (3.14), we find 

(4.17) ^ = pi( |- (" 1 -"^+"I ), 

which expresses the Rankine-Hugoniot conditions on the reflected shock of 
state (2) in terms of £. We use this equality below. 




■ft (£2, -v 2 ) 



Figure 1: Regular reflection in the new coordinates 

4.2. The equations and boundary conditions in terms of ip = 

if — (f 2 . It is convenient to study the problem in terms of the difference 
between our solution cp and the function ip 2 that is a solution for state (2) 
given by (4.16). Thus we introduce a function 



(4.18) 



ip = ip — (p 2 in Q. 
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Then it follows from (1.8)-(1.10), (3.23), and (4.16) by explicit calculation that 
ip satisfies the following equation in Q: 

(4.19) (c 2 (D^tP,^ V ) - (ipt-0 2 )il>K + (<?{D1>,,J>,Z,r)) - (V>„ -V?)^v 

-2(^-0(^-^ = 0, 

and the expressions of the density and sound speed in O in terms of ip are 

(4.20) p{DiP,iP,^ V ) = (pf l +^ + v ^ - l -\DiP\ 2 -ip)~, 

(4.21) c 2 (L>V,V',e^)=cl + (7-l)(^ + ^-^l^| 2 -V')- 

where p2 is the density of state (2). In the polar coordinates (r, 6) with r = 

\/^ 2 + rj 2 , ip satisfies 

(4.22) 

2 1 1 c 2 1 

{c 2 -{ip r -r) 2 )ip„-^{ip r -r)ipeip r e+^{c 2 -^ipl)ipee+— ip r +^(ip r -2r)ifj 2 = 

with 

(4.23) c 2 = (7 - 1) (pV -^ + r^ r - 1 - (W 2 + ^l) ) ■ 
Also, from (4.11)-(4.12) and (4.16)-(4.18), we obtain 

(4.24) V = on T sonic = dnn dB C2 (0), 

(4.25) ip v = on T wedge = dQ n {r? = £ tan 

(4.26) ^ = -u 2 on dft n {77 = -v 2 }. 

Using (4.15)-(4.16), the Rankine-Hugoniot conditions in terms of ip take 
the following form: The continuity of the pseudo-potential function across (4.9) 
is written as 
(4.27) 

i , ~2 ( l2+i( u i- u 2)+ui(u2-Co) = ((u 1 -u 2 )-r]V2--qi+ui(u2-Co) on T shock , 
that is, 

(4.28) c= 1>&v)+V2TI +i 

Ul - U2 

where £ is defined by (4.4); and the gradient jump condition (4.10) is 
(4.29) 

p{Dip, ip) {Dip - (£, r))) ■ v s = p 1 {m -u 2 -€, -v 2 -rj)-v a on T shock , 

where p{Dip, ip) is defined by (4.20) and v s is the interior unit normal to O on 
r 'shock- If \{u2,V2, Dip)\ < Mi/50, the unit normal v s can be expressed as 

D{(pi -<p) (ni -u 2 - ip£, -v 2 - ip v ) 



(4.30) 



\D(<P1-V)\ y/fui -U2- ipz) 2 + (v 2 + ip v ) 2 ' 
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where we have used (4.15)-(4.16) and (4.18) to obtain the last expression. 

Now we rewrite the jump condition (4.29) in a more convenient form for 
ip satisfying (4.9) when a > and HV'llc 1 ^) are sufficiently small. 

We first discuss the smallness assumptions for a > and ||Y'|lc i (n)- By 
(2.4), (3.20), and (3.24), it follows that, if a is small depending only on the 
data, then 

(4.31) — <c 2 <— , — <P2<— , ^l + vl<-. 

We also require that ||V'llc 1 (Q) i s sufficiently small so that, if (4.31) holds, the 
expressions (4.20) and (4.30) are well-defined in Q, and £ defined by the right- 
hand side of (4.28) satisfies |£| < 7c2/5 for rj G (—v 2 , C2), which is the range of 
T] on r shock- Since (4.31) holds and Q C -B C2 (0) by (4.1), it suffices to assume 

(4.32) U\\ cim <min(^__,min(l, C - 2 )|) =: 6*. 

For the rest of this section, we assume that (4.31) and (4.32) hold. 
Under these conditions, we can substitute the right-hand side of (4.30) for 
v s into (4.29). Thus, we rewrite (4.29) as 

(4.33) F(Dip,ip,U2,v 2 ,£,ri) = on T shock , 
where, denoting p = (pi,P2) G R 2 and z G R, 

(4.34) F(p, z, u 2 ,v 2 , C, V) = {Pip - (£, rj)) ~ Pi {ui ~ u 2 - £, -v 2 - rj) ) ■ v 
with p := p(p, z, £, rj) and v := u(p, u 2 ,v 2 ) defined by 

(4.35) pip, z, £, r?) = (p^ + fri + W2 - ^ - z\ ^ , 

fAQR\ "( \ iu\-U 2 -pi,-V 2 -p 2 ) 

(4.36) u(p, u 2 , v 2 ) = 

V(ui - u 2 -piY + [y 2 +p 2 ) 2 
From the explicit definitions of p and z>, it follows from (4.31) that 



p G C°°(B S ,(0) x (-5*, 5*) x S 2c - 2 (0)), v G C°°(B*.(0) x B Ul/50 (0)), 

where £?r(0) denotes the ball in R 2 with center and radius R and, for k G N 
(the set of nonnegative integers), the C fc -norms of p and v over the regions 
specified above are bounded by the constants depending only on 7,1/1, p 2 ,c 2 , 
and k, that is, by §3, the C fc -norms depend only on the data and k. Thus, 



(4.37) F G C°°(B s .(0) x (-8*, 5*) x B Ul/50 (0) x 5 2c - 2 (0)), 

with its C fc -norm depending only on the data and k. 

Furthermore, since ip satisfies (4.9) and hence (4.28), we can substitute 
the right-hand side of (4.28) for £ into (4.33). Thus we rewrite (4.29) as 

(4.38) V(Dip,if},u 2 ,V2,rj) = onT shock , 
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where 

(4.39) •$(p,z,u 2 ,v 2 ,ri) = F(p, z,u 2 ,v 2 ,(z + v 2 rj)/(u 1 - u 2 ) +£,r]). 

If rj € (-6c 2 /5,6c 2 /5) and \z\ < 5*, then, from (4.8) and (4.31)-(4.32), it 
follows that \(z + v 2 rf)/(u\ — u 2 ) + £| < 7c 2 /5. That is, ((z + v 2 r])/(ui — u 2 ) + 
i, V) G B 252 (0) if r? € (-6c 2 /5, 6c 2 /5) and |z| < 5*. Thus, from (4.37) and 
(4.39), ^ G C°°(X) with H^H^r^ depending only on the data and k £ N, 
where A = B s *(0) x (-<5*,<5*) x B Ul/50 (0) x (-6c 2 /5, 6c 2 /5). 

Using the explicit expression of ^ given by (4.34)-(4.36) and (4.39), we 
calculate 

*((0,0),0,«2,V2,»7) 

(«1-«2)M / JiZ_ , „.2 (Ml-«2K 



f A »7~i 2 ( U l ~ U 2)? \ 



a/ (mi - u 2 ) 2 + v{ V ~ a/ (mi - u 2 ) 2 

Now, using (4.17), we have 

*((0,0),0,u 2 ,u 2 ,r?) = for any (u 2 ,v 2 ,rj) £ B Ul/50 (0) x (-6c 2 /5, 6c 2 /5). 
Then, denoting po = ^ and A? = ((pi,p 2 ),Po,u 2 ,v 2 ,rj) £ A, we have 

2 2 

(4.40) = ^pil^tf ((0,0), 0,^2,^,7?) + PiPj^-W, 

i=0 i,j=0 

where ffij(A') = /^(l - t)D^ P] ^((tp 1 ,tp 2 ),tp ,u 2 ,v 2 ,r])dt for i,j = 0,1,2. 

Thus, <7jj G C°°(„4) and Hs^ij llc? fe (lA) — H^llc^C?) depending only on the data 
and k £ N. 

Next, denoting p' 2 := /5'(p 2 _1 ) = /> 2 /c2 > 0, we compute from the explicit 
expression of * given by (4.34)-(4.36) and (4.39): 

D (PiZ) vI/((0,0),0,0,0,r ? ) = (p 2 ( C 2-| 2 ), (^-^i-p^)^, ^-^__^1). 

Note that, for i = 0,1,2, 

fi^tf ((0, 0), 0, u 2 , u 2 , r?) = fi^tf ((0, 0), 0, 0, 0, n) + hi{u 2 ,v 2 , V ) 
with || || ck ( Bui /so (Q) x ( _ 6 - 2/5 6 - 2/5) ) < ||*|| c » +a(3) for fceN, and |^(« 2 ,«2, < 

lew 



C(|m 2 | + |v 2 |) with C = ||L> 2 ^|| rf -j v Then we obtain from (4.40) that, for all 



X = (p,z,u 2 ,v 2 ,r]) £ A, 



(4.41) *(X) = p' 2 (c 2 2 - C 2 ) Pl + (^— ^ - p' 2 i) ( VP2 - z ) + E 1 (X).p + E 2 (X)z, 

U\ 

where E x £ C°°(:4;R 2 ) and E 2 £ C°°(A) with 
1 1 Bi 1 1 c k (A) — ll^ / llc fc+2 (34)' * = f;2, A; £ N, 

|£i(p, z,u 2 ,u 2 ,r/)| < C(|p| + \z\ + |m 2 | + \v 2 \) for all (p, z,u 2 ,v 2 ,rj) £ A, 
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for C depending only on ||D 2, I'|| c , ( -^- ) . 

Prom now on, we fix (u 2 , v 2 ) to be equal to the velocity of state (2) obtained 
in §3.2 and write Ei(p,z,rj) for Ei(p, z, u 2 , v 2 , rj). We conclude that, if (4.31) 
holds and ijj € satisfies (4.32), then ^ = ip — ip 2 satisfies (4.9)-(4.10) on 

r shock if an d only if ip satisfies conditions (4.28) on r s ^ oc fc, 
(4.42) 

and the functions Ei(p, z, rj),i = 1, 2, are smooth on 

5,5.(0) x (-6*, 5*) x (-6c 2 /5,6c 2 /5) 
and satisfy that, for all (p,z,rj) G .6,5.(0) x (-5*, 5*) x (— 6c 2 /5, 6c 2 /5), 
(4.43) \Ei(p,z,rf)\<C(\p\ + \z\ + <7) 

and, for all (p,z,rj) G B s ,(0) x (-<5*,<5*) x (-6c 2 /5, 6c 2 /5), 
(4-44) Dj,,,,)^)! < C, 

where we have used (3.24) in the derivation of (4.43) and C depends only on 
the data. 

Denote by vq the unit normal on the reflected shock to the region of state 
(2). Then uq = (sin6> s , — cos6* s ) from the definition of 9 S . We compute 

(4-45) (p 2 (^-e 2 ),(^^-p 2 l)r/)^o 

= p' 2 (4 - | 2 ) smO s - (^^i - p 2 £) V co S 6 s 

>^ 2 (<l-l 2 )>o, 

if 7r/2 — 9 S is small and r] £ Proj v (T s h oc k). From (3.14) and (4.30), we obtain 

ps - ^o||L~(r shocfc ) < C\\ D *P\\c(n)- Thus > if a > and W D ^\\c(n) are sma11 
depending only on the data, then (4.42) is an oblique derivative condition on 

r shock- 

4.3. The equation and boundary conditions near the sonic 
circle. For the shock reflection solution, equation (1.8) is expected to be 
elliptic in the domain fl and degenerate on the sonic circle of state (2) which 
is the curve T son i c = 50 n dB C2 (0). Thus we consider the subdomains: 

u m) n':=nn {(£, rj) : dist((£, r?), V sonic ) < 2e}, 

{ ' ' n" :=nn{(£,ri) : dist ((£,?/), r somc ) > e}, 

where the small constant e > will be chosen later. Obviously, £1' and fl" 
are open subsets of ft, and fl = fl' U fl". Equation (1.8) is expected to be 
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degenerate elliptic in f2' and uniformly elliptic in Q," on the solution of the 
shock reflection problem. 

In order to display the structure of the equation near the sonic circle where 
the ellipticity degenerates, we introduce the new coordinates in Q! which flatten 
T S onic and rewrite equation (1.8) in these new coordinates. Specifically, denot- 
ing (r,9) the polar coordinates in the (£, 77)-plane, i.e., (£,17) = (r cos#, r sin0), 
we consider the coordinates: 

(4.47) x = c 2 -r, y = -0 W on Q' . 

By §3.2, the domain V does not contain the point (£, 77) = (0, 0) if e is small. 
Thus, the change of coordinates (£,f/) — > (x,y) is smooth and smoothly in- 
vertible on 0'. Moreover, it follows from the geometry of domain especially 
from (4.2)-(4.7) that, if a > is small, then, in the (x, y)-coordinates, 

tt' = {(x,y) : < x < 2e, < y < tt + arctan (rj(x) / f (r](x))) - 9 W }, 

where r](x) is the unique solution, close to rji, of the equation if + f(r}) 2 = 

(C2-X) 2 . 

We write the equation for ip in the (x, y)-coordinates. As discussed in 
§4.2, ip satisfies equation (4.22)— (4.23) in the polar coordinates. Thus, in the 
(x, y)-coordinates in Q', the equation for ip is 

(4.48) {2x-{ 1 +l)ij x +0 1 )ij xx +02^ X y+{-+0^ yy -{l+O i )^ x +0^y = 0, 

where 
(4.49) 

O^Dijj, 1>,x)=-— + ^(2x - ip x )ip x - ^— ^ + 1 tf) , 
C2 2c 2 c 2 2(c 2 - x) z y/ 

2 

2 {Dip, ip, x) = -^(ip x + c 2 - x)ip y , 

c 2 (c 2 - x) z 

3 (Zty, i/,, x) = ^ (x(2c 2 - x) - (7 - l)(if> + (C2 - x)^) x + 

Co(Co — xr V I 



7 + 1 
"2(c 2 -x 



/ 2 

O4 V, Z) = (x - ^ (V + (C 2 " X)^ x + ^2 + g )) , 

C2 — X V C2 2 2(C2 — x) 2/ / 

<D 5 (Dip, ip, x) = — — ^(^x + 2c 2 - 2x)ip y . 

C2{c 2 -X)- iy ' 

The terms Ok^Dip, ip, x) are small perturbations of the leading terms of equa- 
tion (4.48) if the function ip is small in an appropriate norm considered below. 
In order to see this, we note the following properties: For any (p, z, x) G 
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R 2 x R x (0,c 2 /2) with \p\ < 1, 

|Oi(p,z,x)| <C(|p| 2 + |z| + |x| 2 ), 

(4.50) |0 3 (p,z,x)| + |0 4 (p,*,x)| < C(|p| + \z\ + |x|), 
|0 2 (p,z,x)| + \0 5 ( P ,z,x)\ < C(\p\ + |x| + l)\p\. 

In particular, dropping the terms O k , k = 1, . . . , 5, from equation (4.48), 
we obtain the transonic small disturbance equation (cf. [44]): 

(4.51) (2x - (7 + l)tp x )tp xx + —4> yy - tp x = 0. 

Now we write the boundary conditions on T son i c , r s ^ oc fc, and T we( [g e in the 
(x, y)-coordinates. Conditions (4.24) and (4.25) become 

(4.52) v = o on r sonic = dnn{x = o}, 

(4.53) ^ = ^ = on r W(i3e = ann{!/ = o}. 

It remains to write condition (4.42) on r s ^ oc fc in the (x, y)-coordinates. 
Expressing ^ and ^ in the polar coordinates (r, 0) and using (4.47), we write 
(4.42) on T shock n {x < 2e} in the form: 
(4.54) 

(-P 2 (c| - I 2 ) cos(y + - (^ - p'^)( C2 - x) sin 2 (y + 9 W )) ^ x 

+ sin(y + 9 W ) (-^(c 2 - £ 2 ) + («i - ^) C os(y + w j) ^ 

~ (^ST 1 ^ + E i( D (x, y )^,^,x,y) ■ D( x>y) ip + E 2 (D {Xty) ip,ifj,x,y)7p = 0, 

where z,x,y),i = 1, 2, are smooth functions of (p, z,x,y) € R 2 x R x R 2 
satisfying 

|.Ej(p,z,x,y)| < C (|p| + |z| +cr) for \p\ + \z\+x < e (ui,p 2 ). 

We now rewrite (4.54). We note first that, in the (£, r/)-coordinates, the 
point Pi = T son i c H r s ^ oc fc has the coordinates (£1,771) defined by (4.6). Using 
(3.20), (3.22), (4.3), and (4.6), we find 

< |||- 161 <CVr. 

In the (x, y)-coordinates, the point P\ is (0,yi), where yi satisfies 

(4.55) c 2 cos(yi + 9 W ) = £i, c 2 sin(yi + 6» w ) = r/i, 

from (4.6) and (4.47). Using this and noting that the leading terms of the 
coefficients of (4.54) near Pi = (0,yi) are the coefficients at (x,y) = (0,yi), 
we rewrite (4.54) as follows: 
(4.56) 

-f£nfr* ~ (p' 2 " 6) ^ " " P 2 £i) V 

+£i(£>( :E ,y)V ; , ^, x, y) ■ D^ip + E 2 {D( Xyy) ip, tp, x, y)V> = on T shock n {x < 2e}, 
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where the terms Ei(p, z,x,y),i = 1, 2, satisfy 

(4.57) \Ei(p,z,x,y)\ < C {\p\ + \z\ + x + \y - yi \ + a) 

for (p, z, x, y) G T := {(p, z, x, y) G R 2 x R x R 2 : |p| + \z\ < eo(ui, p?)} and 

(4.58) \\(D imy) E u Dl^EJh^ < C. 

We note that the left-hand side of (4.56) is obtained by expressing the left- 
hand side of (4.42) on T shock l~l {c2 — r < 2e} in the (x, y)-coordinates. Assume 
e < C2/4. In this case, transformation (4.47) is smooth on {0 < C2 — r < 2e} 
and has nonzero Jacobian. Thus, condition (4.56) is equivalent to (4.42) and 
hence to (4.29) on T s h oc k n {x < 2e} if a > is small so that (4.31) holds and 
if HV'llc 1 ^) ^ s sman depending only on the data such that (4.32) is satisfied. 

5. Iteration Scheme 

In this section, we develop an iteration scheme to solve the free boundary 
problem and set up the detailed steps of the iteration procedure in the shifted 
coordinates. 

5.1. Iteration domains. Fix 9 W < tt/2 close to ir/2. Since our problem 
is a free boundary problem, the elliptic domain 0, of the solution is apriori 
unknown and thus we perform the iteration in a larger domain 

(5.1) V = V 9m := B C2 (0) n {r? > -v 2 } D {l(rj) <£< r/cos^}, 

where l{rj) is defined by (4.3). We will construct a solution with C D. 
Moreover, the reflected shock for this solution coincides with {£ = l(r])} outside 
the sonic circle, which implies dV n dB c . 2 (0) = dU PI dB C2 (0) =: T son i c . Then 
we decompose V similar to (4.46): 

V':=Vn{(Z, V ) : dist((e,r ? ),r sonic ) <2e}, 
1 " J V":=Vn{(£,Ti) : dist((^r ? ),r sonic )> £ /2}. 

The universal constant C > in the estimates of this section depends only on 
the data and is independent on 9 W . 

We will work in the (x, y)~coordinates (4.47) in the domain T> n {02 — r < 
Ko}, where kq G (0, C2) will be determined depending only on the data for the 
sonic speed C2 of state (2) for normal reflection (see §3.1). Now we determine kq 
so that ipi — (f2 in the (x, y)-coordinates satisfies certain bounds independent 
of 6 W in V n {c2 — r < ko} if a = tt/2 — 9 W is small. 

We first consider the case of normal reflection 9 W = tt/2. Then, from 
(1.15) and (3.2) in the (x, y)-coordinates (4.47) with C2 = C2 and 9 W = ir/2, 
we obtain 

¥1 — ¥2 = —ui(c2 — x) siny — u±£, for < x < C2, < y < tt/2. 
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Recall £ < and |£| < c 2 by (3.25). Then, in the region V := {0 < x < 
C2, < y < 7r/2}, we have ipi — ip 2 = only on the line 

y = fo,o( x ) '■= arcsin ( - ~ ) for x G (°> ^2 - |£|)- 

Denote kq := (c 2 — |£|)/2. Then kq G (0, C2) by (3.5) and depends only on 
the data. Now we show that there exists <7o > small, depending only on the 
data, such that, if 6 W G (tt/2 — ao, vr/2), then 

(5.3) C _1 < d x (ipi - ip 2 ),-d y (ipi - (p 2 ) < C 

/o,o(0) /o 1 o(«o)+^/2 1 



on [0, /to] x 

(5.4) - > C _1 > on [0, kq] x [0, 



2 ' 
70,0(0), 



(5.5) ^ - ^ < -C7- 1 < on [0, k ] x { /o '° ( ^ + ^ /2 }, 

where /o ' o(k ° 2 )+ " /2 < vr/2. 

We first prove (5.3)-(5.5) in the case of normal reflection 6 W = 7r/2. We 
compute from the explicit expressions of ipi — p 2 and /o,o given above to obtain 

< arcsin (^L) < /o,o(s) < arcsin (P) < C' 1 < f' Gfi (x) < C 

for x G [0, kq], 

dxifi — ^2) = ^isiny, and d y ((fi — <p 2 ) = —u\(c 2 — x)cosy, which imply 
(5.3). Now, (5.4) is true since £ = — c 2 sin(/o 5 o(0)) and thus tp\ — ip 2 = 
u\ (c 2 sin(/o,o(0)) — (c 2 — x)siny), and (5.5) follows from (5.3) since (ipi — 
¥2)(ko, /o,o(ko)) = and (/o,o(«o) + V 2 )/ 2 ~ /o,o(«o) > C~ l . 

Now let 6 W < tt/2. Then, from (3.14)-(4.16) and (4.47), we have 

ip 1 -(p 2 = -(c 2 - x) sin(y + W - 9 s )yJ (m - u 2 ) 2 + v% - (ui - u 2 )£. 

By §3.2, when W — ► 7r/2, we know that (2x2,^2) — ► (0,0), 6* s — ► 7r/2, £ — ► £, 
and thus, by (4.4), we also have £ — > £. This shows that, if do > is small 
depending only on the data, then, for all 9 W G (tt/2 — do, vr/2), estimates (5.3)- 
(5.5) hold with C that is equal to twice the constant C from the respective 
estimates (5.3)-(5.5) for W = tt/2. 

From (5.3)-(5.5) for 6 W G (tt/2 — cro,7r/2) and since 

vn{c 2 -r< Ko } = {if! > <p 2 } n {0 < x < kq, < y < /o ' o( ^ + 7r/2 }, 
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there exists fo ■= fo,ir/2-6 w £ C°°(R+) such that 

(5.6) Vn{c 2 -r < k } = {0 < x < no, 0<y<f (x)}, 

(5.7) fo(0)=y Pl , C' 1 < %{x) < C on[0,«o], 

(5.8) /o,o(0)/2 < /o(0) < / («o) < (/o,o(«o) + tt/2)/2. 

In fact, the line y = fo(x) is the line £ = l(rj) expressed in the (x,y)- 
coordinates, and thus we obtain explicitly with the use of (3.14) that 

(5.9) f (x) = arcsin ( p mg Q _ Q w + 9 S on [0, ko]. 

(c 2 - X) 

5.2. Holder norms in f2. For the elliptic estimates, we need the Holder 
norms in weighted by the distance to the corners P 2 = r s hocfc n {n = — v 2 } 
and P3 = (—1*2, — V2), and with a "parabolic" scaling near the sonic circle. 

More generally, we consider a subdomain 0, C T> of the form f2 := T>n{£ > 
/(?7)} with / G C^R) and set the subdomains fi' :=finr>' and 0" :=flnD" 
defined by (4.46). Let S C 50" be closed. We now introduce the Holder norms 
in fi" weighted by the distance to S. Denote by X = (£, 77) the points of fi" 
and set 

5 X := dist(X, £), <5 x ,y := mm(5 x , &y) for Q". 
Then, for k G R, a G (0, 1), and m G N, define 

mi£oV= E sup (dr (l/,|+fc ' V«wi) , 

(5.10) mSJo" == E -up fe~' 0) l ^ ( g- y ^ (y)l ) , 

ll u llm,a,fi" - — ll^llrrt.O,^" ^ L u J m ,a,n"' 

where L" 3 = d^d^ 2 , and /? = (/?i,/3 2 ) is a multi-index with f3j G N and 

\P\ = /3i + /3 2 . We denote by C^^„ the space of functions with finite norm 
II . |i(k,E) 

Remark 5.1. 7/m > —A; > 1 and A; is an integer, then any function 
u G c!£'a\i„ is C^l" 1,1 up to S, oui noi necessarily up to S. 

In the equation is degenerate elliptic, for which the Holder norms 
with parabolic scaling are natural. We define the norm llV'llsfa'n' as follows: 
Denoting z = (x, y) and z = (x, y) with x, x G (0, 2s) and 

S^ ar) (z,z) := (\x - x\ 2 + mm(x,x)\y - y\ 2 ) a/2 , 
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then, for u G C 2 (Q') n C 1 ' 1 (0') written in the (x, y)-coordinates (4.47), we 
define 

i (par) - y, su p (^ +i/2_2 i^^wi) , 



^ 12 ft O' 

' ' ' — ' z<=n> 



(5.11) M2 Q a' := 7. SU P ( mm(x,x) ' - — -. — r - J, 

k + l=2 z ^ eQ '' z ^ 2 ^ 5a (z,z) ) 

|L.|i(p ar ) ._ |i 7 ,||(P ar ) 1 r,,i(P ar ) 

\\ u \\2,a,Q' ■— Il u ll2,0,n' \. u i2,a,n'- 

To motivate this definition, especially the parabolic scaling, we consider a 

scaled version of the function u(x, y) in the parabolic rectangles: 

(5.12) 



R 



= j(s,t) : \s - x\ < ^, \t - y\ < ^-} fl Q for z = (x, y) G 9,'. 



Denote Qi := (—1, l) 2 . Then the rescaled rectangle (5.12) is 
(5.13) Q[ z) := {(S,T) € Q, : (x + |s, y + ^T) G fi}. 

Denote by u^ z \S,T) the following function in 



(5.14) «W(5,T) := 4^(* + y + for ^ Q?°. 
Then we have 

n- 1 sun Ih/Wll / x < ll?/ll (par) < (7 sun , x 

2eft'n{:r<3e/2} V^ 1 / 26^' ° i 

where C depends only on the domain Q and is independent of e G (0, kq/2). 

5.3. Iteration set. We consider the wedge angle close to tt/2, that is, 
a = ^ — 9 W > is small which will be chosen below. Set So := dT>D{r] = — ^2}- 
Let e,cr > be the constants from (5.2) and (3.1). Let Mi,M 2 > 1. We define 
K = ]C(a,e,M 1 ,M 2 ) by 
(5.15) 

K : = G C^ a (p)r\C\V) : UfiZl < M i> UtaV'''^ < M 2 *A> in V 



for a G (0, 1/2). Then K, is convex. Also, <f> G /C implies that 

ll^ll C i.i(so ^ M i> UWc^m ^ M ^ 

so that K, is a bounded subset in C 1,a (£>). Thus, /C is a compact and convex 
subset of C l > a l 2 (V). 

We note that the choice of constants Mi, M2 > 1 and e, a > below will 
guarantee the following property: 

(5.16) amax(Mi,M 2 ) + e 1/4 Mi + aM 2 /e 2 < C' 1 
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for some sufficiently large C > 1 depending only on the data. In particular, 

(5.16) implies that a < C^ 1 since max(Mi,M2) > 1, which implies ir/2 — 6 w < 
C~ l from (3.1). Thus, if we choose C large depending only on the data, then 
(4.31) holds. Also, for tp G /C, we have 

\(Dij;,iP)(x,y)\ < M lX 2 + M lX in V, U\\ c ^v») < M 2 a. 

Furthermore, < x < 2e in T>' by (4.47) and (5.2). Now it follows from (5.16) 
that HV'llc 1 < 2/(7. Then (4.32) holds if C is large depending only on the 
data. Thus, in the rest of this paper, we always assume that (4.31) holds and 
that ip G K, implies (4.32). Therefore, (4.29) is equivalent to (4.43)-(4.44) for 
if) G K. 

We also note the following fact. 

Lemma 5.1. There exist C and C depending only on the data such that, 
if a,s > and M\,M 2 > 1 in (5.15) satisfy (5.16), then, for every <f> G KL, 

(5.17) ll^ll^ a ' E0UrsO,,ic) < C(M l£ l ~ a + M 2 a). 

Proof. In this proof, C denotes a universal constant depending only on 
the data. We use definitions (5.10)-(5.11) for the norms. We first show that 

(5.18) ll^llfc 1 ,^' rSOnic) < CM^-\ 

where o~( x , y ) '■= dist((x, y), T son i c ) in (5.10). First we show (5.18) in the (x,y)- 
coordinates. Using (5.6), we have V = {0 < x < 2e, < y < fo(x)} with 
Fsonic = {x = 0, < y < f (x)}, where || /o IU~((o,2 £ )) depends only the data, 
and thus dist((x, y), T sonic ) < Cx in V . Then, since 1 1 <^> 1 1 q, x>' — M\, we obtain 
that, for (x, y) G V, 

\<f>(x,y)\ < Mix 2 < M\E 2 , \D<p{x,y)\ < M x x < M x e, 
S\-^\D 2 <l>{x,y)\ = x l ~ a \D 2 <P{x,y)\ < e 1 ~ a M 1 . 

Furthermore, from (5.16) with C > 16, we obtain e < 1/2. Thus, denoting 
z = (x, y) and z = (x, y) with x, x G (0, 2s), we have 

5<r r) (z, z) := (\x - ~x\ 2 + min(x, x)\y - y\ 2 ) a/2 

< (\ x _ S: f + 2e\y-y\ 2 ) a/2 < \z-~z\ a , 
and mm(8 z ,Sz) = mm(x,x), which implies 

min(, 2 , Ss) \*M-*m\ < C£ ,- a min(X) \D^-D*m\ 

<Ce 1-a Afi. 

Thus we have proved (5.18) in the (x, y)-coordinates. By (4.31) and (5.16), 
we have s < c 2 /50 if C is large depending only on the data. Then the change 
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(£> V) ~^ ( x ) y) in an d its inverse have bounded C 3 -norms in terms of the 
data. Thus, (5.18) holds in the (£, r/)-coordinates. 

Since </> G /C, then H^ll^ 1 ^"' 2 ^ < M 2 a. Thus, in order to complete the 
proof of (5.1), it suffices to estimate {mm(8 z , 8 Z ) D ^r^k } i n the case 
z G V \ V" and z G T>" \ V for 8 Z = dist(z, T sonic U E ). From z G V \ V" 
and z G P" \ P', we obtain < C2 — |z| < e/2 and C2 — \z\ > 2e, which 
implies that \z — z\ > 3e/2. We have C2 — |z| < dist(z, r so „j c ) < C(c2 — \z\), 
where we have used (4.31) and (5.1). Thus, mm(8 z ,8z) < C(c 2 — \z\) < Ce. 
Also we have \D 2 (f)(z)\ < M ± by (5.11). If 5~ z > 8 Z , then 8 Z > e/2 and thus 
\D 2 (f>(z)\ < (e/2y 1+a M 2 a by (5.10). Then we have 

U8z < 5 2 , then dist(z, So) < dist(z, T son i c ), which implies by (4.8) that \z—z\ > 
1/C if e is sufficiently small, depending only on the data. Then \D 2 4>(z)\ < 
5~ 1+a M 2 a and 

mm(8 z , 8 Z ) < C{8 Z M X + 5 z 8~ 1+a M 2 a) < C(eM 1 + M 2 a). 

□ 

5.4. Construction of the iteration scheme and choice of a. In 

this section, for simplicity of notations, the universal constant C depends only 
on the data and may be different at each occurrence. 

By (3.24), it follows that, if a is sufficiently small depending on the data, 

then 

(5.19) q 2 < m/10, 

where q 2 = y 'u 2 + v 2 . Let <j) G /C. From (4.15)-(4.16) and (5.19), it follows 
that 

(5.20) (if! -w- <f>)z(£, rj) > m/2 > in P. 

Since (pi — (p 2 = on {£ = l(rj)} and <j) > in P, we have <f> > (pi — ip 2 on 
= /(jy)} nap, where l(rj) is defined by (4.3). Then there exists / G C 1 ' a (R) 
such that 

(5.21) {0 = ^ _ ^1 n P = {(f4,{v),v) ■ v e (-«2, r? 2 )}. 

It follows that /^(ry) > ^ (77) for all rj G [— v 2 ,r] 2 ) and 

(5.22) := > U(v)} r\ v = {(/)< — if 2 } n p. 

Moreover, <9f}+(0) = r s/locfc U T sonic U r wed9e U S , where 

( , 2 ^ r S hock{<P) ■= {£ = U(v)} n 50 + (<A), r soraic : = dp n d£? C2 (o), 

l " j r u ,ed 9e :=52?n{»7 = etajiM. S o (0) :=^ + (0)n{r/ = -i>2}. 
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We denote by Pj, 1 < j < 4, the corner points of Q + (</)). Specifically, P2 = 
r s / lO cfc(0) H So(0) and P3 = (—112, —V2) are the corners on the symmetry line 
{r] = -v 2 }, and Pi = r so „ ic nr s/locA .(>) and P 4 = T sonic C\T wedge are the corners 
on the sonic circle. Note that, since (f> G /C implies <f> = on Tronic; it follows 
that Pi is the intersection point (£1, 771) of the line £ = and the sonic circle 
£ 2 +t? 2 = c 2 ., where (£1,771) is determined by (4.6). 

We also note that _fo = I for G /C. From G /C and Lemma 5.1 with 
a G (0, 1/2), we obtain the following estimate of on the interval (— t>2,?7i): 

(5.24) 11/, - ± °( m ^ 1/2 + ± £l/4 > 

where the second inequality in (5.24) follows from (5.16) with sufficiently large 
C. 

We also work in the (x, y)-coordinates. Denote n := ko/2. Choosing C 
in (5.16) large depending only on the data, we conclude from (5.3)-(5.5) that, 
for every <fi G /C, there exists a function / = /^£C| ^(Ok) suc ^ *kat 

(5.25) Q + (cf>)n{c2-r<K} = {0<x<K, < y < U(x)}, 

with 
(5.26) 

/ (O) = /o(0) > 0, /;>0on (0,«), ||^-/o||^o?i° }) ^ ^(Mx^^+Maa), 

where we have used Lemma 5.1. More precisely, 

(5.27) 

2 

sup (x k - 2 \D k (U-f )(x)\) 

k=0 xe(0,2e) 

+ sup (mm(xi,x 2 )) a — , £ )<CM 1 , 

Xl ^x 2 e(0,2e) V |X!-X 2 r 7 

II/* - /o||2,a,( £ /2, K ) < CM 2 a. 

Note that, in the (£, r;)-coordmates, the angles 6*p 2 and 0p 3 at the corners 
P2 and P3 of respectively satisfy 

(5.28) \0 Pi -l\<I_ fori = 2, 3. 

Indeed, 9p 3 = ir/2 — 6 W . The estimate for 6p 2 follows from (5.24) with (5.16) 
for large C. 
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We now consider the following problem in the domain f2 + (<^>): 



(5.29) 
(5.30) 




in n + (</>), 



P 2 £) (wPv ~ ^) 



on r sfeocfe (0), 



(5.31) V = on T sonic , 

(5.32) ip u = on T^e^e, 

(5.33) Vr, = -^2 on dn + ((f))n{ri = -v 2 }, 

where = Aij(Dip,£,r}) will be defined below, and equation (5.30) is ob- 
tained from (4.42) by substituting <p into E^i = 1,2, i.e., 



Note that, for <p G K and (£,77) G £>, we have (!)(?!>(£, 77), </>(£, 77), 77) G -6,5.(0) x 
(-5*,(5*)x(-6c 2 /5,6c 2 /5) by (4.31)-(4.32). Thus, the right-hand side of (5.34) 
is well-defined. 

Also, we now fix a in the definition of JC. Note that the angles 9p 2 and 6p 3 
at the corners P 2 and -P3 of Q + (4>) satisfy (5.28). Near these corners, equation 
(5.29) is linear and its ellipticity constants near the corners are uniformly 
bounded in terms of the data. Moreover, the directions in the oblique derivative 
conditions on the arcs meeting at the corner -P3 (resp. P 2 ) are at the angles 
within the range (77r/16, 97r/16), since (5.30) can be written in the form ip£ + 
ei/j^-dip = 0, where \e\ < Co near P 2 from rj(P 2 ) = -v 2 , (3.24), (4.43)-(4.44), 
and (5.16). Then, by [35], there exists ao £ (0, 1) such that, for any a 6 (0, cto), 
the solution of (5.29)-(5.33) is in C 1,a near and up to P 2 and P3 if the arcs 
are in C l,a and the coefficients of the equation and the boundary conditions 
are in the appropriate Holder spaces with exponent a. We use a = uq/2 in 
the definition of K, for ao = ao(97r/16, 1/2), where ao(9o,£) is defined in [35, 
Lemma 1.3]. Note that a G (0, 1/2) since G (0, 1). 

5.5. An elliptic cutoff and the equation for the iteration. In 

this subsection, we fix G /C and define equation (5.29) such that 

(i) It is strictly elliptic inside the domain £l + (<fi) with elliptic degeneracy 
at the sonic circle r sonic — 

dn + ((f,) ndB C2 (0); 

(ii) For a fixed point ip = (f> satisfying an appropriate smallness condition 
of \Dip\, equation (5.29) coincides with the original equation (4.19). 

We define the coefficients Aij of equation (5.29) in the larger domain T>. 
More precisely, we define the coefficients separately in the domains V and V" 
and then combine them. 



(5.34) 



E?(Z,T l )=E i (D<!>($,r,),<l>($,r,),T)). 
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In V", we define the coefficients of (5.29) by substituting <p into the coef- 
ficients of (4.19), i.e., 
(5.35) 

^h(£, n) = c 2 m, (/>, e, ri) - (k - 2 > ^ 2 (e, ri) = <?(D<fi, (/>, e, ??) - (</>„ - r/) 2 , 

where 0, and ^ are evaluated at (£,77). Thus, (5.29) in ^ + (0) HP" is a 
linear equation 

A^iptf + 2A\ 2 ip£ v + A l 22 ip m = in ft + (0) n V" . 

From the definition of T>", it follows that \/ £ 2 + r] 2 < C2 — £ in T>" . Then cal- 
culating explicitly the eigenvalues of matrix (^4|j)i<ij<2 defined by (5.35) and 
using (4.31) yield that there exists C = C(j, C2) such that, if e < min(l, C2)/10 
and ll^llc* 1 < £/C 5 then 
(5.36) 

2 

^M 2 < £ A^fo riWi ^ 4 c 2 2 \fi\ 2 for any (£, V ) G 2?" and /* G R 2 . 

The required smallness of e and || 0|| c" 1 is achieved by choosing sufficiently large 
C in (5.16), since <p £ JC. 

In V, we use (4.48) and substitute (j> into the terms 0\, . . . , O5. However, 
it is essential that we do not substitute into the term (7 + l)ip x of the 
coefficient of ip xx in (4.48), since this nonlinearity allows us to obtain some 
crucial estimates (see Lemma 7.3 and Proposition 8.1). Thus, we make an 
elliptic cutoff of this term. In order to motivate our construction, we note 
that, if 

x 4x 

\Ok\ < -r -, -77 — — r, i) x < —, — — r in V, 

10max(c 2 , 1)(7 + 1) 3(7 + 1) 

then equation (4.48) is strictly elliptic in V'. Thus we want to replace the term 

(7 + l)ip x in the coefficient of ip xx in (4.48) by (7 + l)xCi(-), where Ci(-) is 

x 

a cutoff function. On the other hand, we also need to keep form (5.29) for the 
modified equation in the (£, r/)-coordinates, i.e., the form without lower-order 
terms. This form is used in Lemma 8.1. Thus we perform a cutoff in equation 
(4.19) in the (£, r/)-coordinates such that the modified equation satisfies the 
following two properties: 

(i) Form (5.29) is preserved; 

(ii) When written in the (x, y)-coordinates, the modified equation has the 
main terms as in (4.48) with the cutoff described above and corresponding 
modifications in the terms Oi, . . . , O5 of (4.48). 

Also, since the equations in T>' and T>" will be combined and the specific 
form of the equation is more important in T>', we define our equation in a larger 
domain T>' 4e := V n {c 2 - r < Ae}. 
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We first rewrite equation (4.19) in the form 

h + h + h + h = 0, 



where 



h 
h 
h 



(c 2 (£V,^£,?7)-(£ 2 + r/ 2 ))A^ 



/ 4 :=-^(^(W| 2 ) € + ^(W| 2 ),)- 
Note that, in the polar coordinates, Ii, . . . ,I 4 have the following expressions: 

Ji = (c| - r 2 + ( 7 - l)(r^ r - ^|Z?Vf - ^)) A^, 

^2 = + WV, 



(5.37) Ci(s) 



h = r(\Dip\ 2 ) r = 2rip r ip rr + ^ipetpre ~ ^''Ph 

h = -\{Mm?)r + ^u\m 2 )9) 

with \Dip\ 2 =4> 2 + ^ip 2 e and Aip = VVr + ^V'ee + fYv- 

From this, by (4.47), we see that the dominating terms of (4.48) come only 
from and the term 2rvp r ip rr of 1$, i.e., the remaining terms of I3 and I4 

affect only the terms Oi, . . . , O5 in (4.48). Moreover, the term (7 + l)ip x m the 
coefficient of ifj xx in (4.48) is obtained as the leading term in the sum of the 
coefficient (7 — l)rtp r of ip rr in h and the coefficient 2rtp r of ip rr in I3. Thus we 
modify the terms I\ and I3 by cutting off the ^-component of first derivatives 
in the coefficients of second-order terms as follows. Let (1 € C°°(R) satisfy 

8, if |s| < 4/(3(7 + 1)), 

5sign( S )/[3( 7 + l)], if \s\ > 2/(7 + 1), 

so that 

(5.38) C'i(s) > 0, Ci(-s) = -Ci(s) on R; 

(5.39) Ci («) < on {s > 0}. 

Obviously, such a smooth function £1 G C°°(R) exits. Property (5.39) will 
be used only in Proposition 8.1. Now we note that ip£ = ^ip r — ^ipo and 
ipr) = ^ipr + ^ipe, and define 

h ■= (4 ~ r 2 + (7 - l)r(c 2 - 0Ci( ^ 2 + _^ ) " (7 " l)(^Vf + V>)) A^, 
*s == 2(^(c 2 - r)Ci( ^ 2 + _^ ) " %m ~ rttj) (^« + ^ € „) 
+ 2 g(c 2 - r)Ci( ^ + _^ ) + " ^)) ^ + 
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The modified equation in the domain T>' ie is 
(5.40) h + I 2 + h + h = 0. 

By (5.37), the modified equation (5.40) coincides with the original equation 
(4.19) if 



^ 4(c 2 - r) 



3(7 + 1)' 

i.e., if \if} x \ < 4x/(3(7 + 1)) in the (x, y)-coordinates. Also, equation (5.40) is 
of form (5.29) in the (£, ^-coordinates. 

Now we define (5.29) in V' Ae by substituting cf> into the coefficients of (5.40) 

except for the terms involving £i( — -5 ^). Thus, we obtain an equation of 

r(c 2 - r) 

form (5.29) with the coefficients: 
(5.41) 

Aj.iD^tv) = 4 ~ (7 - l)(r(pa - r)Ci(^0f) + i|^| 2 + <t>) 

+ e) + 2*(f (<* - r)Ci(fifg^) - £(^ - , 
^ £, 17) = 4 ~ (7 - 1) (V(c 2 - r)Ci(?g^) + l\D4>\* + 0) 
-(^ + V 2 ) + 2r/(2( C2 - r)Ci(^g^) + M^v ~ 
Aj 2 (Di,,Z, V ) = -(^ + ft) + 2(^(c 2 - r)Ci(?g^) + ^(ta, - f^)), 

where (f>,(f>£, and 4>ri are evaluated at (£,r)). 

Now we write (5.40) in the (x, y)-coordinates. By calculation, the terms 
I\ and ^3 in the polar coordinates are 

h = (c 2 - r 2 + ( 7 - l)(r(c 2 - r)&(^:) ~ \m\ 2 ~ ^)) A^, 

f 3 = 2r(c 2 - r)Cl( ^ )VVr + \i>e^rB ~ 4^0- 
c 2 — r r z r z 

Thus, equation (5.40) in the (x, y)-coordinates in V' 4e has the form 
(5.42) 

{2x-{ 1 +l)xU—)+Oi)^ x +Oi^ xy +(- + Ot] 4>yy-{l+Ot)A+Oii> y = 0, 



X \C2 



with 0^(p,x,y) defined by 
(5.43) 

otiv^y) = -S + ^(2x 2 Ci(f) - ^ (<t> + , 

Ot{x,y) = d k (D(f)(x,y),4>(x,y),x) for i = 2, 5, 

-(7-l)(0 + (c 2 -xK! (£) + 
O* (p, x, y) = ^ (x - ^ (0 + (c 2 - xKi + 3^)) , 
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where p = (pi,P2), and (D(j>, cj>) are evaluated at (x, y). The estimates in (4.50), 
the definition of the cutoff function d, and <\> G K. with (5.16) imply 

(5.44) \6t(p,x,y)\ <C\xf/\ \6f (x, y)\ < C\x\ for k = 2, . . . , 5, 

for all p G R 2 and (x,y) G P^. Indeed, using that (f> G /C implies H^ll^'p; < 
Mi, we find that, for all p G R 2 and (x, y) £ V = V' 2e , 

\6i(p,x,y)\ <C(M? + l)\x\ 2 <C7|x| 3 / 2 , 

(5.45) |0£(x,y)| < C(l + Mi|x|)Mi|x| 3/2 < C\x\ forfe = 2,5, 
|Of(p,x,y)| < C(\x\ + M 1 2 |x| 2 ) < C\x\ for fc = 3, 4. 

In order to obtain the corresponding estimates in the domain V 4e \ V' 2£ , we 
note that V' 4e \ V 2e C V" . Since 2e < x < Ae in D^ e \ £> 2e and G JC implies 
ll^ll27a,©"' Eo) ^ M 2°"> we find that > for an y PGR 2 and (x,y) G V' Ae \ V' 2e , 
\df(p,x,y)\ < C(l + M^a 2 + M 2 a)e 2 < Ce 2 < C\x\ 2 , 

(5.46) \Ot(x,y)\ <C(l + M 2 a)M 2 a <Ce 2 <C\x\ 2 for k = 2, 5, 

|Of(p,x,y)| < C*(£ + M 2 V + M 2 a) < Ce < C|x| for = 3,4. 

Estimates (5.45)-(5.46) imply (5.44). 

The estimates in (5.44) imply that, if <f> G K. and e is sufficiently small 
depending only on the data (which is guaranteed by (5.16) with sufficiently 
large C), equation (5.42) is nonuniformly elliptic in T>' . First, in the (x,y)- 
coordinates, writing (5.42) as 

a-Ulpxx + 2a\2lpxy + 0,22^yy + 0,\l\) x + a 2 ^ y = 0, 

with a,ij = ciij(Dip, x,y) = aji and a% = ai(Dip , x , y) , and using (4.31), we have 
x 2 2 

d^ 2 < V a ij (p,x,y)mfj, j < — \n\ 2 for any (p,x,y) G R^xP^and/x G R 2 . 

D . . . C2 

1,3=1 

In order to show similar ellipticity in the (£, r/)-coordinates, we note that, by 
(4.31), the change of coordinates (£,?7) to (x, y) in T>' Ae and its inverse have C 1 
norms bounded by a constant depending only on the data if e < c 2 /10. Then 
there exists A > depending only on the data such that, for any (p, £, 77) G 
R 2 x V i£ and fi G R 2 , 

2 

(5.47) \(c 2 -r)\n\ 2 < ^jfaMWj < ^Vl 2 , 

where A 2 j(p, £, 77), i, j = 1,2, are defined by (5.41), and r = \f~£?+~rj 2 ~. 

Next, we combine the equations introduced above by defining the coeffi- 
cients of (5.29) in V as follows. Let C2 G C°°(R) satisfy 

£,(*) = { £ i| * I J and < < 10/e on R. 
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Then we define that, for p G R 2 and (£, rj) G T>, 

(5.48) ^-(p,e,»7) = C 2 (c 2 -r)4(e,r/) + (l - C 2 (c 2 - r))4(p, £, r?). 

Then (5.29) is strictly elliptic in X> and uniformly elliptic in T>" with ellipticity 
constant A > depending only on the data and e. We state this and other 
properties of in the following lemma. 

Lemma 5.2. There exist constants A > 0, C, and C depending only on 
the data such that, if M\, M 2 ,e, and a satisfy (5.16), then, for any <fi G K., the 
coefficients Aij(p,£,rj) defined by (5.48), i,j = 1,2, satisfy 

(i) For any (£, rj) 6 V and p, fx G R 2 , 

2 

A(c 2 — r-) l/^l 2 < ^2 A ij(P,£,v)fJ-iH < A~V| 2 withr = \/i 2 + rf- 

(ii) Aij(p,Z,r)) = Aj^r]) for any (£,t?) G V n {c 2 - r > 4e} and p G R 2 , 
where A|-(£, 77) are defined by (5.35). Moreover, 

A\ G C 1 ' a (Dn{c 2 -r> 4e}) 

Wife ||^ij|li )a (x>n{c 2 -r>4e}) - ^ 

(hi) I Ay I + ID^^Ayl < C for any (£, 77) G V n {0 < c 2 - r < 12e} ana 1 
p G R 2 . 

Proof. Property (i) follows from (5.36) and (5.47)— (5.48). Properties (ii)— 
(hi) follow from the explicit expressions (5.35) and (5.41) with <f> G /C. In 
estimating these expressions in property (hi), we use that |sCi( s )l < C which 
follows from the smoothness of Ci and (5.37). □ 

Also, equation (5.29) coincides with equation (5.42) in the domain V. 
Assume that e < Ko/24, which can be achieved by choosing C large in (5.16). 
Then, in the larger domain D n {c 2 — r < 12e}, equation (5.29) written in 
the (x, y)-coordinates has form (5.42) with the only difference that the term 
a?Ci(if) i n the coefficient of ip xx of (5.42) and in the terms Of, O3, and Of 
given by (5.43) is replaced by 

^(C2(x)Cl(^) + (l-C2(x))Cl(^)). 

V x x J 

From this, we have 

Lemma 5.3. There exist C and C depending only on the data such that 
the following holds. Assume that M\,M 2 ,e, and a satisfy (5.16). Let <f> G /C. 
Then equation (5.29) written in the (x,y) -coordinates in V n {c 2 — r < 12e} 
has the form 

(5.49) A n ^j xx + 2A l2 ^ xy + A 22 ^ yy + A^ + A 2 ^ y = 0, 
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where Aij = Aij(ip x ,x,y), Ai = Ai(ip x ,x,y), and A21 = A\ 2 . Moreover, the 
coefficients Aij(p,x,y) and Ai(p,x,y) with p = (pi,P2) G R 2 satisfy 

(i) For any (x, y) G V n {x < 12e} and p, jj, G R 2 , 

x A . 2 

(5-50) -|/i| 2 < V A ij (p,x,y)p i p j < —\p\ 2 ; 

D . . . C 2 

(ii) For any (x, y) ePfl{i< 12e} and p G R 2 , 

KAj^fe^)^) I + \(A, D (p,x,y)A)\ < C; 

(iii) An, A22, and Ai are independent of P2; 

(iv) A12, A 2 i, and A 2 ore independent of p, and 

|(Ai 2 ,l 2 i,A 2 )(x,y)| < C\x\, \D(A 12 ,A 21 ,A 2 )(x,y)\ < C^ 1 ' 2 . 

The last inequality in Lemma 5.3 (iv) is proved as follows. Note that 

(Ai2,A 2 )(x,y) = (0 2 ,0 5 ){D4>(x,y),(t)(x,y),x), 

where O2 and O5 are given by (4.50). Then, by <p G K. and (5.16), we find that, 
for (x,y) G V, i.e., x G (0,2e), 

\D(A 12 ,A 21 ,A 2 )(x,y)\<C(l + Mie)|D^,(x,y)| + (1 + M 1 )\cf> y (x, y)\ 

< C*(l + M x e)M x x xt2 + (7(1 + M x )M x x zt2 < Cx l/2 ; 

and, for (x,y) £ V n {e < x < 12e} C D", we have dist(x,S ) > c 2 /2 > c 2 /4 
so that 

\D(A 12 ,A 21 ,A 2 )(x,y)\ < C(l + M 2 a)M 2 a < Ce < Cx. 
The next lemma follows directly from both (5.37) and the definition of 

Aij. 

Lemma 5.4. Let ft C V, ip G C 2 (Vl), and ip satisfy equation (5.29) with 
cj) = ip in fl. Assume also that ip, written in the (x,y)- coordinates, satisfies 
\ipx\ < 4x/(3(7 + l)) intl' := ttn{c 2 -r < 4e}. Then ip satisfies (4.19) inQ. 

5.6. The iteration procedure and choice of the constants. 

With the previous analysis, our iteration procedure will consist of the following 
ten steps, in which Steps 2-9 will be carried out in detail in §6— §8 and the main 
theorem is completed in §9. 

Step 1. Fix cf) G /C. This determines the domain f2 + (</>), equation (5.29), 
and condition (5.30) on r shock(<f>) > as described in §5.4— §5.5 above. 
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Step 2. In §6, using the vanishing viscosity approximation of equation 
(5.29) via a uniformly elliptic equation 

M(tp) + 5AV> = for 5 G (0, 1) 

and sending 5 — > 0, we establish the existence of a solution ip G C(O + (0)) n 
C 1 (n+ ((p)\T~)nC 2 (n + (</>)) to problem (5.29)-(5.33). This solution satisfies 

(5.51) < V < Ccr in ft + (0), 

where C depends only on the data. 

Step 3. For every s G (0,c 2 /2), set Q" := 0, + (<p)n{c 2 -r > s}. By Lemma 
5.2, if (5.16) holds with sufficiently large C depending only on the data, then 
equation (5.29) is uniformly elliptic in Q" for every s G (0, C2/2), the ellipticity 
constant depends only on the data and s, and the bounds of coefficients in the 
corresponding Holder norms also depend only on the data and s. Furthermore, 
(5.29) is linear on {c 2 — r > 4e}, which implies that it is also linear near the 
corners P 2 and P3. Then, by the standard elliptic estimates in the interior 
and near the smooth parts of dQ + ((f)) n Q" and using Lieberman's estimates 
[35] for linear equations with the oblique derivative conditions near the corners 
(-u 2 , -v 2 ) and T shock (<j)) n {q = -v 2 }, we have 

(5-52) IMlS^f ^ C'Wdl^llL-dS") + 1^1), 

if 11^11^00(^77) + l u 2| < 1, where the second term in the right-hand side comes 
from the boundary condition (5.33), and the constant C(s) depends only on 
the ellipticity constants, the angles at the corners P 2 = T s hock{4 1 ) H i 7 ! = -^2} 
and P3 = (— u 2 , —v 2 ), the norm of T"shock(4>) m C 1,a , and s, which implies that 
C(s) depends only on the data and s. 

Now, using (5.51) and (3.24), we obtain HV'IIloo^) + \v 2 \ < 1 if a is 

sufficiently small, which is achieved by choosing C in (5.16) sufficiently large. 
Then, from (5.52), we obtain 

(5.53) IMfe^ < C(s)a 

for every s G (0, c 2 /2), where C depends only on the data and s. 

Step 4- Estimates of if) in &'(4>) := Q + (4>) n {c 2 — r < e}. We work in the 
(x, y)-coordinates, and then equation (5.29) is equation (5.42) in Q'. 

Step 4-1- L°° estimates of ip in tl + (4>) HP'. Since <f> G fC, the estimates in 
(5.44) hold for large C in (5.16) depending only on the data. We also rewrite 
the boundary condition (5.30) in the (x, y)-coordinates and obtain (4.56) with 
Ei replaced by Ef(x,y) := Ei(D<t>(x,y),4>(x,y),x,y). Using (f> G /C, (4.57), 
(4.58), and (5.27) with / (O) = / (0) = yi, we obtain 

(5.54) \E?(x,y)\ < C(M l£ + M 2 a) < C/C, i = 1, 2, 
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for (x,y) £ ^ s hock{4>) H {0 < x < 2e}. Then, if C in (5.16) is large, we find 
that the function 

is a supersolution of equation (5.42) in fi'(</>) with the boundary condition 
(5.30) on T shock (4>) n {0 < x < 2e}. That is, the right-hand sides of (5.30) 
and (5.42) are negative on w(x,y) in the domains given above. Also, w(x,y) 
satisfies the boundary conditions (5.31)-(5.32) within f2' (</>). Thus, 

(5.55) 0<V(x,y)<— — - in fi'^, 

5(7 + 1) 

if w > ip on x = e. By (5.51), w > ip on x = e if 

Cct < e 2 , 

where C is a large constant depending only on the data, i.e., if (5.16) is satisfied 
with large C. The details of the argument of Step 4.1 are in Lemma 7.3. 

Step 4-2. Estimates of the norm IIV'II^*"^,^- We use the parabolic rescal- 
ing in the rectangle R z defined by (5.12) in which Q' is replaced by Q'(4>). 
Note that R z C Q! for every z = (x,y) G £l' ((/>). Thus, ijj satisfies (5.42) in 
R z . For every z G £l' (</>), we define the functions and (f>^ by (5.14) in the 

(z) 

domain Q\ defined by (5.13). Then equation (5.42) for ip yields the following 
equation for ^ z) {S,T) in Q {z) : 

(5-56) ((1 + | )(2 - ( 7 + mi^fi)) + xO^Ws + xO^^i 

+ (- + x0^ z) )4% -{\ + xO¥> z) )4 z) + x 2 of z) ^ z) = 0, 
where the terms 0^ >,Z \S J T,p), k = 1, . . . , 5, satisfy 



Estimate (5.57) follows from the explicit expressions of OjJ' obtained from 
both (5.43) by rescaling and the fact that 

H (Z) \\ — (7T < CM U 

which is true since ll^ll^fi'^) — Now, sm ce every term O^'^ in (5.56) 
is multiplied by x l3k with 0h > 1 and x G (0, e), condition (5.16) (possibly 
after increasing C depending only on the data) implies that equation (5.56) 

(z) 1 

is uniformly elliptic in Q\ and has the C ' a bounds on the coefficients by a 
constant depending only on the data. 
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Now, if the rectangle R z does not intersect dQ + (<p), then = Q±, where 
Q s = (—s,s) 2 for s > 0. Thus, the interior elliptic estimates in Theorem A.l 
in Appendix imply 

(5-58) II^II^^^C 

where C depends only on the data and \\ip( z ) Hl^qT)- From (5.55), we have 

ll^ll^^) < 1/(7 + 1)- 

Therefore, we obtain (5.58) with C depending only on the data. 

Now consider the case when the rectangle R z intersects dQ + (<fi). From its 
definition, R z does not intersect T son i c . Thus, R z intersects either T s h oc k or 
the wedge boundary T we d ge . On these boundaries, we have the homogeneous 
oblique derivative conditions (5.30) and (5.32). In the case when R z intersects 
T we d ge , the rescaled condition (5.32) remains the same form, thus oblique, and 
we use the estimates for the oblique derivative problem in Theorem A. 3 to 
obtain 

(5.59) ll^'U^G 

where C depends only on the data, since the L°° bound of ip^ in Q[ z) follows 
from (5.55). In the case when R z intersects r s ^ oc fc, the obliqueness in the 
rescaled condition (5.30) is of order x 1 / 2 , which is small since x G (0, 2e). 
Thus we use the estimates for the "almost tangential derivative" problem in 
Theorem A. 2 to obtain (5.59). 

Finally, rescaling back, we have 

(5.oo) M £k w < a 

The details of the argument of Step 4.2 are in Lemma 7.4. 

Step 5. In Lemma 7.5, we extend tp from the domain to T> working 

in the (x, y)-coordinates (or, equivalently in the polar coordinates) near the 
sonic line and in the rest of the domain in the (£, r/)-coordinates, by using the 
procedure of [10]. If C is sufficiently large, the extension of ijj satisfies 

(5-61) < C, 

(5.62) IHfc!^ ^ C ^ 

with C depending only on the data in (5.61) and C(e) depending only on the 
data and e in (5.62). This is obtained by using (5.60) and (5.53) with s > 
determined by the data and e, and by using the estimates of the functions 
and in (5.22), (5.26), and (5.27). 

Step 6. We fix C in (5.16) large depending only on the data, so that 
Lemmas 5.2-5.3 hold and the requirements on C stated in Steps 1-5 above are 
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satisfied. Set Mi = max(2C, 1) for the constant C in (5.61) and choose 
(5.63) , ' 



10max((CM 1 ) 4 ,C) 



This choice of e fixes C in (5.62) depending only on the data and C. Now set 
M 2 = max(C, 1) for C from (5.62) and let 

(C' 1 - e - e 1 / 4 M 1 )e 2 

< a < a Q •- 



2(e 2 max(Mi,M 2 ) + M 2 )' 

where do > since e is defined by (5.63). Then (5.16) holds with constant C 
fixed above. 

Note that the constants ao,e,M±, and M 2 depend only on the data and 

C. 

Step 7. With the constants a,e,M±, and M 2 chosen in Step 6, estimates 
(5.61)-(5.62) imply 

Thus, ijj G fC(a, e, Mi, M 2 ). Then the iteration map J : K, —>■ K. is defined. 

Step 8. In Lemma 7.5 and Proposition 7.1, by the argument similar to 
[10] and the fact that K, is a compact and convex subset of C a ' a / 2 (V), we 
show that the iteration map J is continuous, by uniqueness of the solution 
Y> € C^CD) n C 2 {V) of (5.29)-(5.33). Then, by the Schauder Fixed Point 
Theorem, there exists a fixed point ip G /C. This is a solution of the free 
boundary problem. 

5fep 5. Removal of the cutoff: By Lemma 5.4, a fixed point ijj = 4> satisfies 
the original equation (4.19) mO, + (tp) ii\ip x \ < 4x/(3(7+l)) in f} + (V>)n{c 2 -r < 
4e}. We prove this estimate in §8 by choosing C sufficiently large depending 
only on the data. 

Step 10. Since the fixed point ip £ K. of the iteration map J is a solution 
of (5.29)-(5.33) for <\> = -ip, we conclude 



(i) v e c 1 ' a (n+(^))nc 2 ' a (n + ^)); 

(ii) ip = on r sori j c by (5.31), and ip satisfies the original equation (4.19) in 
n+(V0 by Step 9; 

(iii) = on r so „ ic since H^H^^ < Mr, 

(iv) V = <Pi ~ ^2 on T shock (ip) by (5.21)-(5.23) since = ip; 
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(v) The Rankine-Hugoniot gradient jump condition (4.29) holds on T shock (tp) ■ 
Indeed, as we have showed in (iv) above, the function ip = ip + ip2 satisfies 
(4.9) on T s hock{'4')- Since ip G /C, it follows that tp satisfies (4.28). Also, 
ip on T s h oc k(ip) satisfies (5.30) with <p = ip, which is (4.42). Since ip £ K, 
satisfies (4.28) and (4.42), it has been shown in §4.2 that <p satisfies (4.10) 
on T shock (ip), i.e., ip satisfies (4.29). 

Extend the function ip = ip + <f2 from f2 := £l + (ip) to the whole domain A 
by using (1.20) to define ip in A\Q. Denote Ao := {£ > £o}nA, Ai the domain 
with £ < £o and above the reflection shock P0P1P2, and A2 := A \ (Ao U Ai). 
Set So := {£ = £0} H A the incident shock and Si := P0P1P2 fl A the reflected 
shock. We show in §9 that Si is a C 2 ~curve. Then we conclude that the 
domains Ao, Ai, and A2 are disjoint, OAq n A = So, dAi fl A = So U Si, and 
9A2 n A = Si. Properties (i)-(v) above and the fact that if> satisfies (4.19) in 
O imply that 



(p satisfies equation (1.8) a.e. in A and the Rankine-Hugoniot condition (1.13) 
on the C 2 -curves So and Si , which intersect only at Po 6 dA and are transversal 
at the intersection point. Using this, Definition 2.1, and the remarks after 
Definition 2.1, we conclude that 99 is a weak solution of Problem 2, thus of 
Problem 1. Note that the solution is obtained for every a G (0, ao], i.e., for 
every 6 W £ [n/2 — ao, 7r/2] by (3.1), and that ao depends only on the data since 
C is fixed in Step 9. 

6. Vanishing Viscosity Approximation and Existence of Solutions of 

Problem (5.29) (5.33) 

In this section we perform Step 2 of the iteration procedure described in 
§5.6. Through this section, we keep <p G K, fixed, denote by V := {Pi , P2, P3, P4} 
the set of the corner points of Q + (cf>), and use a € (0, 1/2) defined in §5.4. 

We regularize equation (5.29) by the vanishing viscosity approximation 
via the uniformly elliptic equations 



(6.1) Ns{ip) ■= {An + 5)if)tf + 2Ai 2 ipz v + (A 22 + 6)if> m = in n+(<f>). 



In the domain Q' in the (x, y)-coordinates defined by (4.47), this equation has 



<p€W£°(A), ^GC' 1 (A i )nC7 1 ' 1 (A l ) fori = 0,1, 2, 



N{ip) + 5Aip = 
That is, we consider the equation 



for 5 £ (0,1). 



the form 



(6.2) (S + 2x - ( 7 + (— ) + Of) ip xx + Of iP, 



xy 
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by using (5.42) and writing the Laplacian operator A in the (x, y)-coordinates, 
which is easily derived from the form of A in the polar coordinates. The terms 
Of in (6.2) are defined by (5.43). 

We now study equation (6.1) in Q + (4>) with the boundary conditions 
(5.30)-(5.33). 

We first note some properties of the boundary condition (5.30). Using 

Lemma 5.1 with a G (0,1/2) and (5.16), we find |M|^ a ' E ° Ur ' onie) < C, 
where C depends only on the data. Then, writing (5.30) as 

(6.3) M(^,v)--=b 1 ^,r l )^ + b 2 ^,r l )i; ri + b 3 ^,r l )i; = on T shock (cf>) 

and using (4.43)-(4.45), we obtain 

(6-4) MwZZwZC for , = 1,2,3, 

where C depends only on the data. 

Furthermore, (ft € IC with (5.16) implies that 

| C i < Mie + M 2 a < e 3/A /C. 



Then, using (4.43)-(4.45) and assuming that C in (5.16) is sufficiently large, 

we have 

(6.5) 

(hit V), h(t V)) ■ "(S,V) > lP' 2 (ci " I 2 ) > for any (£,7?) G r shock (<t>), 

h(t, v) > - i 2 ) > for an y (& v) e r shock {(ft), 

HZ, V) - - P 2 i) < e 3/4 ^ any (£, r,) G T shock ((ft), 

HZ, V) + - p'J) | < e 3/4 for any (£, v ) G V shock {<ft). 

Now we write condition (5.30) in the (x, y)-coordinates on T s h ock ((ft) n V. 
Then we obtain the following condition of the form 

( 6 - 6 ) . . . _ 

M(ip)(x,y) =b 1 (x,y)ift x +b 2 (x,y)ift y + b 3 (x,y)ip = on T shock ((ft) n V, 

where h(x,y) = rtfff + & 2 (£, , b 2 (x, y) = r^ff + & 2 (£,??)g, and 
bs(x,y) = &3 Condition (5.30) is oblique, by the first inequality in (6.5). 

Then, since transformation (4.47) is smooth on {0 < c 2 — r < 2e} and has 
nonzero Jacobian, it follows that (6.6) is oblique, that is, 

(6.7) (b 1 (x,y),b 2 (x,y))-u s (x,y)>C- 1 >0 on T shock {(ft) n V 7 , 

where v s = O s (x,y) is the interior unit normal at (x,y) G T s h oc k{(ft) H P' to 
0(0). 

As we have showed in §4.3, writing the left-hand side of (4.42) in the (x, y)— 
coordinates, we obtain the left-hand side of (4.56). Thus, (6.6) is obtained from 
(4.56) by substituting (ft(x,y) into E\ and E 2 . Also, from (5.27) with /^(O) = 
/o(0) = 2/i, we estimate \y-yi\ = |/^(x)-/^(0)| < CM x e on r s/locfc n{x < 2e}. 
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Then, using (4.56)-(4.58) and £1 < 0, we find that, if C in (5.16) is sufficiently 
large depending only on the data, then 

"i,a,r shock (<p)nT>' 

(g8) Si(x,y)<-i^g<0 for (x,y) £ T shock (4>) nV 7 , 

h(x,y) < -\m{P2 + ^flCil) < for (x,y) G T ahock (<f>) (YD', 
h(x,y) < -\{p'^i\ + e 4f L ) < for (x,y) G T shock (<P) f\V , 

where C depends only on the data. 

Now we state the main existence result for the regularized problem. 

Proposition 6.1. There exist C,C,5q > depending only on the data 
such that, if a, e > and M\,Mi > 1 in (5.15) satisfy (5.16), then, for every 
5 G (0,5q), there exists a unique solution ib G C^ato+ty) of (6.1) and (5.30)- 
(5.33), and this solution satisfies 

(6.9) < V(£, V) < Co for (£, V ) G il + (</>), 

(6.10) W>{x,y)\ <C°-x for (x,y) G Q>, 

where we have used coordinates (4.47) in (6.10). Moreover, for any s G 
(0, C2/4), there exists C(s) > depending only on the data and s, but in- 
dependent of 5 £ (0, So), such that 

where Qf((f>) := ^+(0) n {c 2 - r > s}. 

Proof. Note that equation (6.1) is nonlinear and the boundary conditions 
(5.30)-(5.33) are linear. We find a solution of (5.30)-(5.33) and (6.1) as a fixed 
point of the map 

(6.12) J : C 1 ' a / 2 (fi+(^)) C 1 ' a / 2 (fi+(^)) 



defined as follows: For -0 G C 1,Q! / 2 (f2 + (</>)), we consider the linear elliptic 
equation obtained by substituting ip into the coefficients of equation (6.1): 

(6.13) 011^ + 2012^ + 022^ = in ft + (0), 

where 
(6.14) 

a ij (£,Ti) = A ij (Dij>(S,Ti),£,ri) + 68i j for (£,r?) G i,j = 1,2, 

with <5jj = 1 for i = j and for i 7^ j, i,j = 1,2. We establish below the 
existence of a unique solution ib G ^2a/2f2+(<£) ^° ^ ne nnear elliptic equation 
(6.13) with the boundary conditions (5.30)-(5.33). Then we define J(V>) = Y>. 
We first state some properties of equation (6.13). 
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Lemma 6.1. There exists C > depending only on the data such that, if 
a, e > and M\,M<i > 1 in (5.15) satisfy (5.16), and 5 G (0, 1), then, for any 
tp G C ,1 ' Q / 2 (r2 + ((f>)), equation (6.13) is uniformly elliptic in f2 + (</>). • 

2 

(6.15) 5\n\ 2 < OijiMlHtij < 2A-VI 2 for ^,rj) G ^ + (0), A* G R 2 , 

where A is /rom Lemma 5.2. Moreover, for any s G (0, C2/2), i/te ellipticity 
constants depend only on the data and are independent of 5 in £lf(<fi) = + (^)n 
{c 2 - r > s}: 
(6.16) 

2 

A(c 2 -s)M 2 < X) < 2A~V| 2 At * = G ^ A* G R 2 

Ftirt/iermore, 

(6.17) oij G C7 Q / 2 (lF(0)). 

Proof. Facts (6.15)-(6.16) directly follow from the definition of Ojj and 
both the definition and properties of Aij in §5.5 and Lemma 5.2. 

Since Aij(p, £, 77) are independent of p in ri + (0) n {02 — r > 4e}, it follows 
from (5.35), (5.41), and (/> €. K, that aij G C[~*f^ +{(f))nV ,, C C a (f}+(0) n 2?"). 

To show ctij G C Q/2 (O+(0)), it remains to prove that G C Q / 2 (fi(</>) n £>')• 
To achieve this, we note that the nonlinear terms in the coefficients Ay(p, £, 77) 
are only the terms 

(c 2 -r)Ci — 7 7-). 

r(c 2 - r) 

Since Ci is a bounded and C°°-smooth function on R, and ([ has compact 
support, then there exists C > such that, for any s > 0, q G R, 



(6.18) ad(^) < (sup|Ci(t)|) S , %,)«!(-)) 



< C. 



Then it follows that the function 

r(c 2 - r) 

satisfies \F(p,£,r))\ < ||Ci||l°°(r)(c2 - r) for any (p,£,rf) G R 2 x V , and 
\D(p£^F\ is bounded on compact subsets of R 2 x T>' . From this and if) G 
C 1,a / 2 (n+(^)), we have G C* Q / 2 (O+(0)). □ 

Now we state some properties of equation (6.13) written in the (x,y)- 
coordinates. 

Lemma 6.2. There exist A > and C, C > depending only on the 
data such that, if a,e > and M\,Mi > 1 in (5.15) satisfy (5.16), and 5 G 
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(0,1), then, for any ip G C 1 ' a ^ 2 (0, + ((f))), equation (6.13) written in the (x,y)- 
coordinates has the structure 

(6.19) awipxx + 2ai2^ xy + a 2 2^yy + + a 2 ^ y = in Q, + (4>) r\V' u , 
where = ciij(x,y) and cii = &i(x,y) satisfy 

(6.20) a ij: at G C a /\n+(cf>)nV> A£ ) for i,j = 1, 2, 

and the ellipticity condition 
(6.21) 

2 

<5A|/j| 2 < ^2 a,ij(£,r])niiJ,j < A _1 |^| 2 for any (x,y) G & + {(f)) r\T>' 4e , n G R 2 . 
Moreover, 

12 1 

^ < au(x,y) < 5 + 2x, — < a 22 (x, y) < — , -2 < y) < --, 

(6.22) |(ai2,a2i,a 2 )(x,y)| < C|x|, \D(a 12 , a 21 , a 2 )(x, y)\ < C\x\ 1/2 , 
|«m(^,2/) - Oii(0,y)| < C\(x,y) - {0,y)\ a for i = 1,2, 

/or a// (^^(Oijefi+Wn^. 

Proof. By (4.31), if e < C2/IO, then the change of variables from (£,7/) to 
(x, y) in X>4 £ is smooth and smoothly invertible with Jacobian bounded away 
from zero, where the norms and lower bound of the Jacobian depend only on 
the data. Now (6.21) follows from (6.16). 

Equation (6.13) written in the (x, y)-coordinates can be obtained by sub- 

stituting tp into the term x(i( — ) in the coefficients of equation (6.2). Us- 

x 

ing (6.18), the assertions in (6.20) and (6.22), except the last inequality, 
follow directly from (6.2) with (5.43) and (4.50), <f> G /C with (5.16), and 
4> G C7 1 ' Q / 2 (0+ (</>)). 

Then we prove the last inequality in (6.22). We note that, from (6.2) and 
(5.43), it follows that au(x, y) = Fu(D(p, cp, x, y) + Gu(x)x(i(^-) , where Fa and 
Gu are smooth functions, and <p and ip are evaluated at (x,y). In particular, 
since d(-) is bounded, 6^(0, y) = Fu(D(f)(0,y),(j>(0,y),0,y). Thus, assuming 
x > 0, we use the boundedness of Ci and Gu, smoothness of Fa, and <p G K, 
with Lemma 5.1 to obtain 

\a u (x,y) - a u (0,y)\ 

< \Fu(D<t>(x, y), <P(x, y), x, y) - F«(£ty(0, y), tf>(0, y),0, y)\ 

. \r* ( \/- ,$x{x,y) 
+x\Gu(x)( 1 { ) 

x 

< Cx + C(M 1 e 1 - a + M 2 a)\{x,y) - (0,y)| a < C\{x,y) - (0,y)| Q , 

where the last inequality holds since a G (0, 1/2) and (5.16). If x = 0, the only 
difference is that the first term is dropped in the estimates. □ 
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Lemma 6.3 (Comparison Principle). There exists C > depending only 
on the data such that, if a,e > and Mi,M2 > 1 in (5.15) satisfy (5.16), 
and 5 £ (0, 1), the following comparison principle holds: Let if) 6 C(£l + ((ft)) PI 
C 1 (^+(0) \ T~ c ) n C 2 (n+(<p)), let the left-hand sides of (6.13), (5.30), and 

(5.32) -(5.33) are nonpositive for -ip, and let ip > on T son i c . Then 

if> > in n + ((f)). 

Proof. We assume that C is large so that (5.19)-(5.22) hold. 
We first note that the boundary condition (5.30) on T s h oc k (</>) , written as 
(6.3), satisfies 

(6i, b 2 ) ■ v > 0, b 3 < on T shock (4>), 

by (6.5) combined with £ < and />2 > Pi- Thus, if ^ is not a constant in 
f2 + (0), a negative minimum of ip over £l + ((f)) cannot be achieved: 

(i) In the interior of 17 + ((/>), by the strong maximum principle for linear 
elliptic equations; 

(ii) In the relative interiors of shock(<i>) , ^wedge, and dil + ((f) n {rj = —V2}, by 
Hopf's Lemma and the oblique derivative conditions (5.30) and (5.32)- 
(5.33); 

(iii) In the corners P2 and P3, by the result in Lieberman [32, Lemma 2.2], 
via a standard argument as in [20, Theorem 8.19]. Note that we have to 
flatten the curve T s h oc k in order to apply [32, Lemma 2.2] near P2, and 
this flattening can be done by using the C l,a regularity of V shock- 
Using that ^>0on T son i c , we conclude the proof. □ 

Lemma 6.4. There exists C > depending only on the data such that, 
if a,e > and M\,M2 > 1 in (5.15) satisfy (5.16), and 5 £ (0,1), then any 

solution^ g c(n+{(P))nc\n^<p)\T~)riC 2 (n + (<p)) 0/(6.13) and (5.30)- 

(5.33) satisfies (6.9) -(6. 10) with the constant C depending only on the data. 
Proof. First we note that, since ^ + (^) C {q < C2}, the function 

= -V2(v - c 2 ) 

is a nonnegative supersolution of (6.13) and (5.30)-(5.33): Indeed, 

(i) w satisfies (6.13) and (5.33); 

(ii) w is a supersolution of (5.30). This can be seen by using (6.3), (6.5), 
P2 > Pi, ui > 0, p' 2 > £ < 0, and \n\ < C2 to compute on Tg^ock'- 

M{w) = -b2V2-b 3 V2(r]-c 2 ) < -V2(p 2 \i\ + P2 ~ Pl -£ 3/4 (l + 2c 2 ) N ) < 
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if e is small depending on the data, which is achieved by the choice of C 
in (5.16); 

(iii) w is a supersolution of (5.32). This follows from Dw • v = — C2 cos#„, < 
since the interior unit normal on T we( i ge is v = (— sm6 w ,cos9 w ); 

(iv) w > on T sonic . 

Similarly, w = is a subsolution of (6.13) and (5.30)-(5.33). Thus, by the 
Comparison Principle (Lemma 6.3), any solution tp G C(Q + ((f>)) n C 1 (Q + ((/)) \ 
T~) n C 2 (Q + (</>)) satisfies 

< < w(Z,v) for any (£,77) G fl+ (</>). 

Since |^ 2 1 < Ccr, then (6.9) follows. 

To prove (6.10), we work in the (x, y)-coordinates in V n + (</>) and 
assume that C in (5.16) is sufficiently large so that the assertions of Lemma 
6.2 hold. Let v(x, y) = Lax for L > 0. Then 

(i) v is a supersolution of equation (6.19) in Q' n {x < e}: Indeed, the 
left-hand side of (6.19) on v(x,y) = Lax is ai(x,y)La, which is negative in 
2?'nfi+(0) by (6.22); 

(ii) w satisfies the boundary conditions (4.52) on <9S1 + (0) n {x = 0} and 
(4.53) on ^+(0) n{y = 0}; 

(iii) The left-hand side of (6.6) is negative for v on T s h oc k r\ {x < e}: 
Indeed, M(v)(x,y) = La(b\ + b^x) < by (6.8) and since x > in $7'. 

Now, choosing L large so that Le > C where C is the constant in (6.9), 
we have by (6.9) that v > ip on {x = e}. By the Comparison Principle, 
which holds since equation (6.19) is elliptic and condition (6.6) satisfies (6.7) 
and 63 < where the last inequality follows from (6.8), we obtain v > tp in 
n + (4>)n{x < e}. Similarly, -i/t > -v in n + (4>)n{x < e}. Then (6.10) follows. 
□ 

Lemma 6.5. There exists C > depending only on the data such that, if 
a, e > and M\,M2 > 1 in (5.15) satisfy (5.16), and S G (0, 1), any solution 
ip € C{n+{^j) n C^W^P) \ T~) n C 2 (ft+(</>)) of (6.13) and (5.30)-(5.33) 
satisfies 

/or any s £ (0, C2/2), where the constant C(s,ip) depends only on the data, 
ll^llcw2(fi+7>))> an,i s - 

Proo/. From (5.22), (5.24), (6.4)-(6.5), (6.16)-(6.17), and the choice of a 
in §5.4, it follows by [35, Lemma 1.3] that 

(6.24) ||^||M-gur 3 , ocfcW ur_ 9e ) < c( a ^)(||^|| 0(n+W) + \v 2 \) < C(s^)a, 
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where we have used (3.24) and Lemma 6.4 in the second inequality. 

In deriving (6.24), we have used (5.24) and (6.4) only to infer that T^ock ty) 
is a C 1,Q; -curve and 6, G C a (T s hockty)) ■ To improve (6.24) to (6.23), we use 
the higher regularity of T s h oc kty) an d bi, given by (5.24) and (6.4) (and a 
similar regularity for the boundary conditions (5.32)-(5.33), which are given 
on the flat segments and have constant coefficients), combined with rescal- 
ing from the balls B d/2 (z) n £l + ty) for any z G fifty) \ {^2,^3} (with d = 
dist(z, {P2, P3} U So)) into the unit ball and the standard estimates for the 
oblique derivative problems for linear elliptic equations. □ 

Now we show that the solution ip is C 2 '"/ 2 near the corner P4 = T son i c n 
I 'wedgety) ■ We work in V' in the (x, ?/)-coordinates. 

Lemma 6.6. There exists C > depending only on the data such that, if 
a, e > and M±,M 2 > 1 in (5.15) satisfy (5.16), and 5 G (0,1), any solution 

y> g c(n +tyj)nc 1 (n+ty )\T~)nc 2 (n + ty)) 0/(6.13) and (5.30) -(5.33) is 

in C 2 ' a / 2 (B g {P A ) n fl+ty)) for sufficiently small g>0. 

Proof. In this proof, the constant C depends only on the data, 5, and 
||(aij,Oi)|| Ca/ 2(n+(0)) for i,j = 1,2, i.e., C is independent of q. 

Step 1. We work in the (x, y)-coordinates. Then P4 = (0,0) and fl + ty) fl 
#2 e = {a; > 0, y > 0}) n £ 2e for £> G (0, e). Denote 

B+ := B,(0) n{x> 0}, B++ := B e (0) n{x>0,y> 0}. 

Then tp satisfies equation (6.19) in B 2 ^~ and 

(6.25) 4> = on V sonic n B 2e = 5 2e n {x = 0, y > 0}, 

(6.26) ^ = ^ = on T wedge nB 2g = B 2g n{y = 0,x>0}. 

Rescale if) by 

v(z) = ip{gz) for z = (x,y) G -B^ -1- . 
Then t> G <?(£++) n C7 1 (S++ \ {x = 0}) fl C 2 {B+ + ) satisfies 

(6-27) IMIl°o(b++) = IMIl~(b++)> 

and v is a solution of 

(6.28) a^f^ + 2af 2 v xy + a 2 e 2 v yy + ^\ >v x + a 2 ^v y = in i? 2 f+ ) 

(6.29) v = on (95++ n {x = 0}, 

(6.30) u„ = «„ = on 3.B++ n {y = 0}, 

where 
(6.31) 

$fty,y)=aij{Q x iQy)i af\x,y) = Qai(gx,Qy) for (x,y) G B 2 + , i,j = 1,2. 



GLOBAL SOLUTIONS OF SHOCK REFLECTION BY LARGE-ANGLE WEDGES 51 

Thus, a\j satisfy (6.21) with the unchanged constant A > and, since g < 1, 
(6.32) ||(aW,a(«>)|| Co/2 __ < \\^.^ ) \\ ca/a — m) for ij = i )2 . 

Denote Q := {z G : dist(z, dB^) > 1/50}. The interior estimates 

for the elliptic equation (6.28) imply ||u||c?2,a/2(Q) < (£++)• The local 

estimates for the Dirichlet problem (6.28)-(6.29) imply 

(6-33) W v Wc^(B 1/10 (z)nB+ + ) ~ C W v h~(B++) 

for every z = (x,y) G {x = 0,1/2 < y < 3/2}. The local estimates for the 
oblique derivative problem (6.28) and (6.30) imply (6.33) for every z G {1/2 < 
x < 3/2, y = 0}. Then we have 

(6-34) ll w ll C .../.(^^)^ C 'll«llL-m- 

Step 2. We modify the domain + by mollifying the corner at (0, 1) and 
denote the resulting domain by D ++ . That is, D ++ denotes an open domain 
satisfying 

D ++ C B+ + , D ++ \ B 1/w (0, 1) = B++ \ B 1/10 (0, 1), 

and 

dD ++ n%(0,l) is a C 2 '"/ 2 -curve. 

Then we prove the following fact: For any g G C a / 2 (D ++ ), there exists a 
unique solution w G C 2 ' a ^ 2 (D ++ ) of the problem: 



(6.35) 



a[\ w xx + CL22 w yy + a^w x = g in D ++ , 

10 = on dD ++ n {x = 0, y > 0}, 

Wu = w y = on <9£> ++ n {x > 0, y = 0}, 



w = v on n {x > 0, y > 0}, 

with 

(6.36) IMIcw»(5++) ^ C '(IMIl~(b++) + llf llc«/ 2 (D++))- 

This can be seen as follows. Denote by D + the even extension of D ++ 
from {x,y > 0} into {x > 0}, i.e., 

D + := L> ++ U {(x, 0) : x G (0, 1)} U 

where D + - := {(x,y) : (x, -y) G £>++}. Then B± % C D + C B+ and <9L>+ is 
a C7 2 '°/ 2 ~curve. Extend F = (u, 5, ajf , a^, a^) from 5^+ to ^ by setting 



F(x, ~y) = F(x, y) for (x, y) G BJ 
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Then it follows from (6.29)-(6.30) and (6.34) that, denoting by v the restriction 
of (extended) v to dD + , we have v G C 2 ' a / 2 (dD + ) with 

(6.37) \\v\\c 2 -"i 2 (dD+) < C\\ v \\l°°{b++)- 

Also, the extended g satisfies g G (7 a / 2 (Z) + ) with 11511^0/2(75+-) = 1 1 g \ \ Ca /2 (p++)- 
The extended (a[f , d^, d^) satisfy (6.21) and 

|| (d^, a 22 , ^1^)11 0/2 (5+) = IK^ll^' ^22^' ^1^)11 C a ' 2 (B++) 

2 

< X/ IK^i' II c°/ 2 (n+ (</>))' 

Then, by [20, Theorem 6.8], there exists a unique solution w G C 2 ' Q / 2 (-D + ) of 
the Dirichlet problem 

(6.38) a[fw xx + w yy + af^Wx = g in D + , 

(6.39) u; = u on <9L> + , 
and w satisfies 

(6.40) \\w\\c 2 .«/ 2 (D+) ^ C (\\v\\c 2 >«/ 2 (dD+) + \\9\\c-/'(p+))- 

From the structure of equation (6.38) and the symmetry of the domain and the 
coefficients and right-hand sides obtained by the even extension, it follows that 
w, defined by w(x,y) = w(x, —y) in D + , is also a solution of (6.38)-(6.39). By 
uniqueness for (6.38)-(6.39), we find 

w{x,y) = w(x, —y) in D + . 

Thus, w restricted to D ++ is a solution of (6.35), where we use (6.29) to see 
that w = on &D++ n {x = 0,y > 0}. Moreover, (6.37) and (6.40) imply 
(6.36). The uniqueness of the solution w G C 2 > a / 2 {D++) of (6.35) follows from 
the Comparison Principle (Lemma 6.3). 

Step 3. Now we prove the existence of a solution w G C 2 ' a / 2 (D ++ ) of the 
problem: 

anWxx + 20-12 w xy + ®22 w yy + af^w x + a^Wy = in D ++ , 
(6 41) W = ° n dD++ n {x = 0, y > 0}, 

w u = w y = on dD ++ n {y = 0, x > 0}, 
w = v on dD ++ n {x > 0, y > 0}. 

Moreover, we prove that w satisfies 

(6-42) IMIcw2 ( i5++) < c 'lblL~(s++)- 

We obtain such w as a fixed point of map K : C 2 > a / 2 {D++) -► C 2 ' Q / 2 (L>++) 
defined as follows. Let G C 2 ' Q / 2 (IF+). Define 

(6.43) g = -2a[ s 2 ) W xy -a i 2 s) Wy. 
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By (6.22) and (6.31) with g G (0, 1), we find 

(6-44) Ua^la^^—^Cg 1 / 2 , 

which implies 

g G C a / 2 (D++). 

Then, by the results of Step 2, there exists a unique solution w G C 2 ' a / 2 (£>++) 
of (6.35) with g defined by (6.43). We set K[W] = w. 

Now we prove that, if g > is sufficiently small, the map K is a contraction 
map. Let WW G C 2 ' a / 2 (IF+) and u>« := for i = 1,2. Then w := 

_ y;(2) jg a solution of (6.35) with 

g = -2aW(Wg) - Wg>) - af^j 1 ) - Wj 2 >), 
u = 0. 

Then 5 G C a / 2 {D++) and, by (6.44), 

C 2 >»/ 2 (D++)' 

Since u = satisfies (6.29)-(6.30), we can apply both (6.36) and the results of 
Step 2 to obtain 

<^ w m _ w w\\ c2a/2 — y 

where the last inequality holds if g > is sufficiently small. We fix such g. 
Then the map K has a fixed point w G C 2 ' a / 2 (D ++ ) which is a solution of 
(6.41). 

Step 4- Since v satisfies (6.28)-(6.30), it follows from the uniqueness of 
solutions in C{D++) n C 1 (D++ \{x = 0}) n C 2 (L>++) of problem (6.41) tha t 
w = v in Thus v G C 2 < a / 2 (L>++) so that V G C^/^P^Pi) n fi+ ((£)). 

□ 

Now we prove that the solution ip is C l,a near the corner Pi = T son i c n 
rshocfc(0) if ^ is small. 

Lemma 6.7. T/iere exist C > and <5o G (0, 1) depending only on the data 
such that, if a, e > and M\,Mi > 1 in (5.15) satisfy (5.16), and 5 G (0, Sq), 
then any solution if> G C(»+ (^)) n C^+W \ ^WO n C 2 (ft + (</>)) o/ (6.13) 
and (5.30)-(5.33) is in C 1 '" (B e (Pi ) n fi+ (</>)) n C 2 ' a / 2 (P e (Pi) n fi + (<£)), for 
sufficiently small g > depending only on the data and 5, and satisfies 

(6-45) W\\E^$<C(6,$)°, 

where C depends only on the data, 5, and \\ift\\c 1 ><*/ 2 ( n+(<t>) ) • Moreover, for 5 as 
above, 

(6.46) <C(£)(dMt(a;,Pi)) 1+a /or any x G ft + (0), 
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where C depends only on the data and 6, and is independent of ijj. 

Proof. In Steps 1-3 of this proof below, the positive constants C and 
Li, 1 < i < 4, depend only on the data. 

Step 1. We work in the (x, y)-coordinates. Then the point Pi has the 
coordinates (0,y Pl ) with yp 1 = tt/2 + arctan (|£i|/t/i) — 9 W > 0. From (5.25)- 
(5.26), we have 

n + (4>) n B K {Py) = {x > o, y < U(x)} n B e (i\), 

where U(0) = y Pl , /£(0) > 0, and / > y Pl on R+ by (5.7) and (5.26). 

Step 2. We change the variables in such a way that P\ becomes the origin 
and the second-order part of equation (6.13) at P\ becomes the Laplacian. 
Denote 

(6.47) n = Vo^)K(Aj' 
Then, using (6.22) and xp 1 =0, we have 

(6.48) Vc2<V2 < H < V^crf- 
Now we introduce the variables 

(X,Y) := (x/n,y Pl -y). 

Then, for g = e, we have 

(6.49) Q+(<f>) r\B e = {X > 0, Y > F(X)} n B e , 

where F(X) = y Pl - f^fiX). By (5.26), we have < f'^X) < C for all 
X G [0, 2e] if C is sufficiently large in (5.16) so that 2e < k. With this, we use 
/^(O) = yp 1 and (6.48) to obtain 

(6.50) F(0) = 0, -L x \Tb < F'(X) < for X G [0, g]. 

We now write ip in the (X, y)-coordinates. Introduce the function 

v(X, Y) := ^(x, y) = i/>(jiX, y Pl - Y). 

Since tp satisfies equation (6.6) and the boundary conditions (5.32) and (6.19), 
then v satisfies 

12 1 

(6.51) Av := — ^auvxx ai 2 vxY + a22^yy H — a\vx - a 2 VY = 

in {X > 0, Y > F(X)} n B B , 

(6.52) Bv := -M* - W + hv = on {X > 0, F = F(X)} n 

A 1 

(6.53) u = on {X = 0, Y > 0} n B e , 
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where 

dij(X,Y) = a ij (fxX,y Pl — Y), di{X,Y) = Oi(fiX,y Pl -Y), 
b i (X,Y)=b i (jiX,y Pl -Y). 

In particular, from (6.20), (6.22), and (6.47), we have 

(6.54) aij , d t e C a l 2 ({X>0, Y>F(X)}nB e ) , 

(6.55) a 22 (0, 0) = -^a u (0, 0), fii 2 (0, 0) = a 2 (0, 0) = 0, 

A 4 

(6.56) \au(X,Y) -a^(0,0)| <C\{X,Y)\ a fori = 1,2, 

(6.57) \a 12 (X,Y)\ + \d 21 (X,Y)\ + \a 2 (X,Y)\ <C\X\ 1 ' 2 , \~ ai {X,Y)\ < C. 

From (6.8), there exists L 2 > such that 
(6.58) 

-L 2 1 < k(X, Y) < -L 2 for any (X, Y) G {X > 0, Y = F(X)} n 5 e . 
Moreover, (6.7) implies 

(6.59) {b 1 ,b 2 )-VF>0 on {X >0, Y = F(X)}nB e , 

where vp = ^f(X,Y) is the interior unit normal at (X, Y) G {X > 0, V = 
F(X)} n J3 e . Thus condition (6.52) is oblique. 

Step 3. We use the polar coordinates (r,9) on the (X, Y)-plane, i.e., 

(X, Y) = (r cos 0, r sin 0). 

From (6.50), we have F,F'<0 on (0, g), which implies that (X 2 + F(X) 2 )' > 
on (0,g). Then it follows from (6.50) that, if 5 > is a small constant 
depending only on the data and g is a small constant depending only on the 
data and S, there exist a function Op € C 1 (R + ) and a constant L3 > such 
that 

(6.60) {X > 0, Y > F(X)} n B e = {0 < r < <?, F (r) < < ^/2} 
with 

(6.61) < 9 F {r) < 0. 

Choosing sufficiently small So > 0, we show that, for any S £ (0, So), a 
function 

(6.62) w(r,0) =r 1+a cosG(6), with G(6) = ^-^(0 - j), 

is a positive supersolution of (6.51)-(6.53) in {X > 0, Y > n B r 

By (6.49) and (6.60)-(6.61), we find that, when < S < S < ( 8 (3+")I 3 ) 2 » 

~ + < G(9) < I - ^vr for all (r, 0) G fi+(<£) fl 
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In particular, 
(6.63) 

\ — a 

cos(G(0)) > sin {—^t-k) > for all (r, 9) efl+(^)nB f \{I = F = 0}, 

which implies 

w > in {X > 0, Y > F(X)} n B e . 

By (6.60)-(6.61), we find that, for all r G (0, £>) and <5 G (0,S ) with small 
5o > 0, 

cos(9 F (r)) > 1 - C5 > 0, | sin(0 F (r))| < Cv^o- 

Now, possibly further reducing <5o, we show that w is a supersolution of 
(6.52). Using (6.48), (6.52), (6.58), the above estimates of (9 F ,G(9 F )) derived 
above, and the fact that 9 = 9 F on {X > 0, Y = F(X)} n -Bp, we have 

t Q I 

£w < — r a ((a + 1) cos(e F ) cos(G(9 F )) + — ]p sin(0 F ) sm(G(9 F )j) 

+Cr Q |6 2 | +Cr Q+1 |6 3 | 

L 2 sin(^7r) C 



,_ r . ((1 _ CTo)( ^™_^_ c)< o, 



Cv^O ^2 

if 5o is sufficiently small. We now fix 5q that satisfies all the smallness assump- 
tions made above. 

Finally, we show that w is a supersolution of equation (6.51) in (X, Y) G 
{X > 0, Y > F(X)}ni? e if g is small. Denote by Aq the operator obtained by 
fixing the coefficients of A in (6.51) at (X,Y) = (0,0). Then A = a 22 (0,0)A 
by (6.55). By (6.22), we obtain a 22 (0,0) = a 22 (0,y Pl ) > l/(4c 2 ) > 0. Now, by 
an explicit calculation and using (6.48), (6.55)-(6.57), (6.60), and (6.63), we 
find that, for 5 G (0, 5 ) and {X, Y) G {X > 0, Y > F(X)} n B e , 

Aw(r, 9) = a 2 (0, 0)Aw(r, 9) + (A — A )w(r, 9) 

< a 22 (0, 0)^" 1 ((a + l) 2 - (^) 2 ) cos(G(0)) 

+Cr a - 1 (^|a 11 (X,y) - on (0,0)| + |a 22 (X,Y) - a 22 (0,0)| 
+-r a - 1 |oi2(X,y)| + -r a |ai(X,y)| + Cr a |a 2 (X,y)| 

< ,.<,-, / (l-a)(5 + 3a) l-« «£\ 

for sufficiently small g > depending only on the data and 5. 

Thus, all the estimates above hold for small So > and g > depending 
only on the data. 

Now, since 

min w(X,Y) = L 4 > 0, 

{A>0, Y>F(X)}ndB e 
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we use the Comparison Principle (Lemma 6.3) (which holds since condition 
(6.52) satisfies (6.59) and b 3 < by (6.58)) to obtain 

l-(si+m)_ w ^ y i n {x > o, y > F(X)} n B e . 



u 

Similar estimate can be obtained for —v. Thus, using (6.9), we obtain (6.46) 
in B g . Since g depends only on the data and S > 0, then we use (6.9) to obtain 
the full estimate (6.46). 

Step 4- Estimate (6.45) can be obtained from (6.8), (6.20), and (6.46), 
combined with rescaling from the balls B d j L (z) n Q. + {<j)) for z G Clf(4>) \ {Pi} 
(with d z = dist(z, Pi) and L sufficiently large depending only on the data) into 
the unit ball and the standard interior estimates for the linear elliptic equations 
and the local estimates for the linear Dirichlet and oblique derivative problems 
in smooth domains. Specifically, from the definition of sets JC and £l + (4>) and 
by (5.16), there exists L > 1 depending only on the data such that 

B d /h{z) n (dn + ((f>) \ T shock ) = for any z G T shock D Q e , 

and 

B d/L {z) n (dQ + ((i)) \ T sonic ) = for any z G T sonic n Sl e . 

Then, for any z G ^ + (<^>) H B e (Pi), we have at least one of the following three 
cases: 

(1) Bx^Cfi+Wi 

10L 

(2) z G Bd^(zi) and £>- G (±,2) for some z\ G T sonic ; 

(3) z G %(zi) and jk- G (|,2) for some zi G r s/locfc . 

2i Z1 

Thus, it suffices to make the C 2 ' a -estimates of tp in the following subdo- 
mains for zq = (xo,yo) : 

(i) Bd I£L (z ) when ^^(zo) C Q + ((/>); 

201, 101/ 

(ii) b^(z ) n for z g r sonic n B e (Pi); 

2i 

(iii) B^zq) n for z g r s/locfc n 

2L 

We discuss only case (iii), since the other cases are simpler and can be 
handled similarly. 

Let zo G T s hock H Bg(Pi). Denote d = ^ > 0. Without loss of generality, 
we can assume that d < 1. 
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We rescale z = (x,y) near zq: 

Z = (X, Y) := \(x - x Q , y - y ). 
a 

Since B^(zq) n (dQ + ((j)) \ T s h oc k) = ; then, for p G (0, 1), the domain obtained 
by rescaling Q + (4>) Pi B^(zq) is 

ft* := B p n {Y < F(X) := U** + ~ U*o) ^ 

where is the function in (5.25). Note that yo = /^(xo) since (xo, yo) S r s ^ oc fc. 
Since L > 1, we have 

WHc^i-m < \\fd ( CS m 

and H/0II2 jj^'^' is estimated in terms of the data by (5.26). 
Define 

(6.64) v(Z) = -J—^{z + dZ) for Z G Of. 

d 1+a 

Then 

(6.65) IMLoo^o) < C 

by (6.46) with C depending only on the data. 

Since ip satisfies equation (6.19) in Q, + (4>) f^V'^ and the oblique derivative 
condition (6.6) on T s h oc k H T> ' 4e , then v satisfies an equation and an oblique 
derivative condition of the similar form in £1*° and on <9f2^° n {Y = F(X)}, 
respectively, whose coefficients satisfy properties (6.8) and (6.21) with the same 
constants as for the original equations, where we have used d < 1 and the 
C a / 2 -estimates of the coefficients of the equation depending only on the data, 
8, and ip. Then, from the standard local estimates for linear oblique derivative 
problems, we have 

|| 7 , 11 < n 

l|t7|l C=W*(ft'0 2) - ^ 

with C depending only on the data, 5, and ip. 

We obtain similar estimates for cases (i)-(ii), by using the interior esti- 
mates for elliptic equations for case (i) and the local estimates for the Dirichlet 
problem for linear elliptic equations for case (ii). 

Writing the above estimates in terms of ip and using the fact that the 
whole domain n B s (Pi) is covered by the subdomains in (i)-(iii), we 

obtain (6.45) by an argument similar to the proof of [20, Theorem 4.8] (see 
also the proof of Lemma A. 3 below). □ 
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Lemma 6.8. There exist C > and 5q G (0, 1) depending only on the data 
such that, if a,e > and M±,M2 > 1 in (5.15) satisfy (5.16), and 5 G (0,#o), 
i/iere exists a unique solution tp G C2 a/2 °/(6-13) and (5.30)-(5.33). T/ie 
solution ip satisfies (6.9) -(6. 10). 

Proof. In this proof, for simplicity, we write f2 + for + (</>) and denote by 
Ti, T2, T3, and the relative interiors of the curves r s ^ oc fe(0), £o(</>), r«>e<2ge> 
and r son j c respectively. 

We first prove the existence of a solution for a general problem V of the 

form 

2 2 
^ OijDj^ = / in n + ; b^Dii; = g t on T k , k = 1, 2, 3; V = on T D , 

i,j=l i=l 

where the equation is uniformly elliptic in f2 + and the boundary conditions on 
Tfc, fc = 1,2,3, are uniformly oblique, i.e., there exist constants Ai,A2,A3 > 
such that 

2 

A1M 2 < ^ diji^rfmnj < A^Vl 2 for a11 i€,v) g /i g R 2 , 
2 



i=i 

> A 3 for fc = 2,3, 



(k) 

and + Wi a(f^) — ^ f° r some L > 0. 

First we derive an apriori estimate of a solution of problem V . For that, 
we define the following norm for tp G C fc '^(0 + ), A; = 0, 1, 2, ... , and (5 G (0, 1): 

3 

*,k,p := E ll^llfc,/3^B 2e (p,)nf2+ + ll^llfc,/3!e 2e (p l )n^ + + ''^''c ,fc .' 3 (n+\(uf =1 B £1 (P l )))' 

i=2 i=l,4 

where g > is chosen small so that the balls B2 Q (Pi) for i = 1,...,4 are 
disjoint. Denote C* M := {ip G C7*' fc,/3 : \\^\U, k ,/3 < 00}. Then C*'*^ with 
norm || • ||* fc^g is a Banach space. Similarly, define 

3 

II II - II \\-P>{ p *} _lV^ 11 ii 1 -/ 3 -!^} ,11 11 

\\9k\\*,p ■- W gk Wk^,B 2e {P,)nv k + ll5,fc|| fc,/3,B 2e (P0nr fc + ll5' fc llc i .' 3 (r fc \(utiS e («)))' 

i=2 i=l,4 

where the respective terms are zero if B2 Q {Pi) fl r& = 0. Using the regularity 
of boundary of from the localized version of the estimates of [33, Theorem 
2] applied in B2 r (Pi) n + , z = 1,4, and of the estimates of [35, Lemma 
1.3] applied in B2 r (Pi) H i = 2,3, and the standard local estimates for 
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the Dirichlet and oblique derivative problems of elliptic equations in smooth 
domains applied similarly to Step 4 in the proof of Lemma 6.7, we obtain that 
there exists (3 = f3(£l + , A 2 , A 3 ) G (0,1) such that any solution V G C^(fi+) n 
0^(0+ \ T D ) n C 2 (Q + ) of problem V satisfies 

3 

(6.66) ||Y>||*,2,/3 < C(\\f\U t0 ,p + Y, \\9k\kf3 + M\o,n + ) 

fe=i 

for C = C(Q + , Ai, A2, A3, L). Next, we show that if) satisfies 

3 

(6-67) HV'II^^Cdl/IU^ + E Hs*M 

fc=i 

for C = C(f2 + , Ai, A2, A3, L). By (6.66), it suffices to estimate ||"0||o,n+ by the 
right-hand side of (6.67). Suppose such an estimate is false. Then there exists 
a sequence of problems V m for m = 1, 2, . . . with coefficients a™ and bf^ ,m , the 
right-hand sides f m and 5™, and solutions tp" 1 G C*' 2 '^, where the assumptions 
on and 6^' m stated above are satisfied with uniform constants Ai, A2, A3, 
and L, and ||/ m ||*, ,/3 + ELi lldk/3 as m -> 00, but ||Y> m || ,n+ = 1 for 

m = 1, 2, Then, from (6.66), we obtain ||V' m ||*,2, / 3 < C with C independent 

of m. Thus, passing to a subsequence (without change of notation), we find 
a% - a°- in C^/ 2 (tt+), &f )l,n - bf } '° in C 1 ^ 2 ^), and V> m - V> in C*> 2 ^ 2 , 
where ||V ;0 ||o,n+ = 1, and a°- and b\ k ^'° satisfy the same ellipticity, obliqueness, 
and regularity conditions as a™!- and b^' m . Moreover, ip° is a solution of the 
homogeneous Problem V with coefficients a^. and Since ||^ ||o,rj+ = 

1, this contradicts the uniqueness of a solution in C*' 2, @ of problem V (the 
uniqueness for problem V follows by the same argument as in Lemma 6.3). 
Thus (6.67) is proved. 

Now we show the existence of a solution for problem V if C in (5.16) is 
sufficiently large. We first consider problem Vq defined as follows: 

AV' = /in^ + ; D u ip = g k on T k , k = 1,2, 3; Y> = on T D . 

Using the fact that T2 and T3 lie on 77 = and 77 = £ tan 6 W respectively, and 
using (3.1) and (5.24), it is easy to construct a diffemorphism 

F: n+^Q:={(X,Y)e(0,l) 2 } 

satisfying 

H F llc^(f! + ) - C ' ll F_1 Hci.«(Q) ^ C > 
F(r D ) = S D := {1 = 1,7 6(0,1)}, 

and 

(6-68) \\DF~ 1 - Id\\ Ca{Qn{x<Vl/2}) < Ce 1 ^, 



GLOBAL SOLUTIONS OF SHOCK REFLECTION BY LARGE-ANGLE WEDGES 61 



where C depends only on the data, and (£1,771) are the coordinates of P\ 
defined by (4.6) with 771 > 0. The mapping F transforms problem Vq into the 
following problem Vq- 

2 

Di(dijDju) = f in Q; 

2 

y aijDjU Vi = g k on I k , k = 1, 2, 3; 
u = on Ed, 

where I k = F(Gk) are the respective sides of dQ, v is the unit normal on 
Ik, \\0"ij\\c^(Q) — C, and dij satisfy the uniform ellipticity in Q with elliptic 
constant A > 0. Using (6.68), we obtain 

(6-69) \\d i:j - Sf ||o(Qn{x<W2}) ^ ° £l/4 > 

where 8\ = 1 and 5\ = for i ^ j, and C depends only on the data. If e > 
is sufficiently small depending only on the data, then, by [13, Theorem 3.2, 
Proposition 3.3], there exists (3 G (0,1) such that, for any / G C l3 (Q) and 
c/k £ C {Ik) with k = 1,2,3, there exists a unique weak solution u G if 1 (Q) 
of problem Vq, and this solution satisfies u G C l3 (Q) D C 1,l3 (Q \ Sd)- We 
note that, in [13, Theorem 3.2, Proposition 3.3], condition (6.69) is stated in 
the whole Q, but in fact this condition was used only in a neighborhood of 
h = {0} x (0,1), i.e., the results can be applied to the present case. We 
can assume that (3 < a. Then, mapping back to Q + , we obtain the existence 
of a solution ib G C^(J2+) n C 1 'P(tt+ \ T D ) n C 2 (ft + ) of problem V for any 
/ G C 3 ^" 1 ") and 5^ G C^(Tk), k = 1,2,3. Now, reducing f3 if necessary and 
using (6.67), we conclude that, for any (f,g±, 92,93) G 3^ ■= {(/, 51,52,53) ■ 
||/||*,0,/3 + J2k=i \\9k\\*,/3 < °°}> there exists a unique solution ip G C*' 2 '^ of 
problem Vq, and ^ satisfies (6.67). 

Now the existence of a unique solution ib G C*' 2 ^ of problem "P, for any 
(7)51)52,53) G 3^ with sufficiently small (3 G (0, 1), follows by the method of 
continuity, applied to the family of problems tV + (1 — t)Vo for t G [0, 1]. This 
proves the existence of a solution ib G C*' 2 ' 13 of problem (6.13) and (5.30)- 
(5.33). 

Estimates (6.9)-(6.10) then follow from Lemma 6.4. The higher regularity 
ib G ci 1 /n follows from Lemmas 6.5-6.7 and the standard estimates for 
the Dirichlet problem near the flat boundary, applied in a neighborhood of 
r 'sonic \ (B g/2 (Pi) U B e/2 {Pi)) in the (x, y)-coordinates, where g > may 
be smaller than the constant g in Lemmas 6.6-6.7. In fact, from Lemma 
6.6, we obtain even a higher regularity than that in the statement of Lemma 
6.8: V £ C*2 a/2 < n+ I (4>) P3 ' P4 ^ • T ne uniqueness of solutions follows from the 
Comparison Principle (Lemma 6.3). □ 
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Lemma 6.8 justifies the definition of map J in (6.12) denned by J{ip) := ip. 
In order to apply the Leray-Schauder Theorem, we make the following apriori 
estimates for solutions of the nonlinear equation. 

Lemma 6.9. There exist C > and So G (0, 1) depending only on the data 
such that the following holds. Let a, e > and M±,M 2 > 1 in (5.15) satisfy 
(5.16). Let 5 G (0, <5o) and fi G [0,1]. Let tp G a / 2 n+(<^) ^ e a so ^ u ^ on °f 
(6.1), (5.30)-(5.32), and 

(6.70) Vr, = -M^2 on So(0) == d^ + (0) n {r? = -v 2 }. 

Then 

(i) There exists C > independent of ip and fi such that 

ll^llc i .«(n+(0)) - C ' ; 

(ii) ^ satisfies (6.9) -(6. 10) mi/i constant C depending only on the data; 

(hi) ^ G ^2ail+U) ■ Moreover, for every s G (0, c 2 /2), estimate (6.11) ao/cfe 
uraffo constant C depending only on the data and s; 

(iv) Solutions of problem (6.1), (5.30)^(5.32), and (6.70) satisfy the following 
comparison principle: Denote by J\f$(?p), Bi(ip), B2(ip), and Bz(ijj) the 
left-hand sides of (6.1), (5.30), (5.32), and (6.70) respectively. If 'ipi,ip 2 £ 

C 2~a,n+W satis fy 

B k {ipi) < B k (tp 2 ) on T shock (4>), T wedge , and £ (</>) /or k = 1, 2, 3, 

^1 > ^2 on r soraic , 

£/ien 

01 > ^2 ^ + ((/>). 

In particular, problem (6.1), (5.30)-(5.32), and (6.70) aas ai most one 
solution ^Cj^. 

Proof. The proof consists of six steps. 

Siep i. Since a solution V G C^^+'^j of (6.1), (5.30)-(5.32), and (6.70) 
with jj, G [0, 1] is the solution of the linear problem for equation (6.13) with 
■0 := if; and boundary conditions (5.30)-(5.32) and (6.70). Thus, estimates 
(6.9)-(6.10) with constant C depending only on the data follow directly from 
Lemma 6.4. 

Step 2. Now, from Lemma 5.2(h), equation (6.1) is linear in f2 + (0) H 
{c 2 - r > 4e}, i.e., (6.1) is (6.13) in ft + (0) n {c 2 - r > 4e}, with coeffi- 
cients aij(^,n) = AjA£,rf) + 5Sij for A\- defined by (5.35). Then, by Lemma 
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5.2(H), aij G C a (Q + ((p) n {c 2 — r > Ae}) with the norm estimated in terms of 
the data. Also, T s h oc k{<p) and the coefficients 6j of (6.3) satisfy (5.24) and 
(6.4)— (6.5). Then, repeating the proof of Lemma 6.5 with the use of the 
L°° estimates of tp obtained in Step 1 of the present proof, we conclude that 

for C depending only on the data. 

Step 3. Now we prove (6.11) for all s G (0,c 2 /2). If s > 6e, then (6.11) 
follows from (6.71). Thus, it suffices to consider the case s G (0, 6e) and show 
that 

( 6 - 72 ) ll^llc2, a( n+Wn{s/2<C2 _ r<6£+V4} ) < C(s)a, 

with C depending only on the data and s. Indeed, (6.71)-(6.72) imply (6.11). 

In order to prove (6.72), it suffices to prove the existence of C(s) depending 
only on the data and s such that 

( 6 - 73 ) ll^llo".-(S^)^ C 'WI^IU-(* / .W) 

for all z := (£,77) G + ((/)) n{s/2 < c 2 - r < 6e + s/4} with 6M(z,dSl + (4>)) > 
s/8 and that 



(6.74) IMIc*.»(B. /8 (*)nn+(*)) - ^WllVllL-CB.^Wnn+W) 

for all z G (r sftocfc (4>) U r werf9e ) n {s/2 < c 2 - r < 6e + s/4}. Note that all 
the domains in (6.73) and (6.74) lie within £l + (<p) n {s/4 < c 2 - r < 12e}. 
We can assume that e < c 2 /24. Since equation (6.1) is uniformly elliptic in 
fl + ((p) n {s/4 < c 2 — r < 12e} by Lemma 5.2(i), and the boundary conditions 
(5.30) and (5.32) are linear and oblique with C 1 '"-coefficients estimated in 
terms of the data, then (6.73) follows from Theorem A.l and (6.74) follows 
from Theorem A. 4 (in Appendix A). Since ||V ; l|L~(n+(vp)) < 1 by (6.9), the 
constants in the local estimates depend only on the ellipticity, the constants 
in Lemma 5.2(iii), and, for the case of (6.74), also on the C 2 ' a -norms of the 
boundary curves and the obliqueness and C 1,Q -bounds of the coefficients in the 
boundary conditions (which, for condition (5.30), follow from (5.24) and (6.4) 
since our domain is away from the points P\ and P 2 ). All these quantities 
depend only on the data and s. Thus, the constant C(s) in (6.73)-(6.74) 
depends only on the data and s. 

Step 4- I n this step, the universal constant C depends only on the data 
and 5, unless specified otherwise. We prove that tp G C 2 ' a (B e (P±) n$7+(</>)) 
for sufficiently small g > 0, depending only on the data and 5, and 

( 6 - 75 ) IMIc»( B.(P 4 )nn+(fl ) ^ C - 
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We follow the proof of Lemma 6.6. Since B Q (P/±) n C T>' for small 

g, we work in the (x, y)-coordinates. We use the notations £>+ and B+ + , 
introduced in Step 1 of Lemma 6.6, and consider the function 

v(x,y) = -tp(gx,gy). 
Q 

Then, by (6.10), v satisfies 



a 



(6.76) ^ 2C ~ £ ^ 

where the last inequality holds if C in (5.16) is sufficiently large. Moreover, v 
is a solution of 

(6.77) A[fv xx + 2A[ e 2 ) v xy + A i 2 e 2 ) v yy + A[ e) v x + A 2 \ 0) v y = O in B++ , 

(6.78) v = on 5 2 n{x = 0,y > 0}, 

(6.79) u„ = v y = on B 2 n {y = 0, s > 0}, 

with (A^\A\ e) ) = (A^\A\ e) )(Dv,x,y), where we use (6.2) to find that, for 

(x,y) e B+ + ,pe R 2 , i,j = 1,2, 

(6.80) 

(p, y) = An (p, £>x, £>y) + <5, 

^(p^.y) = ^(p^y) = ^i2(p,££,£y), 

^(P'^y) = ^22 (p, Qx,gy) + rg, 

(c 2 - px) 2 

A^ } (p,x,y) = pAi(p, £>x,£>y) + -, A 2 e) (p,x,?/) = gA 2 (p, gx, gy), 

(C 2 - £>Xj 

with Ay and Aj as in Lemma (5.3). Since £> < 1, Aj| } and Aj rf satisfy the 
assertions of Lemma 5.3(i)— (ii) with the unchanged constants. Moreover, A[f , 
A 22 , and A[ 8 ^ satisfy the property in Lemma 5.3(iii). The property in Lemma 
5.3(iv) is now improved to 
(6.81) 

\(A[tAitAi s) )(x,y)\ < Cq\x\, \D(A[tA^,Ai s) )(x,y)\ < Clgx^ 2 . 

Combining the estimates in Theorems A.l and A.3-A.4 with the argument 
that has led to (6.34), we have 

where C depends only on the data and 5 > by (6.76), since A\f and Aj e) 
satisfy (A.2)-(A.3) with the constants depending only on the data and 5. In 
particular, C in (6.82) is independent of g. 

We now use the domain D ++ introduced in Step 2 of the proof of Lemma 
6.6. We prove that, for any g G C a (D ++ ) with || #11 0(75++) < 1) there exists a 
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unique solution w G C 2 ' a {D ++ ) of the problem: 

(6.83) A[^w xx + A^w yy + A[ e) w x = g in D ++ , 

(6.84) w = on <9L> ++ n {x = 0, y > 0}, 

(6.85) w u = w y = on <9L> ++ n {x > 0, y = 0}, 

(6.86) w = v on <9D ++ n{x > 0,y > 0}, 

with (Ajf,^) = (A\f,A[ e) )(Dw,x,y). Moreover, we show 
(6-87) lkllc=».«(S++) ^ C > 

where C depends only on the data and is independent of g. For that, similar 
to Step 2 of the proof of Lemma 6.6, we consider the even reflection D + of 
the set D ++ , and the even reflection of (v, g, A[f , , A[^ ) from to 
B2, without change of notation, where the even reflection of (A[f , A$ , A[ e) ) , 
which depends on (p,x,y), is defined by 

Af (p,x,-y)= Af (p, x,y), A[ e) (p, x, -y) = i< e) (p, x, y) for (x, y) G B++ . 

Also, denote by v the restriction of (the extended) v to dD + . It follows 
from (6.78)-(6.79) and (6.82) that v G C7 2 ' Q (<9L> + ) with 

(6.88) ||t>||c".-(8£>+) < C, 

depending only on the data and 5. Furthermore, the extended g satisfies 

g G C a (D+) with \\g\\ Ca(BT) = \\g\\ c ^(D++) < 1- Th e extended 

and satisfy (A.2)-(A.3) in D + with the same constants as the estimates 
satisfied by An and Ai in S7 + (0). We consider the Dirichlet problem 

(6.89) A ( 1 fw xx + A^w yy + A ( f ) w x =g in £>+, 

(6.90) w = v on dL>+, 

with (A\f,A[ g) ) := (A^, A[ e) )(L>u;, x, y). By the Maximum Principle, 

IMIl~(d+) < \\v\\l~(d+)- 

Thus, using (6.88), we obtain an estimate of HHI-L 00 (£>+)• Now, using Theorems 
A.l and A. 3 and the estimates of ||s||cw;d+) an d ||^||c 2 ' q (9D+) discussed above, 
we obtain the a-priori estimate for the C 2 ' a -solution w of (6.89)-(6.90): 

(6-91) IMIc*.«.(5+) - °> 

where C depends only on the data and 5. Moreover, for every w G C 1,a (D + ), 
the existence of a unique solution w G C 2,a {D + ) of the linear Dirichlet problem, 
obtained by substituting w into the coefficients of (6.89), follows from [20, 
Theorem 6.8]. Now, by a standard application of the Leray-Schauder Theorem, 
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there exists a unique solution w G C 2,a (D + ) of the Dirichlet problem (6.89)- 
(6.90) which satisfies (6.91). 

From the structure of equation (6.89), especially the fact that A[f , A^ , 
and are independent of p2 by Lemma 5.3 (iii), and from the symmetry 
of the domain and the coefficients and right-hand sides obtained by the even 
extension, it follows that w, defined by w(x,y) = w(x, —y), is also a solution 
of (6.89)-(6.90). By uniqueness for problem (6.89)-(6.90), we find w(x,y) = 
w(x,—y) in D + . Thus, w restricted to D ++ is a solution of (6.83)-(6.86), 
where (6.84) follows from (6.78) and (6.90). Moreover, (6.91) implies (6.87). 

The uniqueness of a solution w G C 2,a (D ++ ) of (6.83)-(6.86) follows from 
the Comparison Principle (Lemma 6.3). 

Now we prove the existence of a solution w G C 2,a (D ++ ) of the problem: 
(6.92) 

ilf W XX + 2ig W x y + A { £ Wyy + W x + A^ Wy=0 in D + + , 

w = on dD++ n {x = 0, y > 0}, 

w v = w y = on dD ++ n {y = 0, x > 0}, 

w = v on dD ++ n {x > 0, y > 0}, 

where (A^,A^) := (A^ , A^)(Dw,x,y). Moreover, we prove that w satis- 
fies 

(6.93) IMIc».»(5++) - C 

for C > depending only on the data and 5. 
Let iV be chosen below. Define 

(6.94) S(N) := [w G C 2 ' a (D++) : || W|| C2 , a(STT) < iv} . 

We obtain such w as a fixed point of the map K : S(N) — > 5(iV) defined 
as follows (if ii is small and ./V is large, as specified below). For IF G S(N), 
define 

(6.95) 5 = -2ig(x,y)^ - A { 2 Q \x, y)W y . 
By (6.81), 

||S|| C «.(S++) < CNJq<1, 

if q < Qo with = ca72) for C depending only on the data and 5. Then, 
as we have proved above, there exists a unique solution w G C 2,a (D ++ ) of 
(6.83)-(6.86) with g defined by (6.95). Moreover, w satisfies (6.87). Then, if 
we choose N to be the constant C in (6.87), we get w G S(N). Thus, N is 
chosen depending only on the data and S. Now our choice of go = -^jp and 
q < go (and the other smallness conditions stated above) determines g in terms 
of the data and 5. We define K[W] := w and thus obtain K : S{N) -> S(N). 
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Now the existence of a fixed point of K follows from the Schauder Fixed 
Point Theorem in the following setting: From its definition, S(N) is a compact 
and convex subset in C 2 > a l 2 {D++). The map K : S(N) — > S(N) is continuous 
in C 2 ^I 2 (D++): Indeed, if W k G S{N) for k = 1, . . . , and W k -► VF in 
C 2 ' o/2 (L>++), then it is easy to see that VF G 5(7V). Define and 5 by (6.95) 
for VFfe and W respectively. Then g k — > 5 in C a / 2 (-D++) since (^12, ^2) = 
(-A12, ^X^) y) by Lemma 5.3(iv). Let Wk = K[W}.]. Then G S(N), and 
S(N) is bounded in C 2 ' a (Z) ++ ). Thus, for any subsequence u;^, there exists 
a further subsequence converging in C 2 ' Q//2 (D ++ ). Then the limit w 

is a solution of (6.83)-(6.86) with the limiting function g in the right-hand 
side of (6.83). By uniqueness of solutions in S(N) to (6.83)-(6.86), we have 
w = K[W]. Then it follows that the whole sequence A"[Wfc] converges to X[W]. 
Thus K : S(N) — > S(N) is continuous in C 2 ' a / 2 (D ++ ). Therefore, there exists 
w £ S(N) which is a fixed point of K. This function w is a solution of (6.92). 

Since v satisfies (6.77)-(6.79), it follows from the uniqueness of solutions 
in C(D++) n C 1 (D++ \{x = 0}) n C 2 {D ++ ) of problem (6.92) that w = v in 
D++. Thus, v e C 2 ' a (L>++) and satisfies (6.75). 

Step 5. It remains to make the following estimate near the corner Pi: 

where C depends only on the data, a, and 5. 

Since ip is a solution of the linear equation (6.13) for ^ = ip and satisfies the 
boundary conditions (5.30)-(5.33), it follows from Lemma 6.7 that satisfies 
(6.46) with constant C depending only on the data and 5. 

Now we follow the argument of Lemma 6.7 (Step 4): We consider cases 
(i)-(iii) and define the function v(X,Y) by (6.64). Then ip is a solution of 
the nonlinear equation (6.2). We apply the estimates in Appendix A. From 
Lemma 5.3 and the properties of the Laplacian in the polar coordinates, the 
coefficients of (6.2) satisfy (A.2)-(A.3) with A depending only on the data 
and 5. It is easy to see that v defined by (6.64) satisfies an equation of the 
similar structure and properties (A.2)-(A.3) with the same A, where we use 
that < d < 1. Also, v satisfies the same boundary conditions as in the proof 
of Lemma 6.7 (Step 4). Furthermore, since ip satisfies (6.46), we obtain the 
L°° estimates of v in terms of the data and 5, e.g., v satisfies (6.65) in case 
(iii). Now we obtain the C 2 ' a -estimates of v by using Theorem A.l for case 
(i), Theorem A. 3 for case (ii), and Theorem A. 4 for case (iii). Writing these 
estimates in terms of ip, we obtain (6.96), similar to the proof of Lemma 6.7 
(Step 4). 

Step 6. Finally, we prove the comparison principle, assertion (iv). The 
function u = ipi — ip2 is a solution of a linear problem of form (6.13), (5.30), 
(5.32), and (5.33) with right-hand sides A/5 (Vr) —■N'si'4'2) and B k (if)\) — Bk(ip2) 
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for k = 1,2, 3, respectively, and u > on T son i c . Now the comparison principle 
follows from Lemma 6.3. □ 

Using Lemma 6.8 and the definition of map J in (6.12), and using Lemma 
6.9 and the Leray-Schauder Theorem, we conclude the proof of Proposition 
6.1. □ 

Using Proposition 6.1 and sending S — > 0, we establish the existence of a 
solution of problem (5.29)-(5.33). 

Proposition 6.2. Let a,e,M±, and M 2 be as in Proposition 6.1. Then 
there exists a solution tp G C(fi+(0)) n C 1 (fi+(0) \ T^~) n C 2 (Q + (<p)) of 
problem (5.29)-(5.33) so that the solution if) satisfies (6.9)-(6.11). 

Proof. Let S G (0,<5 ). Let i/> s be a solution of (6.1) and (5.30)-(5.33) 
obtained in Proposition 6.1. Using (6.11), we can find a sequence 5j for j = 
1, . . . and V> G C 1 (Q + ((()) \ T son i c ) n C 2 (f2 + (</>)) such that, as j — ► oo, we have 

(i) Sj - 0; 



(ii) ^ -► V in C 1 ^ (</>)) for every s G (0,c 2 /2), where f2+(0) = n 
{c 2 - r > s}; 

(iii) ^ — ► V hi C 2 (K) for every compact K C f2 + (0). 

Then, since each ip$. satisfies (6.1), (5.30), and (5.32)— (5.33), it follows that 
ip satisfies (5.29)— (5.30) and (5.32)— (5.33). Also, since each ipg. satisfies (6.9)- 
(6.11), V a l so satisfies these estimates. From (6.10), we conclude that ip G 
C(n+(4>)) and satisfies (5.31). □ 



7. Existence of the Iteration Map and Its Fixed Point 

In this section we perform Steps 4-8 of the procedure described in §5.6. In 
the proofs of this section, the universal constant C depends only on the data. 

We assume that <f> G /C and the coefficients in problem (5.29)-(5.33) are 
determined by <f>. Then the existence of a solution ijj G (7(0+ (0)) nC 1 (fi+((/>) \ 
T~ c ) n C 2 (n + (</>)) of (5.29)-(5.33) follows from Proposition 6.2. 

We first show that a comparison principle holds for (5.29)-(5.33). We use 
the operators J\f and M. introduced in (5.29) and (5.30). Also, for fi > 0, we 
denote 



Q + ^(<f>) := n {c 2 - r < »}, r» shock {4>) -.= v shock {4>) n {c 2 - r < M }, 

Kedge : = T wedge H {c 2 - r < /*}■ 
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Lemma 7.1. Let a, e, M\, and M2 be as in Proposition 6.2, and fj, G (0, n), 
where k is defined in §5.1. Then the following comparison principle holds: If 
^2 € C(fi+>"(0)) n C 1 \ Q n C 2 (n+^(<f>)) satisfy that 

tftyi) < W2) in n + >»(cj>), 
MM<MM onT» hock (<t>), 
du^i < d v fa on T^ edge , 

ipi > ^2 on T sonic and n {02 - r = fi}, 

then 

Proof. Denote £ M := 0+(</>) n {c 2 - r = /i}. If /x € (0, k), then <9ft + ^0) = 

-^sfeocfc(^) U ^ wedge U r so „j c U S^. 

From N(ipi) < J\T(ip2), the difference ^1 — ^2 is a supersolution of a linear 
equation of form (6.13) in f2 + ' M ((^) and, by Lemma 5.2 (i), this equation is 
uniformly elliptic in Q + '^((f)) n {02 — r > s} for any s G (0,/x). Then the 
argument of Steps (i)-(ii) in the proof of Lemma 6.3 implies that ipi — tp2 
cannot achieve a negative minimum in the interior of n {02 — r > s} 

and in the relative interiors of r^ ocfc (0)n{c2 — r > s} and r^, ed9e n{c2 — r > s}. 
Sending s — ► 0+, we conclude the proof. □ 



Lemma 7.2. A solution V> G C(O+(0)) n C 1 (fi+(«^) \ T sonic ) n C 2 (ft+(^)) 
0/ (5.29) -(5.33) is unique. 

Proof. If V>i and ^2 are two solutions, then we repeat the proof of Lemma 
7. 1 to show that ipi — ip2 cannot achieve a negative minimum in Q + ((f>) and in 
the relative interiors of Tshock{ ( t ) ) an d ^ wedge- Now equation (5.29) is linear, 
uniformly elliptic near Eo (by Lemma 5.2), and the function ipi — ip2 is C 1 up 
to the boundary in a neighborhood of So- Then the boundary condition (5.33) 
combined with Hopf 's Lemma yields that tpi — ip2 cannot achieve a minimum 
in the relative interior of So- By the argument of Step (iii) in the proof of 
Lemma 6.3, ip\ — ip2 cannot achieve a negative minimum at the points P2 and 
P3. Thus, ipi > ip2 in Q + (4>) and, by symmetry, the opposite is also true. □ 

Lemma 7.3. There exists C > depending only on the data such that, 
if a, e, Mi, and M 2 satisfy (5.16), the solution tp G C(Sl+(4>)) n C 1 (^+(0) \ 
T~)nC 2 (n + (4>)) 0/ (5.29)-(5.33) satisfies 

(7.1) < V(x, y) < 3 + x 2 m fi'(^) := ft + > 2e (0). 



70 GUI-QIANG CHEN AND MIKHAIL FELDMAN 

Proof. We first notice that tp > in Q + ((f>) by Proposition 6.2. Now we 
make estimate (7.1). Set 

w(x,y) ^^y/- 

We first show that w is a supersolution of equation (5.29). Since (5.29) 
rewritten in the (x, y)-coordinates in tl'((f>) has form (5.42), we write it as 

where 

MWO = (2x - (7 + l)xCi(— ))V>^ + — - i>x, 

X C2 

M 2 (^) = Ofi> XX + Oti> X y + Ofyyy ~ 0\^ X + O^y. 

Now we substitute w(x,y). By (5.37), 

C l^L\ = t ( 6 ^ = 6 
Ul x > Ul 5( 7 + 1) ; 5(7 + 1)' 

thus 

Using (5.44), we have 

l-^HI = \-^ Y) Oi(Dw J x J y) + -J^O^^x,^! < Cx^ < Ce^x, 

where the last inequality holds since x G (0, 2e) in Q'((f)). Thus, if e is small, 
we find 

M{w) < in n'(4>). 

The required smallness of e is achieved if (5.16) is satisfied with large C. 

Also, w is a supersolution of (5.30): Indeed, since (5.30) rewritten in the 
(x, y)-coordinates has form (6.6), estimates (6.8) hold, and x > 0, we find 

g „ 3 

M{w) = b 1 (x,y)- — — x + b 3 (x,y)- — — x 2 < on V shock ((j)) D V . 
5( 7 + l) 5(7 + 1) 

Moreover, on T we d ge , w u = w y = = ip v . Furthermore, w = = ip on 
r sonic and, by (6.9), V < w on {x = 2e} if 

Ca < e 2 , 

where C is a large constant depending only on the data, i.e., if (5.16) is satisfied 
with large C. Thus, ip < w in by Lemma 7.1. □ 
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We now estimate the norm IIV'II^*"*^/^ m the subdomain &'((/)) := f2 + (</>)H 
{c 2 - r < e} of tf(<f>) := n {c 2 - r < 2e}. 

Lemma 7.4. There exist C, C > depending only on the data such that, 
ifa,£,M 1 , and M 2 satisfy (5.16), the solution ip G C(Q + ((p)) n C 1 (fi+(0) \ 
TW)nC7 2 (0+ (</»)) o/(5.29) -(5.33) safe/?es 

Proof. We assume C in (5.16) is sufficiently large so that <r, e, Mi, and M 2 
satisfy the conditions of Lemma 7.3. 

Step 1. We work in the (x, y)-coordinates and, in particular, we use 
(5.25)-(5.26). We can assume e < k/20, which can be achieved by increasing 
C in (5.16). 

For z := (x,y) G Cl'(4>) an d P e (0> 1)) define 
(7-3)_ 

Rz, P ■= {(«,*) : |« - ar| < |ar, I* - y| < ^v 7 ^} , R z , p := R ZtP nQ + ((/)). 



Since fi'(^) = 0, + (<t>) n {c 2 - r < 2e}, then, for any 2 G fV(0) and p G (0, 1), 
(7.4) ^Cfl + Wn{( S) t) : |x<k|i}C% 

For any z G (</>), we have at least one of the following three cases: 

(i) #2,1/10 = #2,1/10! 

(ii) z G R Zw! \/2 for z w = (x, 0) G r wed9e ; 

(hi) 2 G i? 2s ,i/ 2 for z s = (xj^x)) G r s/locfc (0). 

Thus, it suffices to make the local estimates of Dijj and D 2 ip in the following 
rectangles with z := (x ,yo) : 

(i) #2 ,i/20 for z G n'(^) and i^ ,i/i = R Zo ,i/w; 

(ii) i? Zo ,i/2 for ^0 G r TOe d<?e n {x < e}; 
(hi) # 2o ,i/2 fo r e F shocA .(0) n {x < e}. 

5iep 2. We first consider case (i) in Step 1. Then 

#20,1/10 = {(^o + y 5, 2/0 + ^T) : {S,T)€Q 1/10 }, 
where Q p := (—p,p) 2 for p > 0. 
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Rescale ip in R ZO: i/io by denning 
(7.5) ^(S,T) := \^(x + ^S,y + ^T) for (S,T) G Q 1/w . 

X 4 4 

Then, by (7.1) and (7.4), 

(7-6) ll^ 0) llc(QT7^)< V(7 + l). 

Moreover, since ip satisfies equation (5.42)-(5.43) in -R 20i i/io, then ^ z °^ satisfies 

( 2 o) 

(7-7) ((1 + \S)(2 - ( 7 + l)Ci( j^)) + xoO^y^ + x O^Vst ] 

+ (- + xoO?' 2o) )Vtt " (t + xoO^ Zo) )4 Zo) + x 2 0^ Zo) 4 Zo) = 

in Qi/iO) where 
(7.8) 

dt z °(p,S,T) = -(i±f^ + g( 2 (l + 5/4)^1(1^4) - 16|4 20) | 2 ) 

-?(^ > + (,4 I ;w i^ ) i 2 )- 

Of >, 5, T) = -_ I _-^ W (4* ^ o) + c 2 - x (l + ^/4))^ o) , 
Ot°(p,S,T) 

\l + S/4)(2c 2 -x (l + S/4)) 



c 2 (c 2 -x (l+S/4)) 2 

-(7 - l)(x <f>M + (c 2 - x (l + 5/4))(l + 5/4)Ci( T ^ i ) + 8x |4 Zo) | 2 ) 
8(7+1) -21 ,/,(*>) 12 1 

(c2-x (1+S/4)) 2X 0\ ( Pt I J' 

+ (c 2 - x (l + 5/4)) (1 + S/4)Ci (^j) + { J X x °Ss l J W ) }' 
^(P. 5 T) = C2(C2 _ Xo f 1+6 . /4)) . (4xo^ o) + 2c 2 - 2xo(l + S/4))4 2 °\ 
where <^ 2 °) is the rescaled (f> as in (7.5). By (7.4) and <f> G /C, we have 

ll^ 0) llc-(07^)< CM i' 

and thus 

(7-9) II^1I C1( 7^ XR2) <^(1 + Mf), * = !,.. .,5. 

Now, since every term O^f'^ in (7.7) is multiplied by Xq* with > 1 and 
xo G (0, e), condition (5.16) (possibly after increasing C) depending only on 
the data implies that equation (7.7) satisfies conditions (A.2)-(A.3) in Q1/10 
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with A > depending only on C2, i.e., on the data by (4.31). Then, using 
Theorem A.l and (7.6), we find 

(7-10) H^ o) llc,»(^)<^ 

Step 3. We then consider case (ii) in Step 1. Let zq G T we( i ge H {x < e}. 
Using (5.25) and assuming that a and e are sufficiently small depending only 
on the data, we have R Zo ,i n dtl + (4>) C T we d ge and thus, for any p G (0, 1], 

^ ,P = {(*o + yS,yo + ^T) : (5,T)€Q p n{T>0}|. 

The choice of parameters for that can be made as follows: First choose a small 
so that |£ — £i| < |£|/10, where £ is defined by (3.3), which is possible since 
£i £ as # w — ► 7r/2, and then choose £ < (|£|/10) 2 . 

Define V (2o) (S,T) by (7.5) for (S,T) G Qi n {T > 0}. Then, by (7.1) and 
(7-4), 

(7-H) ll^ o) llc(Qrn { T>o})<V(7 + l)- 

Moreover, similar to Step 2, ?/>( Zo ) satisfies equation (7.7) in Q\ n {T > 0}, and 
the terms Op Za satisfy estimate (7.9) in Qi n {T > 0}. Then, as in Step 2, we 
conclude that (7.7) satisfies conditions (A.2)-(A.3) in Q 1 n {T > 0} if (5.16) 
holds with sufficiently large C. Moreover, since ip satisfies (5.32), it follows 
that 

d T i/,(*°) = Q on {T = 0}nQ 1/2 . 
Then, from Theorem A. 4, 

(7-12) llV< {zo) llc^(^n { T>o})<^ 

Step 4- We now consider case (iii) in Step 1. Let zq 6 ^shocki^) (~){x < e}. 
Using (5.25) and the fact that yo = f<t>{ x o) f° r £ r s / loc fc(0) n < e}, 
and assuming that a and e are small as in Step 3, we have R Zo ,i H Sf2 + (</>) C 
L s / loc fc(0) and thus, for any p G (0, 1], 

R Zo , P = \{x + ^S,y + ^T) : (5, T) G Q p n {T < e^F^S)}} 

with 

F ^ s) = 4 • 

Then we use (5.27) and xq G (0, 2s) to obtain 
^ o )(0) = 0, 

II F II < II4II^°°([0^])^0 , ,1/4 

ll^ )llcM[-i/2,i/2]) < £ i/4^ ^l + MieJe > 

II erf/ || ^ H/^llL°°([0,2e])gQ + [/^]a,(g /2, e )g§ +Q! . , 5/4 

ll^ o) llc-([-i/2,i/2]) < 4£ i/4 ^ ^ C(l + Mi)e , 
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and thus, from (5.16), 

( 7 - 13 ) ll^(z )llc 2 .<*([-l/2,l/2]) < C/C < 1 

if C is large. Define ^ Zo \S,T) by (7.5) for (S, T) G Q x n {T < e 1 / 4 F (zo) (5)}. 
Then, by (7.1) and (7.4), 

(7-14) ll^° ) llc(Q7n { T<^ 0) (5)})< 1 /(7 + l). 

Similar to Steps 2-3, V ( * o) satisfies equation (7.7) in Q 1 n {T < e 1 / 4 F (zo) (S')} 
and the terms Of' 20 satisfy estimate (7.9) in Q x n {T < e 1 / 4 F (zo) (,S')}. Then, 
as in Steps 2-3, we conclude that (7.7) satisfies conditions (A.2)-(A.3) in Qi n 
{T < e 1 / 4 ^)^)} if (5.16) holds with sufficiently large C. Moreover, V 
satisfies (5.30) on T ^^{(f) , which can be written in form (6.6) on Tghocki'P) n 
V. This implies that il)^ satisfies 

dsi) ( Zo ) = £ i/4 (fy&rtW + S 3 V (2o) ) on {T = e 1 / 4 F (zo) ( < S)} n Q 1/2 , 
where 

a fart V%o~ b -2 f , x o Q , v^o^x 
B 2 (S,T) = -- r j I ^-(x + — S,y + —T), 

^T) = - — ^.(xo + T< S,y + —T). 

From (6.8), 

11(^2, ^3)11^^^/^,(5)} < ^ /4 ^1 < < !" 

Now, if e is sufficiently small, it follows from Theorem A. 2 that 
( 7 - 15 ) ll^ ( ^ o) llc*2.-(Q77In{T< e i/4 F(zo) (5) } ) < C. 

The required smallness of e is achieved by choosing large C in (5.16). 

Step 5. Combining (7.10), (7.12), and (7.15) with an argument similar to 
the proof of [20, Theorem 4.8] (see also the proof of Lemma A. 3 below), we 
obtain (7.2). □ 

Now we define the extension of solution ip from the domain to the 

domain D. 

Lemma 7.5. There exist C, C\ > depending only on the data such that, 
if a, e, M±, and M 2 satisfy (5.16), there exists Cz{e) depending only on the data 
and e and, for any 4> £ IC, there exists an extension operator 

: c 1 ' a (n+(^)) n c 2 ' a (n+i^) \ r sonic u s ) - c^(p) n c 2 ' a (v) 

satisfying the following two properties: 
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(i) // V € C7 1 ' Q (^+(0)) n C 2 ' a (Q+(cb) \ T sonic U S ) is a solution of problem 
(5.29)-(5.33), then 

(7.16) ^ ^ 

(7.17) \nn[-^ a) < C2(e)o-, 

(ii) Let (3 G (0, a). If a sequence <pk G /C converges to (f) in C l, ^(V), then 
4> e K. Furthe rmore, if ^ G C 1 -°(n+ ( ^)) n C7 2 '°( Q+(c/> fc ) \ F]~ U S ) 
and V G C 1 - Q (0+((/))) n C 2 ' a {n+(<j>) \ r sonic US ) are i/ie solutions of 
problems (5.29) -(5.33) for 4>k and <fi respectively, then V^ipk ~~ * P^ i n 
C l ^{V). 

Proof. Let k > be the constant in (5.25) and e < n/20. For any 
G /C, we first define the extension operator separately on the domains fii := 
ri + (0) n {c2 — r < k} and O2 := ^ + (0) H {C2 — r > k/2} and then combine 
them to obtain the operator 7^ globally. 

In the argument below, we will state various smallness requirements on a 
and e, which will depend only on the data, and can be achieved by choosing 
C sufficiently large in (5.16). Also, the constant C in this proof depends only 
on the data. 

Step 1. First we discuss some properties on the domains Sl + (0) and V to 
be used below. Recall £ < defined by (3.3), and the coordinates (£1,771) of 
the point Pi defined by (4.6). We assume a small so that |£ — £i| < |£|/10, 
which is possible since £1 — > £ as W — > 7r/2. Then £1 < 0. By (5.24) and 
Pi G T shock (4>), it follows that 

(7.18) r s ^)cPn{(<a + e 1/4 }. 

Also, choosing e 1 / 4 < |£|/10, we have 

(7.19) £1 + e 1/4 < 1/2 < 0. 

Furthermore, when a is sufficiently small, 
(7.20) 

if (£, 77) g P n {£ < £1 + £ V4 }; ( ^ r?) G P; and ^/ > e> ^ en |^| < |£| . 

Indeed, from the conditions in (7.20), we have 

-c 2 <£<a+£ 1/4 <e/2<o. 

Thus, < |£| if £' < 0. It remains to consider the case £' > 0. Since 
P C -B C2 (0) n {£ < 77Cot6» TO }, it follows that < c 2 cos6> w . Thus < |£| if 
C2COS^ w < |£|/2. Using (4.31) and (3.1), we see that the last inequality holds 
if a > is small depending only on the data. Then (7.20) is proved. 
Now we define the extensions. 
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Step 2. First, on fii, we work in the (x, y)-coordinates. Then Q\ = {0 < 
x < k, < y < f^(x)} by (5.25). Denote Q( a> b) '■= (0, k) x (a, 6). Define the 
mapping $ : Q(_oo,oo) * Q(— oo,oo) 

by 

$(x,y) = (x,l-y/ft(x)). 

The mapping $ is invertible with the inverse 3> _1 (x, y) = (x, /<^(x)(l — y)). By 
definition of <£, 

$(fii) = Q (0 ,i), $(r s/locfc (^) n {o < x < «}) = (o, k) x {o}, 

(7.21) $(Dn{0<x< K })cQ K i), 

where the last property can be seen as follows: First we note that /</>(x) > 
> for x G (0, k) by (5.8) and (5.26), then we use V n {0 < x < re} = 
{0<x<re, 0<y< /o(x)} and (5.27) to obtain -J— > on V n {0 < x < re} 
and 

sup T7T = sup T7T - 1 + 7— ^H/*-/ollc(o,«) 

(x,j/)6Z>n{0<x</c} /^(x) x6(0,«) /^(x) /o,o(0) 

< 1 + C{M x e + M 2 ct) < 2, 

if Mi£ and M 2 cr are small, which can be achieved by choosing C in (5.16) 
sufficiently large. 

We first define the extension operator: 

£2 : C7^(^)nC 2 ^(Q^\{x = 0}) - C^Q^^C^Q^Ux = 0}) 

for any j3 G (0, 1]. Let v G C 1 '^(Q^) n C 2 '^(Q^"\ {x = 0}). Define £ 2 v = u 
in Q(o,i)- For 6 ^(-1,0)) define 



(7.22) e 2 v(x,y)=J2aiv(x,-K), 



1=1 



where ai = 6, a 2 = —32, and 03 = 27, which are determined by J2i=i a i (~ j)™ = 

1 form = 0,1,2. 

Now let i) G C 1 ' a (n+(^)) n C7 2 - Q (lF(^) \r sonic us ). Let 

Then v G C7 1 -' 3 (Q^)nC7 2 ' /3 (Q^\{x = 0}). By (7.21), we have Vn{c 2 -r < 
re} C 3> _1 (Q(-i,i))- Thus, we define an extension operator on Q\ by 

= (S 2 v) o $ on V n {c 2 - r < re}. 



Then T^V G C^Pi) D C 2 > a {V x \ T somc ) with V 1 := V (1 {c 2 - r < re}. 

Next we estimate 7-^ separately on the domains 2?' = 2? H {c 2 — r < 2e} 
and Pi n {c 2 - r > e/2}. 
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In order to estimate the Holder norms of on V', we note that $(fi'(<^)) = 
(0,2e) x (0,1) and V C $ -1 ((0,2e) x (-1,1)) in the (x, ^-coordinates. We 
first show the following estimates, in which the sets are defined in the (x, y)- 
coordinates: 

(7.23) U o ^Q^o,!) < CWntl'W fOT ^ * G C tk^V 

(7.24) \\ W o nth ± ^ll-HS,2 £ )x(-i,i) ^ w G < a ? , £)x( _ 1;1) , 

(7.25) || ^2 V 1 1 ^[(0,2e) x (—1, 1) - C 'll t 'll2X(0,2e)x(0,l) f ° r Wy V G ^a^O^e) x (-1,1) " 

To show (7.23), we denote v = ip ° < 1 )_1 and estimate every term in defi- 
nition (5.11) for v. Note that v(x,y) = ij>(x, fcp( x )(l ~ y))- I n the calcu- 
lations below, we denote (v,Dv,D v) = (v, Dv,D 2 v)(x,y), (ip,Dip,D 2 ip) = 
(^,D^,D 2 ^)(x,U(x)(l - y)), and #) = (UJ'^f^)(x). We use that, 

for x G (0,2e), < M\x < 2M l e < 2/C by (5.16). Then, for any (x,y) G 
(0,2e) x (0,1), we have 

\v\ = m < H\\ ( r r L^x\ 



1 2,a,n r w 

kl = IV, + (i - y^yfll < W\\t%(<t>) ( x + * 3/2 (! + M ^)) ^ ciM&w*. 

Kx| = Ife + 2(1 - y)^/; + (1 - yf^ yy {0 + (1 - 2/)^/;'| 



< IU/,l|(P ar ) 



1 + x x / 2 (l + Mix) + x(l + M x xf + Mix 3 / 2 ) 



2,a,Q'(4>) 

<cu\\tl w - 

The estimates of the other terms in (5.11) for v follow from similar straightfor- 
ward (but lengthy) calculations. Thus, (7.23) is proved. The proof of (7.24) is 
similar by using that f$(x) > /o,o(0)/2 > for x G (0, k) from (5.8) and (5.26) 
and that /o,o(0) depends only on the data. Finally, estimate (7.25) follows 
readily from (7.22). 

Now, let i) G C 1 ' a (f2+ ((f))) n C 2 ' a (n+((P) \ T sonic U S ) be a solution of 
(5.29)-(5.33). Then 

WK^PtZl = II Wl* o s- 1 ) o d>||t a :i < c|| o Oll 2 °5U)x ( -i,i) 

^^ll^l^o^H^x^l) < ^ll^llSw < ^ 
where the first inequality is obtained from (7.24), the second inequality from 

(7.25) , the third inequality from (7.23), and the last inequality from (7.2). 
Thus, (7.16) holds for V\. 

Furthermore, using the second estimate in (5.27), noting that M20 < 1 
by (5.16), and using the definition of and the fact that the change of 
coordinates (x, y) — > (£,77) 1S smooth and invertible in V n {e/2 < x < k}, we 
find that, in the (£, r/)-coordinates, 

(7.26) II^V'll^.o^n^^cj-r^K}) - C 'll^llc 2 .<»(n+(^)n{e/2< C2 -r<«})- 
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Step 3. Now we define an extension operator in the (£, r/)-coordinates. 

Let 

£ 2 : C\[0,1] x [-v 2 , m ])nC 2 ([0,l] x (-v 2 , m }) 

^ C \[-l,l] x [-v 2 , m })nc 2 ([-i,i] x {-v 2 , m }) 

be defined by 

3 X 

£ 2 v(X,Y) :=^Oiv(--,Y) for (X,Y) G (-1,0) x (-v 2 , m ), 
i=i * 

where ai, 02, and 03 are the same as in (7.22). 

Let Cl 2 := tt + (<p) n {0 < 77 < 771}. Define the mapping * : (l 2 -> (0, 1) x 
(-u 2 ,r?i) by 

where /,£(•) is the function from (5.21)-(5.22). Then the inverse of ^ is 

*-\X, Y) = (U(Y) + X(Y cot 9 W - U(Y)), Y), 
and thus, from (5.24), 

(7 27) ||^||(- 1 -«.[0,l]x{-^ 2 ,r)i}) ,|^ r -l||(-l-a,[0,l]x{-%,r)i}) < £, 

V ' ' 11 "2,0,^2 11 "2,a,(0,l)x(— u 2 ,»7i) — 

Moreover, by (5.24), for sufficiently small e and a (which are achieved by 
choosing large C in (5.16)), we have V n {—v 2 < rj < 771} C ^ _1 ([— 1,1] x 
[— v 2 ,rii\). Define 

V 2 ^ := 6 2 (ip o o * 011 PR {-v 2 < 7] < 771}. 

Then 7^ G C^P) n C 2 ' a (P \ r sonic U S ) since X> \ C P n {-u 2 < 

7/ < f]i}. Furthermore, using (7.27) and the definition of V^, we find that, for 
any s G (-v 2 ,f] 1 ], 

(7-28) r^llfc^n'St} ^ C ^ ~ 

where C(r]i — s) depends only on the data and 771 — s > 0. 

Choosing C large in (5.16), we have e < k/100. Then (5.25) implies that 
there exists a unique point P' = r s /j OC fc(</>) n {c 2 — r = k/8}. Let P' = (£', rf) 
in the (£, 7/)-coordinates. Then 7/ > 0. Using (7.18) and (7.20), we find 

(V \ Q + (</>)) n{c 2 -r> k/8} CVn{r]< r/}, 

n+(0) n {77 < t?'} c fi+(<£) n {c 2 - r > k/8}. 

Also, k/C <iji-)j'<Ck by (5.22), (5.24), and (4.3). These facts and (7.28) 
with s = rj imply 

y-^l II' </>^ll2,a,X>n{c 2 -r>K/8} - U H ^ll2,a,Q+(^)n{c 2 -r->/t/8}- 
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Step 4- Finally, we choose a cutoff function £ G C°°(R) satisfying 
C= 1 on (-oo,k/4), C = on (3/c/4,oo), ('<0 onR, 
and define 

:= C(c 2 - r)V\i> + (1 - C(c 2 - t))V%$ in V. 

Since V^ip = ip on Q + (cf)) for k = 1,2, so is "P^V- Also, from the properties of "P^' 
above, v^> g C l > a (V)nC 2 > a (V) if ^ g c l > a (n+^j)nc 2 > a (n+^j\r sonic U S ). 
If such -0 is a solution of (5.29)-(5.33), then we prove (7.16)-(7.17): = V\ip 
on V by the definition of £ and by e < k/100. Thus, since (7.16) has been 
proved in Step 2 for V^ip, we obtain (7.16) for V^tp- Also, ip satisfies (6.11) 
by Proposition 6.2. Using (6.11) with s = e/2, (7.26), and (7.29), we obtain 
(7.17). Assertion (i) is then proved. 

Step 5. Finally we prove assertion (ii). Let G /C converge to in 
C^tp). Then obviously <fi £ JC. By (5.20)-(5.22), it follows that 

(7.30) inC^a-^.TTi]), 
where /^ fc , G c\ 

a^-vl^t) ' V ^ are * ne f unc tions from (5.21) corresponding to 
<f> k ,<f>, respectively. Let Y> fe ,Y> G C 1 > a (fi+(^ fc )) n C 2 - ft (n+(4) \r sorac U S ) be 
the solutions of problems (5.29)-(5.33) for 0. Let {V'fc™} be any subsequence 
of {V'fc}- By (7.16)-(7.17), it follows that there exist a further subsequence 
{4>k m J and a function G C^ a (V) n C 2 '°(£>) such that 

^W m „ -0fe m „ ^ in C 2 ' a/2 on compact subsets of V and in C 1 ' / 2 ^). 

Then, using (7.30) and the convergence (f>k — ► in C 1, ^(P), we prove (by the 
argument as in [10, page 479]) that \p is a solution of problem (5.29)-(5.33) for 
4>. By uniqueness in Lemma 7.2, -ijj = ijj in f2 + (0). Now, using (7.30) and the 
explicit definitions of extensions and V 2 , it follows by the argument as in 
[10, pp. 477-478] that 

WO ^ C^(<W)) m C^(P), 
(1 - C)P| fcmn (^ m „ ) - (1 - ()-P 2 Mn + W ) in C^(P). 

Therefore, -ip = ip in V. Since a convergent subsequence {V'fc m „} can be ex- 
tracted from any subsequence {^fc m } of {ipk} and the limit tp = ip is indepen- 
dent of the choice of subsequences {ipk m } and {ipk m }> if follows that the whole 
sequence 0& converges to in C 1,/3 (2?). This completes the proof. □ 

Now we denote by Cq the constant in (5.16) sufficiently large to satisfy 
the conditions of Proposition 6.2 and Lemma 7.5. Fix C > Cq. Choose 
Mi = max(2Ci, 1) for the constant C\ in (7.16) and define e by (5.63). This 
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choice of e fixes the constant C 2 (e) in (7.17). Define M 2 = max(C2(e), 1). 
Finally, let 

_ C- 1 - e - e 1 / 4 M l 2 
a ° ~ 2 (Mf + e 2 max(Mi, M 2 )) & ' 

Then do > 0, since e is defined by (5.63). Moreover, cto, e, Mi, and M 2 depend 
only on the data and C. Furthermore, for any a G [0, o~o\, the constants a, e, 
Mi, and M 2 satisfy (5.16) with C fixed above. Also, ip > on by (6.9) 

and thus 

(7.31) V^ip > onP 

by the explicit definitions of V\,V\, and V</,. Now we define the iteration map 
J by J ((f)) = V^ip. By (7.16)-(7.17) and (7.31) and the choice of a, e, Mi, 
and M2, we find that J : /C — ► /C. Now, /C is a compact and convex subset of 
Ci,a/2(p). The map J : /C -► /C is continuous in C llQ / 2 (P) by Lemma 7.5(h). 
Thus, by the Schauder Fixed Point Theorem, there exists a fixed point (ft £ K, 
of the map J. By definition of J, such ip is a solution of (5.29)^(5.33) with 
cf) = tp. Therefore, we have 

Proposition 7.1. There exists Co > 1 depending only on the data such 
that, for any C > Co, there exist ao,£ > and Mi,M 2 > 1 satisfying (5.16) so 
that, for any a £ (0, o"o], there exists a solution ip £ K,{a, e, Mi, M 2 ) of problem 
(5.29)-(5.33) with cj) = tp (i.e., ip is a "fixed point" solution). Moreover, ip 
satisfies (6.11) for all s £ (0,c 2 /2) with C(s) depending only on the data and 
s. 

8. Removal of the Ellipticity Cutoff 

In this section we assume that Co > 1 is as in Proposition 7.1 which 
depends only on the data, C > Co, and assume that o~q, e > and Mi,M 2 > 1 
are defined by C as in Proposition 7.1 and a £ (0, ao]. We fix a "fixed point" 
solution ip of problem (5.29)— (5.33), that is, ip £ K{a,e,Mi,M 2 ) satisfying 
(5.29)-(5.33) with cp = ip. Its existence is established in Proposition 7.1. To 
simplify notations, in this section we write + , T s h oc k, and So for Q + (ip), 
Tshockitp), and Eo(^>), respectively, and the universal constant C depends only 
on the data. 

We now prove that the "fixed point" solution ip satisfies \ip x \ < Ax/ (3(7 + 
1)) in £l + n {c 2 — r < 4e} for sufficiently large C, depending only on the data, 
so that ip is a solution of the regular reflection problem; see Step 10 of §5.6. 

We also note the higher regularity of ip away from the corners and the sonic 
circle. Since equation (5.29) is uniformly elliptic in every compact subset of 
Sl + (by Lemma 5.2) and the coefficients Aij(p, £, 77) of (5.29) are C 1,a functions 
of (p, £, 77) in every compact subset of R 2 x $7 + (which follows from the explicit 
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expressions of Aij(p,£,r]) given by (5.35), (5.41), and (5.48)), then substituting 
p = Dip(£,rf) with ip e K, into Aij(j?,£,rf), rewriting (5.29) as a linear equation 
with coefficients being C 1 '" in compact subsets of and using the interior 
regularity results for linear, uniformly elliptic equations yield 

(8.1) V e C 3 ' a (n + ) . 

First we bound ip x from above. We work in the (x, y)-coordinates in 
n + n {c 2 - r < 4e}. By (5.25), 

(8.2) n + ((f>) n {c 2 - r < 4e} = {0 < x < k, < y < f<f,(x)}, 
where satisfies (5.26). 

Proposition 8.1. For sufficiently large C depending only on the data, 

(8.3) V* < 3 ( 7 + 1 ) a in n+ n & - 4e ^ 
Proof. To simplify notations, we denote A = 3 ( 7 4 fl ) and 

nf -.= n + n{x< s} for s > o. 

Define a function 

(8.4) v(x,y) :=Ax-ip x (x,y) on fi| £ . 
From i[} £ K, and (8.1), it follows that 

(8.5) t, g c ' 1 (?2+ ) n c 1 (n+ \ {x = o}) n c 2 (n+ ) . 

Since tp G /C, we have |^ x (x,y)| < Mix in fi^. Thus 

(8.6) v = on a^| £ n{x = 0}. 

We now use the fact that ip satisfies (5.30), which can be written as (6.6) in 
the (x, y)-coordinates, and (6.8) holds. Since ip G K, implies that 

\il>{x,y)\ < M lX 2 , \il> v {x,y)\ < M^ 3 / 2 , 

it follows from (6.6) and (6.8) that 

l^xl < C(|^| + IV-I) < CM lX ^ 2 on T shock n {x < 2e}, 

and hence, by (5.16), if C is large depending only on the data, then 

\ip x \ < Ax on T shock n {0 < x < 2e}. 

Thus we have 

(8.7) v>0 on T shock n {0 < x < 2e}. 
Furthermore, condition (5.32) on T we d ge in the (x, y)-coordinates is 

ip y = on {0 < x < 2e, y = 0}. 
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Since ip G K. implies that ip is C 2 up to T we d ge , then differentiating the condition 
on T we d ge with respect to x, i.e., in the tangential direction to T we d ge , yields 
ipxy = on {0 < x < 2e, y = 0}, which implies 

(8.8) v y = on r med9e n{0 < x <2e}. 
Furthermore, since ip £ 

(8.9) \tp x \ < M 2 a < Ae on Q + n {e/2 <x< 4e}, 

where the second inequality holds by (5.16) if C is large, depending only on 
the data. Thus, for such C, 

(8.10) v > on fi£ n {x = 2e}. 

Now we show that, for large C, v is a supersolution of a linear homogeneous 
elliptic equation on Q^e- Since ip satisfies equation (5.42) with (5.43) in Q^ £ , 
we differentiate the equation with respect to x and use the regularity of ip in 
(8.1) and the definition v in (8.4) to obtain 
(8.11) 

anv xx + ai 2 v xx + a 2 2V yy 

+ (A - v x ) ( - 1 + ( 7 + 1){CM -D+CM- |)(| - v x ))) = E(x, y), 

where 

(8.12) an = 2x-(7 + l)xCi(— )+Oi, ai2 = 2 , a 22 = - + 3 , 

X c 2 

(8.13) £(x, y) = ip xx d x 6\ + ip xy d x 6 2 + ip yy d x 6 3 - ip xx 6 A - ip x d x O A 

+1pxyO<j + 1pyd x 5 , 

with 

(8.14) O k (x,y) =0t(D 1 p(x,y),x,y) for fc = 1, . . . , 5, 
for Ojf defined by (5.43) with <j) = ip. From (5.37), we have 

Ci (A) = A. 
Thus we can rewrite (8.11) in the form 

(8.15) a u v xx + a 12 v xx + a 22 v yy + 6^ + cv = -A((j + 1)A - l) + E(x, y), 
with 

(8.16) b(x, y) = 1 - (7 + l)(Ci04 - -) + £(A - -)(- - v x - A)), 



(8.17) c (x,y) = ( 7 + l)-(Ci(A--)- ['([(A-s-^s), 

X X Jo X 

where v and f x are evaluated at the point (x,y). 
Since ip € K, and u is defined by (8.4), we have 

aij ,b,ceC(nl\{x = 0}). 
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Combining (8.12) with (5.16), (5.37), (5.45), and (8.14), we obtain that, 
for sufficiently large C depending only on the data, 

1 1 1 1/2 r,+ 

a\\ > -x, a 22 > — , ai2 < tti= x on n 2e- 

6 2c 2 3^/ci 

Thus, 4ana22 — {an) 2 > t£tX on fi^, which implies that equation (8.15) is 

elliptic on and uniformly elliptic on every compact subset of Q^e \ i x = 

Furthermore, using (5.39) and (8.17) and noting A > and x > 0, we 

have 

(8.18) c(x,y)<0 for every (x,y)eSl^ e such that v(x, y) < 0. 

Now we estimate E(x,y). Using (8.14), (5.43), (4.50), and ip G JC, we find 
that, on $7^., 

|^Oi| < C(x + |V>| + PVI + xl^xxl + + I^Vxwl + |£>Vf ) < CMlx, 

\d x 2t5 \ < C{\D^\ + |DVf + IVvV^I + (1 + 1^1)1^1) < CM 1 x 1 /2( 1 + Ml x), 

\d x 6 3A \ < c(i + |y>| + I^Ci(^) I + (i + I^DI^Vl + m\ 2 ) 

<CM 1 (l + M lX ), 

where we have used the fact that |s£[(s)| < C on R. Combining these estimates 
with (8.13)-(8.14), (5.44), and tp e JC, we obtain from (8.13) that 

\E{x,y)\ < CMfx(l + M lX ) < C/C on Q+ £ . 

From this and (7 + 1)^4 > 1, we conclude that the right-hand side of (8.15) is 
strictly negative in f^ £ if C is sufficiently large, depending only on the data. 

We fix C satisfying all the requirements above (thus depending only on 
the data). Then we have 

(8.19) auv xx + ai 2 v xx + a 2 2V yy + bv x + cv < on Sl% £ , 

the equation is elliptic in Cl^e an d uniformly elliptic on compact subsets of 
^2e \ i x = 0}) an d (8.18) holds. Moreover, v satisfies (8.5) and the boundary 
conditions (8.6)-(8.8) and (8.10). Then it follows that 

v > on Q,2 £ - 

Indeed, let zq := (xo,yo) £ ^ 2e De a minimum point of v over and v(zq) < 
0. Then, by (8.6)-(8.7) and (8.10), either zq is an interior point of Q^e or 
zo G ^wedge n {0 < x < 2e}. If zo is an interior point of fi^, then (8.19) is 
violated since (8.19) is elliptic, t>(zo) < 0, and c(zo) < by (8.18). Thus, the 
only possibility is zq G r we ^ e n {0 < x < 2e}, i.e., zq = (xo,0) with xo > 0. 
Then, by (8.2), there exists p > such that B p (z ) n = b p( z o) H {y > 0}. 
Equation (8.19) is uniformly elliptic in B p / 2 (zo) PI {y > 0}, with the coefficients 
aij,b,c e C(B p / 2 (zo) n {y > 0}). Since v(zo) < and u satisfies (8.5), then, 
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reducing p > if necessary, we have v < in B p (zq) l~l {y > 0}. Thus, 
c < on B p (zq) n {y > 0} by (8.18). Moreover, v(x,y) is not a constant 
in B Xq / 2 (xq) PI {y > 0} since its negative minimum is achieved at (xo,0) and 
cannot be achieved in any interior point, as we have showed above. Thus, 
dyv(zo) > by Hopf's Lemma, which contradicts (8.8). Therefore, v > on 
fi£ so that (8.3) holds on Q+ . Then, using (8.9), we obtain (8.3) on fi^. □ 

Now we bound ip x from below. We first prove the following lemma in the 
(£, ^-coordinates. 

Lemma 8.1. If C in (5.16) is sufficiently large, depending only on the 
data, then 

(8.20) ipr, < m 

Proof. We divide the proof into six steps. 
Step 1. Set w = ip v . From tp £ IC and (8.1), 

(8.21) w G c°' a (n+) nc 1 (n+\ r sonic u s ) n c 2 (n + ) . 

In the next steps, we derive the equation and boundary conditions for w 
in S7 + . To achieve this, we use the following facts: 

(i) If C in (5.16) is sufficiently large, then the coefficient An of (5.29) 
satisfies 

-2 tl 

(8.22) \A U (Zty(£, t,),£, rj) \ > > in fi+ 

where C2 and £ are defined in §3.1. Indeed, since c 2 > |£| by (3.5) and (c2,£) — ► 
(c 2 , £) as 0„, -> vr/2 by §3.2, we have c| - £ 2 > 9(c2 - f 2 )/ 10 > if a is small. 
Furthermore, for any (£, 77) G P, we have C2 cos ^ £ ^ £ and thus, assuming 
that cr is small so that |£| < 2|£| and c 2 < 2c~2, we obtain |£| < C. Now, since 
ip G /C, it follows that, if C in (5.16) is sufficiently large, then yl n defined in 
(5.35) with (f> = V implies A n > (c 2 , -£ 2 )/2 on V, and in (5.41) with 
4> = ip implies A 2 n > \t% - £ 2 )/2 onPn{c 2 -r< 4e}. Then (8.22) follows 
from (5.48). 

(ii) Since ip satisfies equation (5.29) in Q + with (8.22), we have 

(8.23) ^ = - 2il2 ^ + i22 ^ in n+, 

An 

where A i:j (£,r)) = A {j (£>Y>(£, 77), £, 77) in f2+. 

5iep 2. We differentiate equation (5.29) with respect to n and substitute 
the right-hand side of (8.23) for tp^ to obtain the following equation for w: 
(8.24) 

A 11 w^+2A 12 w^ ri +A 22 w vri +2(d v A 1 2- dv f 11 Ai2)w^+(d ri A 2 2 — = 0. 

An An 
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By Lemma 5.2, (8.22), and ip G /C, the coefficients of (8.24) are continuous in 
0+ \ r son i c U So, and the equation is uniformly elliptic on compact subsets of 

iL \ L sonic • 

Step 3. By (5.33), we have 

(8.25) w = -v 2 on S := dtt + n {q = -v 2 }. 

Since ip € K, it follows that |£ty(£,fj)| < CM x (c 2 - r) for all (£,77) G 
n+ n {c 2 - r < 2e}. Thus, 

(8.26) w = on r sora j c . 

5iep ^. We derive the boundary condition for ip on T wedge . Then ^ 
satisfies (5.32), which can be written as 

(8.27) - sin 9 W ip% + cos 9 W ip v = on T wedge . 

Since ip £ JC, we have V £ C 2 (f2+ \ r so „j c U S ). Thus we can differentiate 
(8.27) in the direction tangential to T we dge, i-e., apply d T := cos 6 W 9^+sin 9 W d„ 
to (8.27). Differentiating and substituting the right-hand side of (8.23) for ip^, 
we have 
(8.28) 

(cos(20 w ) + ^-sm(29 w ))w£ + isin(20 w )(l + ^ L )w r) = on T wedge . 
An 1 Mi 

This condition is oblique if a is small: Indeed, since the unit normal on ^ wedge 
is (— sin cos 9 W ), we use (3.1) and (8.22) to find 

(cos(2^) + ^sin(2^),isin(2^)(l+4^))-(-sin^,cos^)) > l-Ca > \. 

An 2 An 2 

Step 5. In this step, we derive the condition for w on T s hock- Since ip is 
a solution of (5.29)-(5.33) for <p = ip, the Rankine-Hugoniot conditions hold 
on r shock- Indeed, the continuous matching of ip with (pi — ip 2 across r shock 
holds by (5.21)-(5.23) since (p = ip. Then (4.28) holds and the gradient jump 
condition (4.29) can be written in form (4.42). On the other hand, ip on r shock 
satisfies (5.30) with <p = ip, which is (4.42). Thus, ip satisfies (4.29). 

Since tp £ K, which implies ip £ C 2 (fi + \ T son i c U So), we can differentiate 
(4.29) in the direction tangential to T s h oc k- The unit normal u s on T s hock is 
given by (4.30). Then the vector 

(8.29) Ts = {T S ,T S ) := { , 1 J 

U\ -u 2 Ui - u 2 

is tangential to T s hock- Note that r s / if C in (5.16) is sufficiently large, 
since 

(8.30) \Dip\ < C(a + e) in Q+, \u 2 \ + \v 2 \ < Ca, 
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and u\ > from ip G K, and §3.2. Thus, we can apply the differential operator 
d Ts = T l s d^ + T 2 s d ri to (4.29). 

In the calculation below, we use the notation in §4.2. We have showed in 
§4.2 that condition (4.29) can be written in form (4.33), where F(p, z, u 2 ,v 2 ,£, if) 
is defined by (4.34)-(4.36) and satisfies (4.37). Also, we denote 

(8.31) t(p,u 2 ,v 2 ) = (t 1 ,t 2 )(p,u 2 ,v 2 ) := ( V2+P2 , 1 



ui - u 2 ui- u 2 

»2 



where p = (pi,p 2 ) G R and z G R. Then f G C°°(B s ,(0) x B Ul/50 (0)). Now, 
applying the differential operator d Ts , we obtain that tp satisfies 

(8.32) $(£>fy, Di>, i>, u 2 ,v 2 ,£, V ) = on T shock , 

where 
(8.33) 

2 2 

*(J?,p, z, U2, v 2 , e, r?) = £ r7-;, /»';,• ^(^Pi+^i) for ^ = (^i)ij=i' 

ij=l i=l 

and, in both (8.33) and the calculation below, D^ ly ^F denotes as D^^F, 
(F p .,F z ,F£.) as (F Pj ,F 2 ,F ? J(p,z,U2,v 2 ,£,f7), {t,v) as (f , £)(p, u 2 , v 2 ), and p 
as p(p, z,£,ri), with p(-) and £(•) defined by (4.35) and (4.36), respectively. By 
explicit calculation, we apply (4.34)-(4.36) and (8.31) to obtain that, for every 
(p,z,u 2 ,v 2 ,£,rj), 

2 

(8.34) fl ^ + F^) = (pi-p)r-u = 0. 

i=l 

We note that (4.28) holds on T s h oc k- Using (8.32) and (8.34) and express- 
ing £ from (4.28), we see that ijj satisfies 

(8.35) $(D 2 il>,Dil>,ip,U2,V2,ri) = onT shock , 
where 

2 

(8.36) ®(R,p,z,u 2 ,v 2 ,r)) = ^ f l ^ Pj {p,z,u 2 ,v 2 ,rf)R ij , 

* is defined by (4.39) and satisfies * G C°°(A) with 11*11^*^ depending 
only on the data and k G N, and .4 = B s *(0) x (— 5*,<5*) x -B Ul / 50 (0) x 
(-6c 2 /5,6c 2 /5). 

Now, from (4.34)^(4.36), (4.39), and (8.31), we find 

f((0,0),0,0) = (0, 1), D p *((0,0), 0,0, 0,17) = (p 2 (c|-f ), (^^i-p' 2 £»- 

Ul 
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Thus, by (8.36), we obtain that, on R 2x2 x A, 
(8.37) $(R,p,z,u 2 ,v 2 ,v) 



2 

■92 ~ Pi 



= P2( c 2 ~ i 2 )R21 + (— - p' 2 £)vR22 + ^2 Eij(p,Z,U2,V 2 ,ri)Rij, 

U\ . . 

1.3=1 

where E {] - £ C°°(A) and 

\E ij (p,z,u 2 ,v 2 ,r))\ < C(\p\ + \z\ + \u 2 \ + \v 2 \) for any (p, z,u 2 ,v 2 ,rj) £ A, 

with C depending only on ||Z? 2 *I'|| C , 0( .^. 

From now on, we fix (u 2 ,v 2 ) to be equal to the velocity of state (2) obtained 
in §3.2 and write Eij(p,z,r]) for Eij(p, z,u 2 ,v 2 ,rj). Then, from (8.35) and 
(8.37), we conclude that Y> satisfies 
(8.38) 

2 

,P2 ~ Pi 



P2( c 2-| 2 )V^r,+ ( P^V^vv+^^Z E ij{ D 4>^,v) D ij^ = on T shock , 

Ul . . 

and Eij = Eij(p, z, r]),i,j = 1, 2, are smooth on 



B := B s ,(0) x (S*,S*) x (-6c 2 /5, 6c 2 /5) 

and satisfy (4.43) with C depending only on the data. Note that 

(Dip(£,rf), ip{€,rf), rf) 66 on T shock , 

since ip G /C and (5.16) holds with sufficiently large (7. Expressing from 
(8.23) and using (8.22), we can rewrite (8.38) in the form 

(P2(4 - i 2 ) + £i(Zty,^M))^„ + ((^^ - p' 2 £> + E 2 (D^,4>,r,))i> m = 

on T s h oc k, where the functions Ei = Ei(p, z,rf),i = 1,2, are smooth on B and 

satisfy (4.43). Thus, w satisfies 

(8.39) 

(P 2 (C2 " f) + E 1 (Zty, ^, »/)) ^ + ((^^ - p 2 |)r? + E 2 (D^, ^ rf)) w v = 

on r s /j OC fc. Condition (8.39) is oblique if C is sufficiently large in (5.16). Indeed, 
we have c 2 > yjjc 2 , which implies — |£| 2 > c 2 C2 ~^ > by using (4.8). Now, 
combining (4.30) and (4.43) with ij) G K, and (3.24), we find that, on T s hock-> 



(P2(4 - e 2 ) + ^, »/), (— - - p' 2 e> + e 2 (d^, ^,r,))- v s 



> ftc 2 ^j@ - C(M l£ + M 2 a) > P 2 c 2 ^-g^. > o. 
Also, the coefficients of (8.39) are continuous with respect to (£,77) £ r s ^ oc fc. 
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Step 6. Both the regularity of w in (8.21) and the fact that w satisfies 
equation (8.24) that is uniformly elliptic on compact subsets of 0+ \ T son i c 
imply that the maximum of w cannot be achieved in the interior of 17 + , unless 
w is constant on by the Strong Maximum Principle. Since w satisfies 
the oblique derivative conditions (8.28) and (8.39) on the straight segment 
Fwedge and on the curve T s h oc k that is C 2 ' a in its relative interior, and since 
equation (8.24) is uniformly elliptic in a neighborhood of any point from the 
relative interiors of T wec ig e and r<,/j OC /%, it follows from Hopf's Lemma that the 
maximum of w cannot be achieved in the relative interiors of T we d ge and T s h oc k, 
unless w is constant on Q + . Now conditions (8.25)-(8.26) imply that w < 
on VL + . This completes the proof. □ 

Using Lemma 8.1 and working in the (x, y)-coordinates, we have 

Proposition 8.2. If C in (5.16) is sufficiently large, depending only on 
the data, then 

(8.40) V* > - 3(7 4 +1) ^ in V + n {x < As}. 

Proof. By definition of the (x, y)-coordinates in (4.47), we have 

/o a-,s , ■ n , cos# , 

(8.41) ip v = -sm 9 ip x -\ — ipy, 

where (r, 9) are the polar coordinates in the (£, r/)-plane. 

From (7.20), it follows that, for sufficiently small a and e, depending only 
on the data, 

7]>7]* for all (£, n) £ V n {c 2 - r < 4e}, 

where {l(r)*),V*) is the unique intersection point of the segment {(l(v),v) '■ 
r\ e (0,771]} with the circle <9-B C2 _4 £ (0). Let f)* be the corresponding point 
for the case of normal reflection, i.e., f)* = \J (c 2 — 4e) 2 — £ 2 . By (3.5), f)* > 



c\-i 2 l2 > if £ is sufficiently small. Also, from (4.3)-(4.4) and (3.24), 

and using the convergence (9 S ,C2,£) — > (7r/2, c 2 ,£) as 9 W — > vr/2, we obtain 
rf > fj* /2 and c 2 < 2c 2 if cr and e are sufficiently small. Thus, we conclude 
that, if C in (5.16) is sufficiently large depending only on the data, then, for 
every (£, 77) € X> n {c 2 — r < 4e}, the polar angle satisfies 



(8.42) sin(9 = ^=2^ > 0, |cot0| 

V? 2 + ?7 2 



- 852 <C 

-e 



From (8.41)-(8.42) and Lemma 8.1, we find that, on f2+ n {c 2 - r < 4e}, 

/ . , 1 cot 9 , cot 9 , , . 

8.43 Vx = — + V»> ^i/>-C^i,- 
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Note that ip G K, implies \ip y (x,y)\ < Mix 3/2 for all (x,y) G J7 + n{c 2 -r < 2e}. 
Then, using (8.43) and (5.16) and choosing large C, we have 

ip x > -r^- — -x in tt + n {x < 2e}. 

Vx ~ 3( 7 + 1) 1 - / 

Also, tp <E K. implies 

IV-xl < M 2 <7 < 3(7 4 + ^ (2e) on n + D{2e<x< 4e}, 

where the second inequality holds by (5.16) if C is sufficiently large depending 
only on the data. Thus, (8.40) holds on Q^. □ 



9. Proof of Main Theorem 

Let C be sufficiently large to satisfy the conditions in Propositions 7.1 and 
8.1-8.2. Then, by Proposition 7.1, there exist <7o,e > and Mi,M 2 > 1 such 
that, for any a G (0, oo], there exists a solution -0 G /C(cr, e, Mi, M2) of problem 
(5.29)-(5.33) with cp = ip. Fix <r G (0, do] and the corresponding "fixed point" 
solution ip, which, by Propositions 8.1-8.2, satisfies 

4 

\il>x\ < ^7 — in n + n {x < 4e}. 
3(7 + 1) 

Then, by Lemma 5.4, ip satisfies equation (4.19) in Q, + (ip). Moreover, ip sat- 
isfies properties (i)-(v) in Step 10 of §5.6 by following the argument in Step 
10 of §5.6. Then, extending the function ip = ip + <p2 from f2 := £l + (ip) to the 
whole domain A by using (1.20) to define ip in A \ 0, we obtain 

<p g w£> (A) n {u* =0 c\K u5)n c^{K)) , 

where the domains Aj, i = 0,1,2, are defined in Step 10 of §5.6. From the 
argument in Step 10 of §5.6, it follows that ip is a weak solution of Problem 2, 
provided that the reflected shock S± = P0-P1-P2 H A is a C 2 -curve. 

Thus, it remains to show that Si = Po-Pi-fbnA is a C 2 -curve. By definition 
of p and since ip G /C(<r, e, Mi, M2), the reflected shock Si = PqP\P2 n A is 
given by Si = {£ = fs 1 {v) ■ Vp 2 < V < VPo}, where t]p 2 = -v 2 , f?p = 
iaS^y>0, and 

[ ' ISlW \l(ri) ii V e(r) Pl ,r, Po ), 

where l{rj) is defined by (4.3), i]p 1 = 771 > is defined by (4.6), and r/p > rjp 1 
if a is sufficiently small, which follows from the explicit expression of ijp given 
above and the fact that (0 S , C2, £) — > (vr/2, c~2, £) as — ► it/2. The function 
is defined by (5.21) for <p = ip. 
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Thus we need to show that fs 1 G C 2 ([r]p 2 ,r]p a ]). By (4.3) and (5.24), it 
suffices to show that fs 1 is twice differentiable at the points rjp 1 and r]p 2 . 

First, we consider fg 1 near r]p 1 . We change the coordinates to the (x,y)~ 
coordinates in (4.47). Then, for sufficiently small e\ > 0, the curve = 
fSriv)} n {c 2 - £\ < r < c 2 + £1} has the form {y = fs 1 (x) : -e 1 < x < £i}, 
where 



(9-2) f Sl (x) 



f^(x) if x G (0,£i), 
/o(x) if x € (-£i,0), 



with /o and defined by (5.9) and (5.25) for <f> = ip. In order to show that fs 1 
is twice differentiable at rjp 1 , it suffices to show that fs 1 is twice differentiable 
at x = 0. 

From (5.26)-(5.27) and (5.9), it follows that / Sl G C 1 ((-£i, £i)). More- 
over, from (5.3), (5.6), (5.22), and (5.27), we write ifi,<p 2 , and Y> i n the (x,y)- 
coordinates to obtain that 



(9-3) f' Sl (x) 



d y ^-V2-ri [XjsAX)) ifXG(0 ' £l) ' 

cUpi - ip 2 ) 



and that /q(x) is given for x G (— £i,£i) by the second line of the right-hand 
side of (9.3). Using (5.3) and ip G /C with (5.16) for sufficiently large (7, we 
have 

(9.4) \f' Si (x) - ft(x)\ < C\D (x ^(x,U(x))\ for all x G (0, £l ). 

Since ip satisfies (5.30) with <p = ip, it follows that, in the (x, y)-coordinates, 
ip satisfies (6.6) on {y = f^(x) : x G (0,£i)}, and (6.8) holds. Then it follows 
that 

\il>x{x,Uix))\ < C(W y (x,U(x))\ + \i;(x,U(x))\) < Cx*l\ 

where the last inequality follows from ip G /C. Combining this with (9.2), (9.4), 
and f Sl ,fo G C 1 ((-£]_,£])) yields 

\f' Si (x) - ft(x)\ < Cx 3 / 2 for all x G (-£i,£i). 

Then it follows that f' Sl {x) — fo(x) is differentiable at x = 0. Since /o G 
C ,oc ((— £i, £i)), we conclude that fs 1 is twice differentiable at x = 0. Thus, /gj 
is twice differentiable at rjp 1 . 

In order to prove the C 2 -smoothness of fs ± up to r]p 2 = —v 2 , we extend 
the solution <p and the free boundary function fs 1 into {rj < —v 2 } by the even 
reflection about the line Eo C {ij = — ^2} so that P 2 becomes an interior point 
of the shock curve. Note that we continue to work in the shifted coordinates 
defined in §4.1, that is, for (£, 77) such that r\ < —v 2 and (£, — 2v 2 — rj) G f2+(^), 
we define (</?, y>i)(£, rf) = (ip,ip 1 )(^,-2v 2 - rj) and fs.iv) = ~^v 2 - rj for 991 
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given by (4.15). Denote ^+(P 2 ) := B £l (P 2 ) n > fsAv)} for sufficiently 
small ei > 0. From y> G C 1 > a (fi+(^)) n C 2 ' a (^+(V')) and (4.13), we have 

^c 1 - a («)nc 2 - a (fl+(P 2 )). 

Also, the extended function ip\ is in fact given by (4.15). Furthermore, from 
(5.20) and (5.22), we can see that the same is true for the extended functions 
and hence 

{£ > fsM}nB £l (P 2 ) = W< p 1 }nB £l (P 2 ), f Sl e c 1 ' a ((-v 2 - £ j-,-v 2 + £ j-)). 

Furthermore, from (1.8)-(1.9) and (4.13), it follows that the extended (p sat- 
isfies equation (1.8) with (1.9) in fi+(P2), where we have used the form of 
equation, i.e., the fact that there is no explicit dependence on (£,77) in the 
coefficients and that the dependence of Dip is only through \D<p\. Finally, 
the boundary conditions (4.9) and (4.10) are satisfied on r £l (P2) := {£ = 
fSiiv)} n B £l {P 2 ). Equation (1.8) is uniformly elliptic in JI+ (P2) for if, which 
follows from 92 = (p 2 +ip and Lemmas 5.2 and 5.4. Condition (4.10) is uniformly 
oblique on r £l (P 2 ) for 92, which follows from §4.2. 

Next, we rewrite equation (1.8) in fi+(P2) and the boundary conditions 
(4.9)-(4.10) on r El (P2) in terms of u := <p\ — (p. Substituting u + ipi for ip into 
(1.8) and (4.10), we obtain that u satisfies 

F(D 2 u,Du,u,£,ri) =0 in ft+(P 2 ), u = G(Du,u,£,r]) = onr £l (P 2 ), 

where the equation is quasilinear and uniformly elliptic, the second boundary 
condition is oblique, and the functions F and G are smooth. Also, from (5.20) 
which holds for the even extensions as well, we find that d^u > on T £l (P2). 
Then, applying the hodograph transform of [28, §3], i.e., changing (£,77) — ► 
(X, Y) = (u(£, 77), r/), and denoting the inverse transform by (X, Y) — > (£, 77) = 
(v(X,Y),Y), we obtain 

v e C^(Bl((o,-v 2 ))) n C 2 '«(P+((0, -t; 2 ))), 

where Bf((0, -v 2 )) := B s ((0, -v 2 ))n{X > 0} for small 5 > 0, u(X,y) satisfies 
a uniformly elliptic quasilinear equation 

F(D 2 v,Dv,v,X,Y) =0 inP+((0,-7j 2 )) 

and the oblique derivative condition 

G(Dv, v, Y) = on 8Bf((0, -v 2 )) n {X = 0}, 

and the functions P and G are smooth. Then, from th e local estimat es near 
the boundary in the proof of [31, Theorem 2], v G C 2,a (B^ 2 ((0, —v 2 ))). Since 
fSiiv) = u (0, 77), it follows that fs 1 is C 2,a near 77p 2 = — v 2 . 

It remains to prove the convergence of the solutions to the normal reflec- 
tion solution as 6 W — ► 7r/2. Let 9 % w — >7r/2asi — > 00. Denote by <p l and f l 
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the corresponding solution and the free-boundary function respectively, i.e., 
PqP\P 2 PI A for each i is given by {£ = f l (rf) : rj € (r/p 2 , ??p )}. Denote by <^°° 
and f°°(rj) = £ the solution and the reflected shock for the normal reflection 
respectively. For each i, we find that tp 1 — ip 2 = ^ i n the subsonic domain flf, 
where ip l is the corresponding "fixed point solution" from Proposition 7.1 and 
Y>* € /C(vr/2 - el^^MlMl) with (5.16). Moreover, /* satisfies (5.24). We 
also use the convergence of state (2) to the corresponding state of the normal 
reflection obtained in §3.2. Then we conclude that, for a subsequence, f l — ► f°° 
in and ip l — > t^ 00 in C 1 on compact subsets of {£ > £} and {£ < £}. Also, 
we obtain \\(D(p l , (p l )\\ (K) < C(-^0 f° r every compact set K. Then ip l — > 9900 
in W^'^A) by the Dominated Convergence Theorem. Since such a converging 
subsequence can be extracted from every sequence l w — > 7r/2, it follows that 
— > Voo as — > 7r/2. 

A. Appendix: Estimates of Solutions to Elliptic Equations 

In this appendix, we make some careful estimates of solutions of bound- 
ary value problems for elliptic equations in R 2 , which are applied in §6-§7. 
Throughout the appendix, we denote by (x, y) or (X, Y) the coordinates in R 2 , 
by := {y > 0}, and, for z = (x,Q) and r > 0, denote B+(z) := B r (z) nR 2 + 
and Y. r {z) := B r (z) n {y = 0}. We also denote B r := B r (0), B+ := B+(0), 
and T, r := S r (0). 

We consider an elliptic equation of the form 

(A.l) A n u xx + 2A 12 u xy + A 2 2U yy + Aiu x + A 2 u y = f, 

where Aij = Aij(Du,x,y), A4 = Ai(Du, x,y), and / = f(x,y). We study 
the following three types of boundary conditions: (i) the Dirichlet condition, 
(ii) the oblique derivative condition, (iii) the "almost tangential derivative" 
condition. 

One of the new ingredients in our estimates below is that we do not assume 
that the equation satisfies the "natural structure conditions" , which are used in 
the earlier related results; see, e.g., [20, Chapter 15] for the interior estimates 
for the Dirichlet problem and [37] for the oblique derivative problem. For 
equation (A.l), the natural structure conditions include the requirement that 
|p| I .Dp^ij I < C for all p £ R 2 . Note that equations (5.42) and (5.49) do not 
satisfy this condition because of the term xQ\{—) in the coefficient of ip xx . 
Thus we have to derive the estimates for the equations without the "natural 
structure conditions" . We consider only the two-dimensional case here. 

The main point at which the "natural structure conditions" are needed 
is the gradient estimates. The interior gradient estimates and global gradient 
estimates for the Dirichlet problem, without requiring the natural structure 
conditions, were obtained in the earlier results in the two-dimensional case; 
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see Trudinger [47] and references therein. However, it is not clear how this 
approach can be extended to the oblique and "almost tangential" derivative 
problems. We also note a related result by Lieberman [34] for fully nonlin- 
ear equations and the boundary conditions without obliqueness assumption in 
the two-dimensional case, in which the Holder estimates for the gradient of a 
solution depend on both the bounds of the solution and its gradient. 

In this appendix, we present the C 2 ' a -estimates of the solution only in 
terms of its C-norm. For simplicity, we restrict to the case of quasilinear 
equation (A.l) and linear boundary conditions, which is the case for the ap- 
plications in this paper. Below, we first present the interior estimate in the 
form that is used in the other parts of this paper. Then we give a proof of 
the C 2 ' Q -estimates for the "almost tangential" derivative problem. Since the 
proofs for the Dirichlet and oblique derivative problems are similar to that for 
the "almost tangential" derivative problem, we just sketch these proofs. 

Theorem A.l. Let u G C 2 (B 2 ) be a solution of equation (A.l) in B2. 
Let Aij(p,x,y), Ai(p,x,y), and f(x,y) satisfy that there exist constants A > 
and a G (0, 1) such that 

n 

(A.2) X\fi\ 2 < ^ AijViVj < A'Vl 2 for all (x,y) G B 2 , p,n& R 2 , 

(A. 3) \\{Aij, ^i)|lca( R 2 X ^) + \\D p (Aij, ^4j)Hc( R 2 X ^) + Il/Ilc«(sl) - ^ 1 - 

Assume that \\u\\c(b^) ^ M. Then there exists C > depending only on (A, M) 
such that 

(A.4) \H\ C ^ m < C (\\ u \\c(B- 2 ) + \\f\\c°(B- 2 )^ 

Proof. We use the standard interior Holder seminorms and norms as 
defined in [20, Eqs. (4.17), (6.10)]. By [20, Theorem 12.4], there exists G 
(0, 1) depending only on A such that 

MUb, < C(A)(||u|| , B2 + ||/ - A 1 D 1 u - A 2 D 2 u\\V> Ba ) 
<C(X,M)(l + \\f\\^ 2 + \\Du\\^ 2 ). 

Then, applying the interpolation inequality [20, (6.82)] with the argument 
similar to that for the proof of [20, Theorem 12.4], we obtain 

\\u\\* 1AB2 <C(\,M)(l + \\f\\^ 2 ). 
Now we consider (A.l) as a linear elliptic equation 

n n 

^ aij(x)u XiXj + ^ j a i (x)u Xi = f(x) in B 3/2 

i,j=l i=l 
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with coefficients a«j(x) = Aij(Du(x),x) and at = Ai(Du(x),x) in C^(B 3 / 2 ) 
satisfying 

11(^^)11^(1^) < C(A,M). 

We can assume (3 < a. Then the local estimates for linear elliptic equations 
yield 

Mc^iB^l) ^ C (\M)(\\ U \\ c(E7J - 2) + ll/Hc^)). 

With this estimate, we have ||(ai.j,ai)||c<a(^I) — C(X,M). Then the local 
estimates for linear elliptic equations in B 5 / 4 yield (A. 4). □ 

Now we make the estimates for the "almost tangential derivative" prob- 
lem. 

Theorem A. 2. Let A > 0, a e (0, 1), and e > 0. Zei $ e C 2 ' a (R) satisfy 
(A.5) ||*||c*.«(R) < A" 1 , 

and denote n+ := 5 R n {y > e$(x)} for R > 0. Ze£ u G C 2 (B+) n C^Kf) 
satisfy (A.l) m fij and 

(A. 6) u x = eb(x,y)u y + c(x,y)u on r$ := i?2 n {y = e$(x)}. 

Lei Aij(p,x,y), Ai(p,x,y), a(x,y), b(x,y), and f(x,y) satisfy that there exist 
constants A > and a £ (0, 1) such that 

n 

(A.7) X\fi\ 2 < ^2 Aijfiifij < A -1 !^! 2 f° r ( x ,y) G ^> P>A* G r2 ' 

(A.8) ^)ll c ^ xR2) + IIA>(^, ^)II C( ^ XR2) + ll/ll^ < A" 1 , 

(A.9) ll(^' c )llcfi,a(n+) — A _1 - 

Assume that ||-u|| ,— + < M. Then there exist £q(X, M, a) > andC(X, M, a) > 
smc/i #iai, if e G (0,£o), we /iave 

(A.10) ll^llc-(^)^ C (^llc(^) + ll/Hc=(^))- 

To prove this theorem, we first flatten the boundary part T$ by defining 
the variables (X,Y) = ^>(x,y) with (X,Y) = (x,y — e&(x)). Then = 
* _1 (X,y) = (X,y + From (A.5), we have 

(A.H) ||* - Id\\ c2 ^ } + II*- 1 - Id\\ C2 , a(5f) < eX-\ 

Then, for sufficiently small e depending only on A, the transformed domain 

:= satisfies 
(A.12) 

5 2 + -2 £ /a c ^ + c 5+ , t>2 c r+ := {y > o}, a2?+n{y = o} = *(r*) ; 
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the function 

v(X, Y) = u(x, y) := u(* _1 (X, Y)) 

satisfies an equation of form (A.l) in with (A.7)-(A.8) and the correspond- 
ing elliptic constants A/2; and the boundary condition for v by an explicit 
calculation is 
(A.13) 

Vx = e(6(^" 1 (X,0)) + &(X))vy + c(^~ 1 (X,0))v on V+ n {Y = 0}, 

i.e., it is of form (A. 6) with (A. 9) satisfied on with elliptic constant A/4. 
Moreover, by (A.11)-(A.12), it suffices for this theorem to show the following 
estimate for v(X,Y): 

(A.14) ||t,|| 2>a>flj>5 < C(\,M,a)(\\v\\ , Bt2s/x + \\f\\ a , Bt2 J. 

That is, we can consider the equation in -B^L 2 e/A an< ^ condition (A.13) on 
S 2 _2e/A or ) by rescaling, we can simply consider our equation in B£ and con- 
dition (A.13) on S2 := B2(~\{Y = 0}. In other words, without loss of generality, 
we can assume <E> = in the original problem. 

For simplicity, we use the original notation (x , y , u(x , y)) to replace the 
notation (X,Y,v(X,Y)). Then we assume that <3? = 0. Thus, equation (A.l) 
is satisfied in the domain B^ ', the boundary condition (A. 6) is prescribed on 
£ 2 = B 2 n {y = 0}, and conditions (A.7)-(A.9) hold in B% . Also, we use 
the partially interior norms [20, Eq. 4.29] in the domain B^ U £2 with the 
related distance function d z = dist(z, dB^ \ £2)- The universal constant C in 
the argument below depends only on A and M, unless otherwise specified. 

As in [20, §13.2], we introduce the functions Wi = Diu for i = 1,2. Then 
we conclude from equation (A.l) that w\ and u>2 are weak solutions of the 
following equations of divergence form: 



^ A22 A22 A22 A22 A22 



J_ 


_A 


A22 


a 22 


f 


A! 




An 



(A.16) D 11 W2+D 2 (^D 1 W2+^D 2 W2) = D 2 (4 ^D x u-^D 2 u). 

A\\ An An An An 

From (A. 6), we have 
(A. 17) w\ = g on £2, 

where 

(A. 18) g:=ebw 2 + cu for B%. 



We first obtain the following Holder estimates of Diu. 
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Lemma A.l. There exist (3 G (0, a] and C > depending only on A such 
that, for any zq G B^ U X2; 
(A.19) 

d i[ w i]oAB dzo/1B (zo)nB+ < C (\\( Du J)\\o,o,B dzo/2 (z )nB+ + d z [9}oAB dzo/2 (z )nB+)- 

Proof. We first prove that, for z\ G £2 and B^ R {z\) C Bg", 
(A.20) j R /3 [^i] 0A b+( 2i ) < C^{Du,Rf)\\ 0iQtB}R{zi) + R f3 [9] 0A Bl R (z 1 ))- 
We rescale u, w\, and / in B^izi) by defining 

(A.21) «(Z) = \u{ Zl + 222Z), /(Z) = 2Rf( Zl + 2i?Z) for Z G B+, 
2il 

and Wi = DziU. Then w\ satisfies an equation of form (A. 15) in Bf with u 
replaced by u whose coefficients Aij and Ai satisfy (A.7)-(A.8) with unchanged 
constants (this holds for (A. 8) since R < 1). Then, by the elliptic version of [36, 
Theorem 6.33] stated in the parabolic setting (it can also be obtained by using 
[36, Lemma 4.6] instead of [20, Lemma 8.23] in the proofs of [20, Theorem 8.27, 
8.29] to achieve a = ao in [20, Theorem 8.29]), we find constants /3(A) G (0, 1) 
and C(A) such that 

[wi] 0AB + /2 < C(\\(Du,f)\\ 00B + + [li*l]o, A B in { F 0}) 

for (3 = min(/3,a). Rescaling back and using (A.17), we have (A.20). 

If z\ G B2 and B2r{z\) C B^ , then an argument similar to the proof of 
(A.20) by using the interior estimates [20, Theorem 8.24] yields 

(A.22) i?>iW R (*0 ^ C\\(Du,Rf)\\ 0fiMzi) . 

Now let zq = (xo,yo) G B^ U £2- When yo < d Za /8, then, denoting 
Zq = (xq , 0) and noting that d z > o > d Zo , it is easy to check that 

B dzJ16 (z ) n B+ C B+ o/8 (z' ) C B+, B+ J8 (z' ) C B dzj2 {z Q ) n B+ 

and then applying (A.20) with z\ = z' Q and R = d Za /8 < 1 and using the 
inclusions stated above yield (A.19). When yo > d Zo /8, B dz i^{zq) C 5^". 
Then applying (A.22) with zi = z and i? = d z J16 < 1 yields (A.19). □ 

Next, we make the Holder estimates for Du. We first note that, by (A. 9) 
and (A. 18), g satisfies 

(A.23) \Dg\<C(e\D 2 u\ + \Du\ + \u\) in B+ 

(A.24) \g] 0> i3,B dx/a (z)nB+ ^ C ( £ l Du \o,f3,B dz/2 (z)nBt + ll«lli,o,B <wa (*)nB+) • 

Lemma A. 2. Let (3 be as in Lemma A.l. Then there exist £o(A) > and 
C(A) > such that, ifO<e<eo, 

d z [ Du \o,/3,B dzo/32 (z )nB+ ^ C (\\ u \\l,0,B dzo/2 (z )nB+ + ed z [ Du \o,f3,B d , o/2 (z )nB+ 

(A- 25 ) +H/llo,OA, o/2 (,o)nB 2 +) 
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for any zq E B^ U T, 2 . 

Proof. The Holder norm of D\u has been estimated in Lemma A.l. It 
remains to estimate D211. We follow the proof of [20, Theorem 13.1]. 

Fix zq G By U £ 2 . In order to prove (A. 25), it suffices to show that, for 
every z G n -B+ and every R > such that C B^/ieC^o), 

we have 

(A.26) / \D 2 u\ 2 dz < 4ji? 2/3 , 

JB R (z)nB+ df 

where L is the right-hand side of (A. 25) (cf. [20, Theorem 7.19] and [36, 
Lemma 4.11]). 

In order to prove (A.26), we consider separately case (i) B 2 r(z) n S 2 / 
and case (ii) B 2 r(z) n £2 = 0- 

We first consider case (i). Let B 2 r(z)C\Y> 2 7^ 0- Since Br{z) C /32(^o)> 
then B 2 r(z) C /i6(zo) so that 

(A.27) 2i?<d 20 . 

Let ?? G Cg(5 2 fl(£)) and C = ?? 2 (^i - £?)• Note that ( G W 1,2 (S 2 fl(£) n 5+) by 
(A. 17). We use ( as a test function in the weak form of (A. 15): 

2 2 

(A.28) f V] AijDiWiDj(dz = [ -7- ( + f)Di(dz, 

Jb+A 2 2^ 1 Jb+A 22 ^ 

and apply (A.7)-(A.8) and (A.23) to obtain 
(A.29) 

y + \D Wl \ 2 ri 2 dz < C j + (^(5 + e)\D Wl \ 2 + e\D 2 u\ 2 ^jn 2 

+(] + 1) (d^l 2 + \f\ii 2 ){w! - g) 2 + {\Du\ 2 + M V) 

where C depends only on A, and the sufficiently small constant 5 > will be 
chosen below. Since 

(A.30) |Z^!| 2 = (D llU ) 2 + (D 12 u) 2 , 

it remains to estimate |-£?22^| 2 - Using the ellipticity property (A. 7), we can 
express D 22 u from equation (A.l) to obtain 

/ \D 22 u\ 2 r] 2 dz<C(X) [ (\D u u\ 2 + \D 12 u\ 2 + \Du\ 2 )r] 2 dz. 
Jb+ Jb} 
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Combining this with (A.29)-(A.30) and using (A. 8) to estimate |/| yield 
(A.31) 

/ \D 2 u\ 2 r) 2 dz<C [ ((e + 6)\D 2 u\ 2 rj 2 

JB+ JB+ V 

+(i + 1) ((\D V \ 2 + V 2 )( Wl - gf + (\Du\ 2 + M V) )dz. 

Choose Eq = 5 = (AC)' 1 . Then, when e G (0, £o)> we have 
(A.32) 

J \D 2 u\ 2 ifdz < C J ((\Dr]\ 2 + r ] 2 )(w 1 - g) 2 + (\Du\ 2 + |u| 2 )?? 2 ) dz. 

Now we make a more specific choice of r\: In addition to rj G Cq(B 2 r(z)) : 
we assume that rj = 1 on Br(z), < 77 < 1 on R 2 , and jDryj < 10/ i?. 
Also, since B 2 r(z) Pi E 2 / 0, then, for any fixed z* G B 2 r(z) Pi £ 2 , we have 
|z-z*| < 2i? for any z G B 2R (z). Moreover, (wi - #) (z* ) = by (A.17). Then, 
since B 2R (z) C 5^ o / 16 (z ), we find from (A.19), (A.24), and (A.27) that, for 
any z G -B 2fi (z) Pi , 

|K - 5 )(z)| = |(wi - 5 )(z) - (wi - g)(z*)\ < \ Wl (z) - Wl (z*)\ + \g(z) - g(z*)\ 
C 

< -w{\\(Du,f)\\ 0>0!Bdzo/2{zo)nB} +di\g] 0AB MnB +)\z - z*f 



+ b]o,/3,B^ o/2 ( Zo )niJ : + l' : ' - - 



^ C {\\\( Du if)\\vfl,B dzo/2 { Zo )nBi + £ [Du] 0ABdzo/2{zo)nB} 

+ \\ u \\o,o,B dzo/2 (z )nB+) Rl3 ■ 
Using this estimate and our choice of rj, we obtain from (A.32) that 

\D 2 u\ 2 dz 

B R (z)nB+ 

< C(-^\\(Du,f)\\ 2 R . , nR+ +e 2 [Du} 2 R , , nR+ )i? 2/3 

— V ,2/3 1 1 v 'JJ»o,0,B d /2 (z )nBt l J0,/3,B d;So /2(2o)n.BJ/ 

+^H-lllo,^ o/2 (, )nB 2+ (^ /3 + ^ 2 ), 

which implies (A. 26) for case (i). 

Now we consider case (ii): z G B^ and R > satisfy Br(z) C -B^ ^(-^o) 
and 5 2iJ (z) n S 2 = 0. Then B 2R (z) C ^ q/16 (z ) D S+. Let r/ G (^(^(i)) 

and C = r? 2 (^i - u>i(z)). Note that ( G Wq' 2 (B£) since B 2 r(z) C £+. Thus 
we can use £ as a test function in (A. 28). Performing the estimates similar to 
those that have been done to obtain (A.32), we have 

(A.33) / \D 2 u\ 2 ifdz <C{\) f ({\Dr]\ 2 + r] 2 )(wi-wi(z)) 2 + \Du\ 2 r] 2 )dz. 

JBt JBt 
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Choose 7] G Cq (B 2 r(z)) so that 7] = 1 on Br(z), < ?j < 1 on R 2 , and 
|Dr?| < 10/R. Note that, for any z G B 2R (z), 

\ Wl (z) - Wl (z)\ < C( w \\(Du,f)\\ QABdzQ/2[zo)nB} +e[Du] 0ABdzo/2(zo)nB +)Rf 5 

by (A.19) since B 2R (z) C B d ^ /w (z )nB^ . Now we obtain (A.26) from (A.33) 
similar to that for case (i). Then Lemma A. 2 is proved. □ 

Lemma A. 3. Let (3 and Eq be as in Lemma A. 2. Then, for e G (0,£o), 
there exists C(A) such that 

(A.34) M^b+ue, - C ^ Wus 2 + ^MUb+ue, + H/llo,o, B +), 

where [•]* and || • ||* denote the standard partially interior seminorms and norms 
[20, Eq. 4.29]. 

Proof. Estimate (A.34) follows directly from Lemma A. 2 and an argument 
similar to the proof of [20, Theorem 4.8]. Let z±,z 2 G B^ with d Zl < d Z2 (thus 
d Zl ,z 2 = d Zl ) and let \z\- z 2 \ < d Zl /64. Then z 2 G B dz j 32 (z ) n -B^ and, by 
Lemma A. 2 applied to zq = z±, we find 

+ ll/llo,o,B ti , l/2 (2i)n J B 2 + ) 

^(IMU.O.B+UE, + £ K,AB 2 +UE 2 + ll/llo,0,B 2 + )' 

where the last inequality holds since 2d z > d Zl for all z G B dz ^j 2 {z\) n £?2~. If 
zi,Z2 £ ^2" w ith 4i < 4 2 and \z\ — z 2 \ > d 2l /64, then 

g ^yf^ < 64(4,1^01+4.1^2)1) < 64||«||; >0> ^ uEa . 
This completes the proof. □ 

Now we can complete the proof of Theorem A. 2. For sufficiently small 
Eq > depending only on A, when e G (0, £0), we use Lemma A. 3 to obtain 

(A.35) MUb 2 +us 2 < C ( A )(IMIWue 2 + Wf\\o,o,B})- 

We use the interpolation inequality [20, Eq. (6.89)] to estimate 

IMIi,o,b 2 + us 2 < C((3,S)\\u\\^ B + + S[u}* 1AB + u x 2 
for 5 > 0. Since (3 = /3(A), we choose sufficiently small 5(A) > to find 
(A.36) IMI^b+ue, < C(A)(hlloA^ + II/Ho,o,b 2 + ) 
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from (A. 35). In particular, we obtain a global estimate in a smaller half-ball: 

(A.37) < C(A)(IMIo,o,b 2 + + imio,o,£ 2 +)- 

We can assume (3 < a. Now we consider (A. 15) as a linear elliptic equation 

2 

(A.38) D i {a ij {x,y)D j w l ) = D\F in B+ /5 , 

where aij(x,y) = (Aij/A 22 )(Du(x,y),x,y) for i+j < 4, 022 = 1, and F(x,y) = 
(A 1 D 1 u + A 2 D 2 u + f)/A 22 with (Ai^Ai) = (Aij,Ai)(Du(x,y),x,y). Then 
(A. 36), combined with (A. 8), implies 

(A.39) \M\oab+ /b <C(X,M). 

From now on, d z denotes the distance related to the partially interior norms in 
B^ 5 U £9/5, i.e., for z G B^ 5 , d z := dist(z, dB^ 5 \ £9/5). Now, similar to the 
proof of Lemma A.l, we rescale equation (A.38) and the Dirichlet condition 
(A.17) from the balls B^(z[) C £+ /5 and B R (zi) C 5+ /5 with R < 1 to 
B = Bf or B = Bi, respectively, by defining 

(w 1 ,g,c Hj )(Z) = (w 1 ,g,a ij )(z 1 +RZ), F(Z) = RF( Zl + RZ) for Z G B. 

Then Yli j=i Di(&ij( x i v)DjW\) = DiF in f?, the ellipticity of this rescaled 
equation is the same as that for (A.38), and ||ay||o,/3,B < C for C = C(X,M) 
in (A.39), where we have used R < 1. This allows us to apply the local C 1 '^ 
interior and boundary estimates for the Dirichlet problem [20, Theorem 8.32, 
Corollary 8.36] to the rescaled problems in the balls B^ d /g( z b) an d B^/^zq) 
as in Lemma A.l. Then, scaling back and multiplying by d Za , applying the 
covering argument as in Lemma A.l, and recalling the definition of F, we 
obtain that, for any zq G B^ 5 U S9/5, 
(A.40) 

4 +/3 [ w l]l,/3,B dzo/16 ( 2 „)nB+ /5 + d2 z o i w i}l,0,B d ^ /16 (z o )nB+ 5 

< C7(rf ato ||J5T*|| 0j0i _ B<ejBo/a(;Bb)rijB + o H- rf^^Mi^s^^^jn^BH-^ + ll/llo,/3,B dzo/2 (^)niJ 9 + /5 
+d 2 zi l3 [g] lif 3 > B dza/2 (z )nB+ & + d t +1 b]fc,o,B £i , o/2 ( 2 o)n J B 9 + /5 )' 

k=0,l 

where we have used d Zo < 2. Recall that Dw\ = (Di±u, D\ 2 u). Expressing 
D 22 u from equation (A.l) by using (A.7)-(A.8) and (A. 36) to estimate the 
Holder norms of D 22 u, in terms of the norms of Duu, D 22 u, and Du, and by 
using (A. 18) and (A. 9) to estimate the terms involving g in (A.40), we obtain 
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^9/5 



from (A. 40) that, for every zq G B^ 5 U £2, 

d 2 z ^[D 2 u] 0t(3 ^ Bd ^ /i6izo)nB + 5 +d 2 Zo [D 2 u] 0fitBd ^ /i6{za)nB ^ 5 
< c(d z jDu\\ ciBd ^ Mza)nB ^ 

+ d abll u lll,0 ) B < ,, o/a («o)nB+ 6 + ll/llo,/3,B d2o/2 ( 2o )nB+ 5 

+e(^ +/3 [ J D 2 n] 0i/3iB ^ /2(zo)nB9+/5 + dl[D 2 U ] 0fi;Bdzo/2{zo)nB +j). 
From this estimate, the argument of Lemma A. 3 implies 
(A.41) |M|^ B + 5US9/5 < ^(ll^llt^^+^us^, + e||«||* 5/3 ^+ 5USg/5 + ||/|lo,/5,s+ 5 )- 
Thus, reducing £0 if necessary and using (A. 37), we conclude 
( A - 42 ) IMIW /5 ue 9/5 < C(X,M)(\\n\\ ^ + H/llo,^)- 

Estimate (A. 42) implies a global estimate in a smaller ball and, in particular, 

\\ u \\l,a,B+ /5 < C ( X i M )(IMIo,.B+ + ll/llo,/3,B 2 + )- N ° W We Cal1 re P eat the ar S ument > 

which leads from (A. 37) to (A. 42) with f3 replaced by a, in B^ 5 (and, in 
particular, further reducing £0 depending only on (A, M, a)) to obtain 

IMIW /5 u£s /5 ^ C(\,M,a)(\\u\\ 0tB + + ||/||o ia ,B+), 

which implies (A. 14) and hence (A. 10) for the original problem. Theorem A. 2 
is proved. 

Now we show that the estimates also hold for the Dirichlet problem. 

Theorem A. 3. Let A > and a G (0,1). Let $ G C 2 '°(R) satisfy (A.5) 
and n+ := B R f]{y > $(x)} /or i? > 0. Let u G C 2 (0[) n C(fif) safe/t/ (A.l) 
in and 

(A.43) u = g on T$ := B 2 n {y = 

where Aij = Aij(Du, x,y) and Ai = Ai(Du,x,y), i,j = 1,2, and / = f(x,y) 
satisfy (A.7)-(A.8), and g = g(x,y) satisfies 

(A.44) llslU«.(n+) < A -1 , 

with (A, a) defined above. Assume that \\u\\ctn+) < M. Then 

(A.45) IM| C2 , a(??) < C7(A,M)(|M| C(? ^ + Il/H^^ + ll^ll^^)- 

Proof. By replacing it with u — g, we can assume without loss of generality 
that g = 0. Also, by flattening the boundary as in the proof of Theorem A. 2, 
we can assume $ = 0. That is, we have reduced to the case when (A.l) holds 
in B^ and u = on £2. Thus, u x = on £2- Then estimate (A.45) follows 
from Theorem A. 2. □ 
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We now derive the estimates for the oblique derivative problem. 

Theorem A. 4. Let A > and a G (0,1). Let $ G C 2 '°(R) satisfy (A.5) 
and £1+ := B R D{y> $(x)} for R>0. Let u G C 2 (0+) n C 1 ^) 

(A. 46) Auu xx + 2Ai 2 u xy + A 2 2U yy + Am x + A 2 u y = in 
(A.47) &i« x + b 2 u y + cu = on F$ := B 2 Pi {y = $(x)}, 

where Aij = Aij(Du,x,y) and Ai = Ai(Du,x,y), i,j = 1,2, satisfy (A.7)- 
(A.8), and 6j = bi(x,y),i = 1,2, and c = c(x,y) satisfy the following oblique- 
ness condition and C 1,a -bounds: 

(A.48) b 2 (x,y) > A /or (x,y) G 

(A.49) 11(61,62,0)11^^ < A" 1 . 

Assume that \\ u \\ c ^+^ < ^f- 27ien i/iere exists C = C(A, M, a) > suc/i i/ia£ 

(A.50) W u Wcr.~(n?) - C ^c^fy 

Proof. Step 1. First, we flatten the boundary T$ by the change of coordi- 
nates (X, Y) = #(x, y) = {x,y- $(x)). Then (x, y) = ^~\X, Y) = (X,Y + 
<S>(X)). From (A.5), \\n C ^t) + H* _1 llc».«.(D+) < C( A )> where 2>+ := 
satisfies D+ c := {F > 0} and r := dV+ n {F = 0} = By a stan- 

dard calculation, F) = u(x,y) := u(^~ 1 (X, V)) satisfies the equation 
of form (A. 46) in and the oblique derivative condition of form (A.47) on 
To, where (A.7)-(A.8) and (A.48)-(A.49) are satisfied with modified constant 
A > depending only on A. Also, \v Hc(x>+) — M. Thus, (A.50) follows from 

(A.51) \H* 2 ,a,viuT ^ C(\,M,a)\\v\\ 0tV} . 

Next we note that, in order to prove (A.51), it suffices to prove that there 
exist K and C depending only on (A, M, a) such that, if v satisfies (A. 46)- 
(A.47) in B+ and Si := B x n {y = 0} respectively, (A.7)-(A.8) and (A.48)- 
(A.49) hold in Bf, and \v\ < M in Bf , then 

(A.52) II^IU-^^^I^Bf)- 

Indeed, if (A.52) is proved, then, using also the interior estimates (A. 4) in 
Theorem A.l and applying the scaling argument similar to the proof of Lemma 
A.l, we obtain that, for any zq G T>2 U £2, 

d2 zt a \\ v \\ c ^(B dzo/(16K) (z )nv+) - C W v Wc(B dzo/2 (z )nv+)- 

From this, we use the argument of the proof of Lemma A. 3 to obtain (A.51). 

Thus it remains to show (A.52). First we make a linear change of variables 
to normalize the problem so that 

(A.53) 6i(0) = 0, 6 2 (0) = 1 
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for the modified problem. Let 

(X, Y) = y) := (b 2 (0)x - h(0)y, y). 

Then 

(x,y) = y~\X,Y) = {X + b 1 (0)Y,b 2 (0)Y), |M| + {D^l < C(A), 

where the estimate follows from (A.48)-(A.49). Then the function w (X, Y) := 
v(x, y) = v(X + 6i(0)Y, b2{0)Y) is a solution of the equation of form (A. 46) in 
the domain ^(B^ ) and the boundary condition of form (A. 47) on the boundary 
part *(Si) such that (A.7)-(A.8) and (A.48)-(A.49) are satisfied with constant 
A > depending only on A, and (A. 53) holds, which can be verified by a 
straightforward calculation. Also, HHIc(>i>(.B + )) — M. 

Note that C R| := {Y > 0} and tt(Ei) = d^(B+) n {Y = 0}. 

Moreover, since + < C(A), there exists i^i = ifi(A) > such 

that, for any r > 0, C ^(-Br) C Thus it suffices to prove 

IHU^) ^ cikllccB+j 

for some r G (0, \/K\). This estimate implies (A. 52) with K = 2K±/r. 

Step 2. As a result of the reduction performed in Step 1, it suffices to prove 
the following: There exist e € (0, 1) and C depending only on (A, a, M) such 
that, if u satisfies (A. 46) and (A. 47) in B^ and on T, 2e respectively, if (A. 7)- 
(A.8) and (A.48)-(A.49) hold in B+, and if (A.53) holds and ||«|| 0iB + < M, 
then 

\\ U \\2,a,B+ < C \\ U \\o,B+- 

We now prove this claim. For e > to be chosen later, we rescale from 
B^ e into B^ by defining 

(A. 54) v(x,y) = -(u(ex,ey) - u(0, 0)) for (i,i/)eB 2 + . 

Then v satisfies 

(A.55) A u v xx + 2A 12 v xy + A 2 2Vyy + AiVx + A 2 v y = f in B 2 + , 
(A. 56) v y = b\v x + b 2 v y + ct> + cu(0, 0) on E27 

where 

Aj(p,x,y) = Aij(p,ex,ey), Ai(p,x,y) = sAi(p,ex,ey), 

h{x, y) = -h(ex, ey), b 2 (x, y) = -b 2 (ex, ey) + 1, c(x, y) = -ec(ex, ey). 

Then A\j and A\ satisfy (A.7)-(A.8) in B^ and, using (A.49), (A.53), and 
£ < 1, 

(A.57) ||(&i,62,c)|| 1)aiB + < Ce for some C = C(A). 
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Now we follow the proof of Theorem A. 2. We use the partially interior 
norms [20, Eq. 4.29] in the domain B% U Y, 2 whose distance function is d z = 
dist(z, <9i?^~ \ S2). We introduce the functions wi = Div, i = 1,2, to conclude 
from (A. 55) that w\ and w 2 are weak solutions of equations 

(A.58) £>i (-f±D lWl + -^Z^u-i) + D 22 w x = -Di {-^Dyv + ~f-D 2 v) , 

A 2 2 A 2 2 A 2 2 A 22 

2v4i9 Aoo Aa Ao 

(A.59) D n w 2 + D 2 (-f^i ^2 + -^^2^2) = -D 2 {^-D 1 v + -y^A^) 

An An An 

in £^f", respectively. From (A. 56), we have 

(A. 60) u>2 = g on £2, 

where g := &i« x + b 2 v y + cv + cu(0, 0) in B^- 

Using equation (A.59) and the Dirichlet boundary condition (A. 60) for w 2 
and following the proof of Lemma A.l, we can show the existence of (3 € (0, a] 
and C depending only on A such that, for any zq G B^ U £2, 
(A.61) 

d foK]o,/3,B £iio/16 ( 2o )n J B+ ^ ^(11^110,^/2(20)05+ + (i 2o[5]o,/3, J B t!zo/2 (2o)niJ + )- 

Next we obtain the Holder estimates of Dv if e is sufficiently small. We 
first note that, by (A. 57), g satisfies 

(A.62) \Dg\ < Ce(\D 2 v\ + \Dv\ + \v\ + ||u|| 0) b+) in 5+, 

(A.63) [g} ,f3,B d ; /2 (z)r>vl ^ C£ (\\ v Wi,f3,B d ; /2 (z)nvi) + \\ u \\o,B+ ) 

for C = C(A). The term £||w|| os + in (A.62)— (A.63) comes from the term 
cu(0, 0) in the definition of g. We follow the proof of Lemma A. 2, but we now 
use the integral form of equation (A.59) with test functions ( = i] 2 (w2—g) and 
£ = r] 2 {w2 — w 2 (z)) to get an integral estimate of |-Du>2| and thus of \Dijv | for 
i + j > 2, and then use (A. 55) to estimate the remaining derivative Duv. In 
these estimates, we use (A.61)-(A.63). We obtain that, for sufficiently small e 
depending only on A, 

(A 64) d z^ Dv ^l 3 > B «*oM z ^ nB Z 

< C {\\ v \\cHB dzo/ 2(z )nB+) + £d zo[ Dv ]o,P,B dzo/2 (z )nT>+ + ed z ll n llo,B+ ) 

for any zq G -B^ U £2, with (7 = C(A). Using (A. 64), we follow the proof of 
Lemma A. 3 to obtain 

(A.65) K^B+UE, ^ C '(IMII,o,B+uE 2 + e Mi, /3 , B2+uE2 +^lkll ,B+)- 

Now we choose sufficiently small e > depending only on A to have 

Hl i/3 ,B 2 + US 2 ^ C ( A )(IMIl,0,B 2 +UE 2 + IMIo,B+ )• 
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Then we use the interpolation inequality, similar to the proof of (A. 36), to 
have 

( A -66) IMIi,/3,b 2 + u£ 2 < C(X){\\v\\ 0B + + ||u|| 0jB +). 

By (A. 54) with e = e(A) chosen above, (A. 66) implies 

( A - 67 ) W U W*l,p,B+UBl ^ C W\\ U \\o,B+ ■ 

Then problem (A. 46)— (A. 47) can be regarded as a linear oblique derivative 
problem in B^ £ .^ whose coefficients a>ij(x,y) := Aij(Du(x,y),x,y) and ai{x,y) 

:= Ai(Du(x, y),x, y) have the estimate in C°'^(-B^ 4 ) by a constant depending 
only on (A, M) from (A. 67) and (A. 8). Moreover, we can assume (3 < a so 
that (A. 49) implies the estimates of (6«,c) in C 1,/3 (.B^ 4 ) with e = e(A). Then 
the standard estimates for linear oblique derivative problems [20, Lemma 6.29] 
imply 

(A.68) \\u\\ 2ABt/2 <C(\,M)\\u\\ 0jBie/i . 

In particular, the C°' Q (i? ; ^. //2 )-norms of the coefficients (a^,aj) of the linear 
equation (A. 46) are bounded by a constant depending only on (A, M), which 
implies 

\\ U \\2,a,B+ < C( A ; M )IMIo,B 3 + /2 > 

by applying again [20, Lemma 6.29]. This implies the assertion of Step 2, thus 
Theorem A.4. □ 
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